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Summary. — General expressions for the transition matrix are given for 
reactions between four particles of arbitrary spins. Restrictions imposed 
by 7, C or P invariance requirements are expressed in terms of simple 
equations between some phenomenological scalar coefficients which occur 
in the general expression for the transition matrix. For any given type 
of reaction one can find in this way definite physical relations which have 
the advantage of being free of any special assumptions and approxi- 
mations. 


1. — Introduction. 


Recent discovery of the non-conservation of parity (!) urged the physicists 


to a thorough revision of some other conservation laws as well. It seems to 
be certain that the invariance under charge conjugation (C) cannot be an 
exact law either. However, it has been suggested (?) that all the interactions 
may be still invariant under the product (CP) of charge conjugation and space 
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inversion. Because of the so-called 7CP theorem (*) the invariance under (CP) 
is rather completely equivalent to the invariance under time reversal transform- 
ation 7. The problem whether all the physical processes are in fact invariant 
or not invariant under time reversal is still open. Solution of this problem 
seems to be much more difficult than in the case of space inversion. There 
the mere presence of a spatial pseudoscalar in the cross-section was an unequi- 
vocal proof of the non-conservation of parity. On the other hand time re- 
versal invariance does not exclude the occurrence of time pseudoscalars in 
the cross-sections. Thus, one cannot prove violation of time reversal invariance 
without a more detailed theory. Several authors have been investigating this 
problem (*’) but most of their results have only a limited value as they are 
based on a whole series of simplyifying assumptions and approximations. 
Most authors investigate just the primary interactions and not their final 
results, i.e. the reaction itself. In order to obtain definite results they usually 
assume that all the interactions are local, do not contain any derivatives ete.; 
sometimes they even assume a particular type of coupling. Furthermore, most 
calculations are based on the perturbation expansion or other approximative 
methods and it is often not clear which results are generally valid and which 
are not. With respect to the B-transitions there are also some doubts whether 
the usual treatment of these processes is correct (5). 

It seems, therefore, advisable to investigate the problem on the basis of 
the general phenomenological S-matrix formalism which can be made free of 
most of the simplifying assumptions and other disadvantages mentioned above. 
The aim of this paper is first to construct a general theory of the transition 
matrix with the least possible number of assumptions and then to investigate 
what restrictions are to be imposed due to P, C, or 7 invariance requirements. 
In the present paper we shall discuss only the general and rather formal re- 
lations. A discussion of physical implications and some examples will be given 
in another paper. 

In order to be able to discuss not only reactions between so-called ele- 
mentary particles but also those which involve some composed particles 
(e.g. nuclei) we shall consider in this paper reactions between four particles of 
arbitrary spins. 
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Notation. — Small Greek indices (e.g. x, f, y, à) are used as ordinary four- 
dimensional tensor indices which run from 1 to 4. Small Latin indices 
(e.g. r, 8, u, v) are used as abbreviations for whole groups of tensor indices 
of some fixed but unspecified rank and order. Capital indices are reserved 
for other purpose. With respect to all tensor indices we adopt the usual sum- 
mation convention: repetition of a tensor index in an upper and a lower po- 
sition means a summation over this index. Lowering or raising of tensor 


indices is performed by means of the metric tensor g’” = Jap (Qu = 923 = 93 = 1, 

Y= +1). Thus we have for the position vector x“ and the momentum 

vector p" 

x == (r,t), Lb (—r, 1), p*=(p,®), Py (—p, w), Dye pe = == pe 
A tensor function G@““(p,«) can be, therefore, denoted by G"(p*). We 


shall have often to do with the transformation 


Yuv... 


Gaia pae) CPAS Caldi piaci ATE (p, ©) 


followed by the transformation 


o, (P; ©) =e GA (— Pp; ) Gi 


Using our conventions we can write the result of these two transformations 
in the following more compact form 


G" (p*) > G,(p, ). 

The symbol y” = (y, ya) denotes four 4x4 Dirac matrices which satisfy 
the following relations 
D TT ey tao, eil e e 
It follows that the matrix y;= éy1y»Y3y1 is Hermitian and anticommutes 
with y,: 

ea i ea Head 

The symbol (4), (*) and (—) denote the operations of Hermitian conjugation, 


complex conjugation and transposition respectively. 
Natural units are used in which c=h=1. 


2. — Definition of P, C and 7 operations. 


Since we are interested in transitions between asymptotic states in which 
only free (i.e. non-interacting) particles are present, we shall use the inter- 
action picture. In the interaction picture all the field operators satisfy free 
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field equations of motion and free field commutation relations. Bosons of any 
particular spin value can be described by a tensor field g'(x°*) which satisfies 
the Klein-Gordon equation 


(1) (0,0% + m*) g(a”) = 0, (0, = 0/00") 


and some other additional equations whose form depends on the spin and 
mass of the particles. Fermions of any particular spin value can be described 
in Rarita-Schwinger formalism by a « spinor-tensor » field y'(x*) which satis- 
fies the Dirac equation 


(2) (—iy*d, + m)p"(a°) = 0 


and some other equations whose form depends on the spin and mass of the 
particles. It is understood that besides the tensor index («) the function y" 
has also a four-spinor index which will be not denoted explicitly. Both y” 
and y" are g-numbers. 

P, C, and T operations will be now defined as follows: 


a) Space inversion operator P is a linear, unitary operator defined by 


the following transformation laws for the fields g” and yp": 


Pp'(e°)P = 9%, (L,), Po" (x) P= nv, (2,), 
(3) 
Ri Gea Wade si; AAA Pip (3). = TA CA A ’ 


where y"=y"*y,. The phase factors 7,, 7, are ordinary c-numbers of mo- 
dulus 1. 
b) Charge conjugation operator C is a linear, unitary operator defined 
by the following transformation laws 
vray Sid TONO T+ [an VAT fin SAL PINE Sed fe 
Coen Os =), Con) = 1) nag (a), 


Oy"(x*) bei (= VPage Cy'(x"), Cy"(x*) 07 L SE. > (— Ly Hep ayes 


(4) 


where the 4x4 matrix © is defined by 
(5) CHL (Ci TR Cyret=—y,- 


For the sake of simplicity we do not consider isotopic spin in this work 
The symbols 7, and 7, denote the ranks of the tensor indices (r) and (w) res- 
spectively. H.g. t,=t,=0 for 0 and } spins, 7,= 7,=1 for 1 and à spins etc. 
The phase factors 7,, 7), are ordinary c-numbers of modulus 1. 
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ce) Time reversal operator 7 is an antiunitary (i.e. anti-linear-unitary) 
operator defined by the following transformation laws (7) 


TEE) FT = nr Le), TP (a) TP" = np; (2), 
(6) 
Dig Gy = Her yp (am) , pit) Te 


| 


17 ue #,)ys€ Ki 


, . . 
The phase factors 7, and 7, are again ordinary c-numbers of modulus 1. 
As P, © are linear and 7 antilinear we have for any c-number 


(7) Pole, OC 6G-*= 6. Lol gt 


Furthermore we require that under any one of the transformations P, C and 7 
the vacuum be transformed into vacuum. 

In the following we shall use mainly momentum representation. Since all 
the (a) and w°(x°) satisfy Klein-Gordon equation we define their Fourier 
transforms as follows: 


where dp =d,pdqm; à is Dirac delta function. Inserting (8) into (3), (4) and 
(6) we find the transformation laws for the momentum representation: 


À PRD) P= 10 (D,) Po'*(p)P = np; (p*) , 

Ss Pyp"(p°) PT = Hp Pula)» Pp"(p°)P = He PP Va » 

Co GC Sey SR LE 
ie Cp) = Le), PP)! =—(Ayity"(pyer, 
ii To (p*)T* = 1,0, (Pq) ; To (PTE HP, (2,5 

si Typ) TI = hye Yad » ARIA II 


Lowering (or raising) of external (7, v) as well as of internal (x) indices at P 
and 7 transformations is to be noted. 


(7) This definition of time reversal is just a simple extension of that used e.g. by 
G. Lüpers (Ann. of Phys., 2, 1 (1957)) and others. S. WATANABE (Rev. Mod. Phys. 
27, 40 (1957)) uses another definition of time reversal which differs from 7 by the 
operation of complex conjugation. 
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Let us define charge conjugate (antiparticle) fields 
(10) PP) = pp); pop”) = Ep (—p°). 


Under proper Lorentz transformations g and y, transform exactly in the 
same way as gy” and y" respectively. Under P, C or 7 transformations q, 
and y transform essentially in the same way as g’ and y" respectively but 
the phase factors can be different in this case. 

From (3), (4) and (6) one can easily find the transformation laws under 
any product of P, C and T. For triple products 


W, = TGR, W, = PTC, W, = GER; 


(11) 
W,=PCT, W,=CPT, W,=—TPC, 
we find 
We) Wet = (1g), 
(12) Wg'* (Mo == (— 1)n5 g'(— DIG 
VEC VERI ZONA 
Weve") Wa = 1atw 0): 
with 
DST een 
de Ma = MMM » Ms = MoMoMes Ne = Mr No Nes 
Hae te Tegel ge iy et 0 
UP = Noe ’ 9; Fe ne Uetie ; n “ VEE . 


From now on we shall forget about the co-ordinate space regarding (7), (9), 
(9°), (9") as basic relations that define the P, C and 7 operations in the mo- 
mentum representation. In the following we shall ascribe quantised Boson 
(g°(p)) or Fermion (y"(p*)) fields not only to the so-called elementary part- 
icles but also to nuclei that are composed particles. This is surely not allowed 
in the Lagrangian or Hamiltonian formalism since the primary interactions 
can depend on the structure of the particles, i.e. on some internal co-ordinates. 
However, in the S-matrix formalism we are interested in the final results of 
these interactions and we have to do with free particles whose internal degrees 
of freedom do not change any more. Thus it seems possible to describe a freely 
moving (non-interacting) nucleus just in the same way as we describe free 
elementary particles. The structure of a particle manifests itself only in the 
interaction and this can be accounted for by taking the most general form of 
the transition matrix. 
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3. — Conditions imposed on the S- matrix by P, C and 7 invariance requirements. 


In this section we shall study the conditions imposed on the S-matrix by 
P, C, T and W, invariance requirements. The S-matrix is given by the 
formula 


(13) Ei 
with U(t, 0) and U(0, —t) defined as solutions of 


d 
dt 


ie NARICI 


Ud) di 


Ad = 2100, =) Het 0. 


H(t) = H(t) is the interaction Hamiltonian. 
The theory is invariant under P, C or 7 transformation if 


(15) PH(t)P-:= H(t), CH()0-:= H(), THWT-*+=H(-1) 


respectively. It follows from (13), (14) and (15) that, when applied to the 
S-matrix, the P, C or 7 invariance conditions have the forms 


(16) IPS Pate 8. OS 0-255", TSE cee 


Consequently we find the following invariance condition under triple pro- 
ducts W, 


(17) | W, SW; = S*. 


Conditions (16) and (17) are quite general but they have no practical value 
unless we know the form of the matrix S. However, we can write our inva- 
riance conditions in another equivalent form which will enable us to draw 
definite physical conclusions for any given type of reaction without any de- 
tailed knowledge of S. 

Let us consider processes of some fixed type specified by the spins, charges 
and masses of all the interacting particles 


(18) A+ B+... 2 4"4 B+... 
A process obtained from (18) by changing one or more of the particles occurring 


in the initial (final) state into antiparticles occuring in the final (initial) state 
belongs by definition to the same type as (18). If we consider processes of 
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some fixed type we can split the S-matrix into two parts 
(19) S=1+i1M +R, 


where R is defined as that part of the S-matrix which does not give any con- 
tribution to any process i +f belonging to the type under consideration: 


(20) SIE 


The part M will be called transition matrix for this type of processes. Now 
we can write our invariance conditions in terms of M instead of S: 
a) P invariance: 
(21) PMP SMI, 
b) C invariance: 
(22) CLM ‘Catt aS 


c) T invariance: 


(23) TAMTSE=RM 
d) W, invariance: 
(24) WoW == 


We shall now forget the fact that we have deduced relations (21)-(24) 
from some symmetry assumptions about the interaction Hamiltonian. We 
Shall regard the invariance conditions (21)-(24) as relations which are to have 
the same form even for theories which are non-local in the interaction or have 
no Hamiltonian. Relations (21)-(24) have the advantage of giving a possi- 
bility for drawing definite physical conclusions which are independent of any 
particular approximation method and of any particular choice of the form of 
coupling. In fact for any given type of reaction we can construct a general 
form of the transition matrix which: 


a) conserves energy and momentum, 
b) conserves electric and nucleonic charges, 
c) is invariant under homogeneous proper Lorentz transformations. 


These are just invariance properties which are still believed to be « absolute ». 

In the next three sections we shall construct general forms of M satisfying 
these «absolute » invariance conditions for all reactions between four particles 
of arbitrary spins. Consequently we shall be able to express our invariance 
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conditions in a more transparent form which can be then interpreted in terms 
of definite relations between physical observables. 

In order to avoid repetitions and detailed discussions of many particular 
cases differing only with respect to the values of charges let us first make some 
general remarks about electric and nucleonic charge conservation. For any 
given type of reaction between » particles the matrix M contains products 
of n (and only n) field operators chosen suitably from the set g, vy, gi: Vi: 
where A=1... n. Instead of w* and y we can also take charge conjugates 
of y and y, given by (10). Electric charge conservation requires that M be 
invariant under the following transformation: 


+ 


(25) p, > exp [ig,el@’, 5 Gi > exp [10 E19)", 


2 


(26) vi > expligelvi, = vi > xp lime; 


whe g,(=0,1,2,....) is the electric charge of the particle A and e is an ar- 
bitrary real constant. Similar condition but with different values of q, is to 
be imposed on M by nucleonic charge conservation. Both charge conservation 
requirements restrict the number of possible combinations of @,, wi, qu: y 
which are allowed to occur in M. However, apart from a difference in the 
constant phase factors, the results to be obtained below are independent of 
any particular combination of the fields and their Hermitian (or charge) 
adjoints. Therefore, for the sake of convenience we introduce two new sym- 
bols ®’(p*) and W"(p*) defined as follows: 


(27) Du) = gp) or Gp) = po) = pp); 
(28) Pp) = py (p*) or to) = pp’) = Cy"(— P°) : 


according to the type of reaction. We shall assume in the following that all 
the M’s contain such combinations of particle and antiparticle fields which 
satisfy both electric and nucleonic charge conservation. Of course in more 
detailed considerations one must state explicitly what fields occur in the actual 
expression for M. 


4. — Processes between four Bosons. 


In this Section we shall discuss reactions between four Bosons of arbitrary 
spins, The particles will be described by four tensor fields Co), Dae) 
D'(p!), Op’) and their Hermitian adjoints. Let us discuss first the case 
when all four Bosons are different. We find for this case the following general 
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on 


expression for M which satisfies all the « absolute » requirements listed at the 
end of the preceding section: 


4 4 
ia fo > p,) IT dp, — mi) dp} {Pr Das Piu OL” + DI OL PIP, vu wo 


A=1 AT 


TUO 


where the tensors /”" and 7 are c-numbers of the following structure 


N N! 
Jirsuv __ YPSUV | aByo TUSv 
Î => 4,67 4 Vas Ge bye 
K=1 R= 


(30) 
N N! 
ps sup ! xpyò firu 
Vi ak >) DER ! > bn € GE, a BYÒ 2 


K=1 R=1 


The coefficients 4,, 4,, b,, 0, are, in general unknown, scalar functions of 
two scalar products of the four-momenta p/ (since because of the energy- 
momentum conservation only two of the scalar products are independent). 


e is the well known antisymmetric pseudotensor satisfying the relation 
€ yo 
(31) Ex By) TEA a à 


The Gi" and Gig, are known tensors constructed from the three independent 
RA and from the metric tensor g”; they do not contain the pseudo- 
tensor e° and any complex numbers. The numbers of N and N' of linearly 
independent (*) tensors of each kind and their actual form depend of course 
on the type of reaction, 7.e. essentially on the spins and masses of the inter- 
acting particles. We have to note some important identities which follow 
from the definition of the G’s 


1114 


Yrsuvee __ firsuv suv S 
Gr + (Ce ’ Tay) = 7 Cree DI, 
9 TSUL L ya Byd £ n 
(32) Ci Can Pa dal 
MEDIE DIRE 


where 7 is the rank of the tensor (rswv). 
Let us now investigate what conditions are to be imposed on (29) due to 
each of the conditions (21)-(24). 


a) TOP invariance. — Transformation of the g’s under any of the triple 
products W, is given by (11). Apart from some constant phase factors of 


(*) This means that any other tensor Q”"° constructed from p# and gl! can be 
written as a linear superposition of N linearly independent tensors G7* with scalar 
coefficients, provided that we take both equations of motion and energy-momentum 


conservation into account. The same is true for the tensors GT. 
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modulus 1 all the ®s are transformed into their Hermitian adjoints with a 
simultaneous change of sign of the arguments pi >— pl. Introducing new 
integration variables pi ——5p#p} we can transfer this change of sign of internal 
arguments from the @’s upon the /"s. Now, in order to be able to compare 
right and left hand sides of (24) we assume that kinematically different Boson 
fields commute. On the other hand in the presence of two or more non- 
commuting Boson fields the matrix M must be symmetrized in the corres- 
ponding field operators (5). Under these two assumptions it is now always 
possible to arrange the operators occuring on both sides of (24) in the same 
order and to compare the coefficients. Remembering that under W, the Js 
transform into /*’s we find that the W, invariance condition W,MW7!= M* 
reduces to 

CODE Sp) Lp) 
(33) 

ue re) pe” 


where 0, is the product of four phase factors 7,, or 77%, corresponding to four 
fields D". It can easily be seen from (30), (32) that relation (33) will always 
be satisfied provided that we take 0, — +1. 

Just the same assumptions about commutation properties of the Boson 
fields and similar arguments with suitable new integration variables are to 
be used in the subsequent analysis of P, C and 7 invariance. For the sake 
of brevity we shall state only the main results. 


b) P invariance condition (21) reduces to 


DEN LS 


SP ot Pau) = 


ARES (pi bi) é 
(34) 
07 OM (PD, lt == = LEE ), 


where 0, is the product of four phase factors corresponding to four fields @. 
Because of (30), (31) and (32) the relations (34) are equivalent to the following 
set of equations 


(8) The assumptions about the commutation relations used in this paper are suf- 
ficient but not necessary for the validity of the TOP theorem. H.g. we can assume 
that kinematically different Fermion fields commute with each other and with all 
Boson fields and the TOP theorem will remain valid with slightly changed definitions 
of the constant phase factors 0,. It is to be noted that in the case when the com- 
mutator (or anticommutator) of two field operators is a ¢-number different from zero 
the matrix M must be symmetrized (or antisymmetrized) in the corresponding field 
operators. Compare: G. LipERS: Vertauschungsrelationen zwischen verschiedenen Feldern 
(to be published). 
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There are two possible sets of solutions of (35): 


AV if Op =n 
(36) n 
ag = by = 0 if ol 


c) C invariance condition (22) reduces to 
(37) (ele ep 


Because of (30) and (32) the relations (37) are equivalent to 


I 


(38) OG pee gal per 


d) 7 invariance condition (23) reduces to 
(39) Orl ww (Pa) = DPI) « 
Because of (30), (32) and (31) this is equivalent to 
(40) Oa De 
0, and ©, denote again the products of the four phase factors corresponding: 
to the fields ® which occur in (29). 

e) M is Hermitian if 
(41) jrs — pr 
This is equivalent to 


(42) b= 4, b=a*. 


5. — Processes between two Bosons and two Fermions. 


In this Section we shall investigate the matrix M for reactions between 
two Bosons and two Fermions of arbitrary spins. It can be shown that in 
the present case the general expression for M which satisfies all the absolute 
invariance requirements has the following form 


(43) M TE tels. Di) Il {d(pi — mi) dp,}- 


A=1 


{OD DEV 


su 1r$ » 


402 


o 


es 
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where the /”s are given by 


N 


perso CRE > (4, st Veg ja uy 9 


K=1 


N 
Pirsuv "+7 YTEUT + I 
ie CE > Gr Valba 7, o) 


K=1 


The coefficients a, An, DER bi. are, in general unknown, scalar functions de- 
pending on two scalar products of the four momenta pf. The Gare known 
covariant forms constructed from the three independent momenta pf, Dirac 
matrices y" and the metric tensor g’’; they do not contain the pseudotensor 
e or any complex numbers. We have to note three important identities 


which follow from the definition of the G’s 


Va len hb "Wa = Creo Vi) = Ge A? y J 


YTSUT 


(45) I ERE es Wal 
(— 1) {GE (— DE, y ae ra Y, IGES pi ,Ÿ By Va 


Let us now investigate what conditions are to be imposed on (43) due to 
W,, P, C or T invariance requirements. In order to be able to arrange the 
field operators occuring on both sides of (21)-(24) in the same order and then 
to compare the coefficients we must make some further assumptions about 
the commutation relations between the field operators. We assume that dif- 
ferent Fermion fields anticommute but all the Boson fields commute with all 
the Fermion fields. In the presence of two or more non-commuting Fermion 
fields the matrix M should be antisymmetrized in the corresponding field 
operators (5). 


a) TCP invariance condition (24) can be now reduced to 


| re 1)"y LT (— pe. y") ys 2 LEO DE 
(46) 


A 


| etary, fee ma = Poy). 


It can easily be seen from (44) and (45) that (46) will be always satisfied pro- 
vided that we can take 0, = +-1. 


b) P invariance condition (21) reduces to 


On val EAU) AIA "4 = Lp Pur 5), 


077: Psp Au? 14 Va = (pi V43 Ÿ ") 
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Because of (44) and (45) this is equivalent to the following set of equations 


* e A NU na At 
(48) OpAx = x) ou eur Op b, = bd ? Op db, ra br : 


There are again two possible solutions of (48): 
C= OF 0) if Oe ike 
(49) 

Opt a) if 0, = —1. 


c) C invariance condition (22) reduces to 


(50) See WER pt, Nid ae = ined ioe y") 7 


Because of (44) and (45) this is equivalent to the following set of equations. 


(51) I 
d) 7 invariance condition (23) reduces to 


(52) Ory GR rsuv (Die ETS = a ad, (Pi: 7°) ihe Va 


Because of (44) and (45) this is equivalent to the following set of equations 


(53) QG by Cr = br 
e) M is Hermitian if 

(54) Pr yeh ey 

This is equivalent to 

(55) Dam RR I 


6. — Reactions between four Fermions. 


Now we shall investigate the form of the transition matrix M for reactions 
between four Fermions of arbitrary spins. It can easily be seen that the most 
general expression for the transition matrix that satisfies all the « absolute » 
invariance requirements has now the following form: 


d 


(56) M = | dla — Pa + Ps — Ps) TT (60° — mi) dp ,}- 
A=1 


2! 17 ie IRA CER Ha Da Ww, Pg, ee ie Re LR ? 
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where the /”s have the following structure 


N 
mTSUu® __ aie Ù Yrsuv 
TR DI (EV) 
Hail 


= 
OL 
— 
Le 


N 
14 


Tp => Va Gate Ve ax —7s dx) i 


K=1 


If two or more Fermion fields that occur in the reaction do not commute the 
matrix M should be antisymmetrized in the corresponding field operators (8). 

In order to be able to write M in a more compact form we have introduced 
two sets of Dirac matrices y“ and y“ with the convention that y acts only 
on ¥,, and ¥,, or their Hermitian adjoints; on the other hand y" act on x 
and Ÿ,, or their Hermitian adjoints only. All the y“ commute with all the >”. 

The coefficients a,., ae b,, b are again, in general unknown, scalar functions 
of two independent scalar products of the four-momenta pi’. The Gs are known 
covariant forms constructed from pi, 7”, y? Vs y, and g” only; they do not 
contain any other complex numbers. The matrices y, and y, can occur in 
the G’s only in the form of a factor Ves We have to note the following three 
identities which follow from the definition of the G’s: 


Vila Grew Pay? y”, Y yaya ae Gees Die DE V1) = Ge ia, W 9 4) ; 
(58) À (—1) ppt, Vs = FE (PEs 79") 
(1) CE Cri Di; vi Ai = Vis Gay VA vara . 


Using just the same assumptions and the same procedure as in the pre- 
ceding section we can again reduce the general invariance conditions (21)-(24) 
to some relations between the /s. Because of the identities (58) we can 
express these relations in terms of some simple equations between the scalar 
coefficients a,, a,, b,, b,. The results obtained in this way are formally the 
same as those of the preceding section. Thus: 


a) TCP invariance condition (24) will be always satisfied provided that 
we take 0, =+1. 


b) P invariance condition (21) reduces to 


| a=, =0 if AE A 


| a=, =0 if dI 
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c) C invariance condition (22) reduces to 


! 


I 
(60) Cox = dx; O50, = — 0, 


d) T invariance condition (23) reduces to 


e) M is Hermitian if 


7. Summary and outlook. 


In Sect. 4-6 we have constructed the most general expression for the 
transition matrix M describing all the possible processes involving four part- 
icles of arbitrary spins. Just the same procedure can be applied to the ge- 
neral case of a reaction between n particles of arbitrary spins. We always 
find that the structure of the transition matrix M is determined first by the 
product of » field operators, second by some covariant forms denoted here by 
the symbols 7, Î. The /°s are c-numbers and can always be written as sums 
of products of unknown scalar coefficients multiplied by known covariant forms 
G7- which can be chosen to satisfy some simple identities. The word « known » 
means that the G’s can always be found without any essential difficulties for 
any given type of reaction (*!). For higher spins, however, the task of finding 
the set of basic independent G’s can be quite laborious. 

It turns out that the assumption of invariance of the theory under proper 
Lorentz transformations together with some assumptions about the com- 
mutation relations between the fields are just sufficient to secure the invariance 
under triple products TOP ete. Strictly speaking the 7CP invariance requi- 
rement imposes some restriction upon the constant phase factors 0, only, but 
does not lead to any other physical relations unless supplemented by some 
additional restrictive assumptions. On the other hand the P, C and T in- 
variance requirements or Hermiticity condition impose much more serious 
restrictions upon the general form of M. Because of the identities satisfied 
by the G’s these restrictions can be always expressed in terms of some simple 


(9) Compare M. KawaGUCHI: Prog. Theor. Phys., 12, 760 (1954), where some re- 
actions between particles of spins 0, 4 and 1 are discussed. 
(1°) Some more detailed considerations about the form of basic independent tensors 


Gr" are given in another paper: J. WERLE: Bull. Pol. Ac. Sc. (to be published). 


2406 


ZHENOMENOLOGICAL THEORY OF THE S-MATRIX AND JT, C AND P INVARIANCE - 1 585 


equations between the phenomenological scalar coefficients that occur in the 
general expression for the transition matrix. 

The problem which stands now before us in that of physical interpretation 
of our formal relations. However, it seems that this cannot be done in such 
a generality as that obtained in this paper and that we must study every 
reaction Separately. Some examples and a discussion of physical implications 
will be given in another paper. 


The author would like to express his thanks to Prof. L. INFELD and Dr. W. 
KROLIKOWSKI and to the members of the Theoretical Division of CERN, 
expecially to Dr. J. PRENTKI and Dr. R. OmxES for many helpful discussions. 


RIASSUNTO (*) 


Si danno espressioni generali per la matrice di transizione delle reazioni fra quat- 
tro particelle di spin arbitrari. Le restrizioni imposte dall’esigenza dell’invarianza di 
T, C 0 P si esprimono in termini di semplici equazioni tra alcuni coefficienti scalari 
fenomenologici che compaiono nell'espressione generale della matrice di transizione. 
Per ogni dato tipo di reazione si possono così trovare relazioni fisiche definite che 
offrono il vantaggio di essere esenti da qualsiasi ipotesi e approssimazione generali. 


(*) Traduzione a cura della Redazione. 


38 - Il Nuovo Cimento. 
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On the Symmetry Law and the Pion Decay. 


H. UMEZAWA and H. MIYAZAWA 


Department of Physics, University of Tokyo - Tokyo 


(ricevuto il 15 Marzo 1958) 


Summary. — Assuming the global symmetry of the strong pion-inter- 
actions, a simple combination of the V- and A-types the (-interaction 
is introduced to explain the problem of the B- and radiative B-decays of 
the pion. This combination is in good accordance with the rule saying that 
the weak interactions are invariant under the transformation, under which 


= 


the nucleon and the E are interchanged and y; appearing in the baryon 
changes its sign. 


It has long been known that the z-nucleon interaction together with the 
interaction for $-disintegration of the nucleon leads us to the B- and radiative 
6-decays of the pion. So far no events of these decay processes, have appeared. 
in experiments. This means that 


[ofr —e +v) or om —e +v+y)]/o(r >u+v) S105. 


Since the life time for the x-u decay is known to be 2.5-10-* s, the life times 
for the {- and radiative 6-decays of the pion are longer than 10-3 s. 

Only the P- and A-types of the 6-interaction contribute to the non-radiative 
B-decay of the pion as far as we disregard the radiative corrections due to the 
electromagnetic and weak interactions. The B-interaction of the P-type gives 
us the life time shorter than 10-58 s. As for the A-type, the lowest order ap- 
proximation in the perturbation theory leads to the life time T= 2:1045s 
with the choice of the cut-off momentum equal to the nucleon mass ('). This 
life time is shorter than the experimental value. 

As for the radiative B-decay (x + e+v-++y), the life time t, in the case of 


(7) S. B. TREIMAN and H. W. Wy tp, jr.: Phys. Rev., 101, 1552 (1956). 
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the V- or A-type for the £-interactions is quite large, t, being 3 s for the A-type 
and 0.48 for the V-type. The £-interaction of the T-type gives a somewhat 
longer life time of — 1073 s (1). 

The recent trend in the theory of 8-decay indicates the combination of 
V-A for the f-interaction. The presence of A-type interaction, however, causes 
difficulties since it gives a too short life time for x —e+y as we have seen. 
It has been known that the difficulties concerned with the 8- and radiative 
B-decays of the pion are hard to be settled by means of. a linear combination 
of the interactions (?). KAWAGUCHI and NISHIJIMA (3) discussed the situation 
by taking into account the effect due to the E-particle, and showed that the 
rate of the radiative B-decay is much reduced in the case of the 7- and V-coup- 
lings when all the interactions concerned are assumed to be invariant under 
the transformation exchanging the = and nucleon. The problem concerned 
with the A-type, however, has been left open. In this note we shall give a 
theory, where both the B- and radiative B-decays are forbidden or much re- 
duced even when the A-type exists. 

Let us discuss the pion decay processes due to the }-interaction of the 
V- and A-types: 


(1) H, uf f(yT--7, (a sn bys)p)L, ’ 


where 


tN 


nu Paule + Pro), - 


Here y is the field operator of the sixteen components describing the &-particle 
and nucleons. It is assumed that a and b can depend only on the matrix u, 
whose eigenvalue is 1 for the nucleon and —1 for the =. We shall further 
assume the global symmetry for the strong z-interactions (*). We shall first 
discuss the process 7 >e+v. A cancellation occurs between the nucleon loop 
and the © loop when b—% is assumed. In this case the transition matrix 
element is roughly proportional to (AM/M), the mass difference between the 
nucleon and =. One then obtains t, = 2 -10-4(M/AM)? = 10-38 which is not 
inconsistent with the experiment. 

In a same way the transition probability for the radiative 8-decay is also 
proportional to (AM/M)?, when b=w and a is independent of wu. It should 
be noted that the situation does not much change even when 4 depends on wu. 
This is because the transition probability of the radiative £-decay due to the 
A- and V-types is quite small whatever the constants a and b are. 


(2) See, however, F. Low and K. HuanG: Phys. Rev., 109, 1400 (1958). 

(3) S. KawacucHI and K. NisHiJima: Phys. Rev., 108, 907 (1957). 

(4) J. ScHwINGER: Phys. Rev., 104, 1164 (1956); M. GELL-MANN: Phys. Rev., 106, 
6 (1957). 
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One can prove that both the B- and radiative f-decays are completely 
forbidden, when b =, a is independent of # and the global symmetry holds 
(i.e. AM and the radiative corrections are neglected). Let us call the charge 
conjugation, T-transformation (y > 11, PP? > 9, OP >— oF, PS > — pr) 
and the exchange of the nucleon and &-particle the C, 7 and @ transformations 
respectively. In the theory of global symmetry the strong x-interactions are 
invariant under the CTG-transformation. Under this transformation 7_, 4, 


and b in (1) change as follows, 
T_T, a>—- 4, b>—D. 


Thus the process 7 >e+v is forbidden in the approximation under consi- 
deration. The current of the nucleon and the S-particle is of the form 


ME — 
‘9 AU ae T3)Y 


and one has 


Vi mac: 2 UM UN Tarde ES 


under the C7G-transformation. The current is thus invariant under this trans- 
formation. Same is true for the current of the pion. In this way one can 
prove that the process 7 >—e+y+y is also forbidden in the approximation 
under consideration. 

However, the decay probabilities are not exactly zero because of the mass 
difference AM. It should be noted, however, that AM would not be neces- 
sarily equal to the observed value; it could be smaller than the observed value 
because high energy baryons are involved in the virtual states. However, there 
is a reason to imagine that AM could not be much different from its observed 
value. Indeed, assuming the observed value for AM, the theory of global 
Symmetry gives the life time 5-10-15 s of the neutral pion (5), which is not 
much different from the maximum value of the experimental life time. 

The type of the interaction (1) suggests us a possible way of extending the 
law of global symmetry into weak interactions. This is to assume that the 
latter interactions are invariant under the transformation, under which the 
nucleon and the & are interchanged and y; appearing in the baryon parts of 
the weak interactions changes its sign. An experimental proof for this sym- 
metry law may be given by the cascade decay of the =-particle: 


HA Lx, A>ptr. 


(5) J. Tromno: Nuovo Cimento, 6, 255 (1957). 
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Indeed, the A produced in the first step should be polarized in the direction 
of its momentum, and therefore, the distribution @(0) of the pion due to the 
subsequent A-decay would show forward-backward asymmetry (°): 


D(0) c (1 +a cos Ü). 


Here 0 is the angle between the momentum of the produced A-particle and 
that of the pion in the rest system of the A-particle. The global symmetry 
law for the weak interactions says that the asymmetry factor « should be 
positive. 

Let us now come to the problem of the process 7 +—u+v. We are forced 
to introduce the direct interaction for this process as far as we assume the 
interaction (3) (with m,—,: the u-meson field) with b =u, because the latter 
interaction leads us to a too long life time (>10-5 s) for the process under 
consideration. We can omit the direct interaction when we assume the fol- 
lowing interaction: 


fpr Vul@ as V5) j (qu Val +s Bys)Py) 


In other words, we can explain both of the (xt ~ pw) decay and (x +e) decay 
by means of the weak fermion-fermion interactions, when b=w and b=1 
are assumed for the (-interaction and for the baryon-u-meson interaction, 
respectively. In this way, however, the universality of the weak interaction 
is lost. 


We are grateful to Prof. S. NAKAMURA and Professor S. MACHIDA for their 
inspiring discussions. We are also indebted to Dr. M. KONUMA, and K. NAKA- 
GAWA for their excellent cooperation. 


(6) H. Umezawa, M. Konuma and K, NAKAGAWA: Nuclear Physics, Vol. 7, No. 2. 


RIASSUNTO (*) 


Assumendo la simmetria globale delle interazioni forti dei pioni, si introduce una 
semplice combinazione dei tipi V ed A dell’interazione $, onde spiegare il problema 
dei decadimenti 8 e @-radiativi del pione. Tale combinazione è in buon accordo con 
la regola che stabilisce che le interazioni deboli sono invarianti rispetto alla trasforma- 


= 


zione che scambia nucleone e E e muta il segno del y; che compare nel barione. 


19 


(*) Traduzione a cura della Redazione. 
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Slow Mesons Accompanying High Energy Interactions (*). 


B. CHAMANY (+) and A. G. BARKOW 


Marquette University - Milwaukee, Wisconsin 


(ricevuto il 5 Aprile 1958) 


Summary. — Four nuclear events were found in G-5 emulsions, exposed in 
Texas at an altitude of 104000 feet, where the breakup of the primary or 
the resulting interaction was accompanied by the emission of slow pions. 
Two of the events could be analysed as nucleon-nucleon interactions, and 
the energies of the slow pion produced were (8.06.09) and (8.6+0.1) MeV 
respectively. On the basis of existing theories on multiple meson pro- 
duction an attempt was made to calculate the probability of observing 
the production of a 10 MeV pion or less. Computations for two energies 
of the incident nucleon were made. When the incident nucleon in the 
laboratory system of co-ordinates was 1.38 GeV, one meson was produced 
and the probability that it had an energy of 10 MeV or less was found 
to be 6 in 100. If the incident nucleon is 11.1 GeV (N=2), the proba- 
bility that one meson had an energy of 10 MeV or less was found to be 
2 in 1000. 


1. — Introduction. 


In the systematic scanning of a stack of G-5 [ford emulsions, four high 
energy interactions accompanied by slow mesons were observed. The stack 
was exposed at San Angelos, Texas, at a magnetic latitude of 31° N for a 
period of better than 10 hours at an altitude of 104000 feet. In a total of 
40000 stars, 37 events were observed involving heavy nuclear collisions or 
breakups. The charged distribution of 37 events is given in Table I. The 
energy of the fragments was obtained from grain count and scattering measu- 


(*) Work supported by the National Science Foundation. 
(~) Present address: Purdue University, West LaFayette, Indiana. 


2412 


ile) 
[nei 


SLOW MESONS ACCOMPANYING HIGH ENERGY INTERACTIONS 5 


rements while the energy of the heavy primary was obtained by the opening 
angle measurement of the «-particles (1). The identification of charge was 
made by 3-ray count and gap count. 


TABLE I. — Charge-energy distribution of the 37 events. 
Z | Number | Energy (GeV/nucleon) 
— - - — — | 
2 12 1.2 — 1.6 
6 + | 0.8 — 2.2 
7 2 0.8 + 2.0 
8 2 1.1 — 1.9 
10 4 1.4 — 2.6 
12 3 Leo | 
12 10 10218 | 


The four interactions that were acempanied by the emission of slow mesons 
are pictured in Fig. 1-4. 


2. — Description of events and measurements. 


Event 1, Fig. 1. A heavy cosmic ray particle identified as a nitrogen nucleus 
enters the stack at an angle of 29° with the vertical. The nitrogen nucleus 
traveled a distance of 6.5 cm in the emulsion and then broke up into three 
«-particles and a proton. One of the «-particles traveled an added distance 
of 2.2 cm and made a nuclear interaction, while the other particles and proton 
left the plate. Their ranges in the emulsion were approximately 5 em. The 
nuclear breakup is accompanied by the emission of a slow positive pion. The 
pion is produced at an angle of 60° with the initial direction of the primary 
particles. The pion had a total range of 1576 um in the emulsion at the end 
of which it underwent a u-e decay in the plane of the emulsion. The energy 
of the pion was (8.06 + .09) MeV.. The basic reactions that participated in 
the production of the z+ can be written as 


P+P-s>P+n+7n* 


with the first one preferred for the balance of charge in the event. 


(1) M. F. Kapcon, B. Peters, G. T. ReyNoLps and D: M. Ritson: Phys. Rev. 
) 
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Event 2, Fig. 2. A primary cosmic ray which was identified as that of 
a carbon nucleus enters the stack at an angle of 30° with the zenith and after 
traveling a distance of 8 em broke up into two x tracks and three protons. 
All of the particles left the emulsion without making any further nuclear inter- 
actions. The interaction was accompanied by the emission of a low energy 
negative pion at angle of 59° to the initial direction of the primary particle 
It traveled a distance of 1765 um, after which it was absorbed by a nucleus 
of the emulsion producing a one prong o-star. The energy of the pion was 
(8.6-- 0.1) MeV. The basic reactions that participated in the production of 
the x- can be written as 


Pene 


LTT 


Ag! 
12 


n nun TTT 


with the latter reaction preferred for the balance of charge in the event. 


Event 3, Fig. 3. This event was produced by a carbon nucleus and re- 
sulted in a jet of 17 minimum tracks, two slow protons, and a slow x meson. 
The slow rt was emitted at an angle of 63° to the original direction of the 
primary, has a range of 6700 um which was followed by the characteristic 
u-e decay. The energy of the pion was (18.3 + 0.2) MeV. 


Event 4, Fig. 4. A high energy nuclear interaction of a primary «-particle 
with one of the nuclei in the emulsion. The interaction was accompanied by 
a cone of lightly ionizing tracks, two gray tracks and three pions, two negative 
ones which are identified by the o-stars they produced, and a positive pion 
which undergoes a 7-u-e decay. In the cone of lightly ionizing particles are 
two electron pairs produced by neutral 7-meson. 


The first two events were observed under very ideal conditions for measu- 
rements, being practically in the plane of the emulsion. 

The identification of charge of the primary particle producing events 1 
2 and 3 was done by è-ray counting on the primaries and also on the «-particles 
emitted in the breakups. Let Z,, Z, and Z; be the charge of the particle pro- 
ducing event (1), (2) and (3) and N,, N, and N; the number of 3-rays per 
100 um along their respective tracks, and N, the d-ray density for an «-particle. 
The Ô-ray density has been shown to be proportional to the square of the 
charge and is constant at relativistic energies The number of 8’s associated 
with the particles in question will be given by 


N,= KZ:+A where ¢==-1, 2,306 


The values of N; for the three particles are N,= 6.3 +32, N, — 4.8 + .62, 
N;= 4.46 + .49 and N,=0.65+.05. The error due to background is A. 
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Fig. 1. — Breakup of a nitrogen nucleus into three alphas, one proton and one x* 
meson of energy (8.06 + .09) MeV. The suggested reaction is iN + p> 3e+ptxn+nt 
and can be considered as a nucleon-nucleon interaction. 
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Fig. 2.— Breakup of a carbon nucleus into three protons and two alphas with a 77 
meson of energy (8.6 + 0.1) MeV. The suggested reaction is $C+n > 2x+3p+xn+7x- 
and can be considered as a nucleon-nucleon, interaction. 
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+ .5) MeV together with a jet of high energy pions and 


Fig. 3. — Collision of a carbon nucleus with a heavy nucleus and the emission of a slow 
5 


m+ meson of energy (18.34 
nucleons. This reaction cannot be considered as a simple nucleon-nucleon interaction. 
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Fig. 4. — A collision of an «-particle with a heavy nucleus resulting in the emission of 

two slow x and one slow x*. The associated jet consisted of seven high energy mesons 
a + a ss > 
four high energy protons, and two sets of electron pairs originating from two 7°. 
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If one takes the ratio of successive differences, we obtain the expressions 


GANN 
Zi—4 N,—N, 
LSA NIN, 


€ 


DATES TI N 
E ca 1360.02 


= 62 lf 0.2 


— 
bo 
— 
= 


Assuming various values of Z,, Z, and Z, the best fit is found four Z,=7 
al Ah Aas Ce 

The complete data for events No. 1 and 2 are given in Table II. 

Energy and momenta in the rest system of the primary are indicated by 
‘the asterisks. The values of f and y were obtained from the relationship 


where P;, and U are the components of momentum and total energy of the 
xs involved in the breakup. 


3. — Theory. Production of mesons. 


HEISENBERG (°); FERMI (*), LANDAU (5), and others (9-19) have designed 
various theories to explain the multiple emission of pions in nucleon-nucleon 
collision. When nucleon-nucleon encounters take place at very high energies, 
7.e., 100 GeV or more, many mesons are produced. Up to the present time 
three somewhat different theories have been proposed by HEISENBERG, FERMI 
and LANDAU. Of these the thermodynamic theory of Fermi and the hydro- 
dynamic theory of Landau appear to have the greatest potentialities. The 
theory proposed by FERMI has certain limitations, but up to the present time 
has met with some experimental verification. 

When two fast nucleons collide the energy is released in a very small volume 
F given by the common volume of two very thin Lorentz contracted discs. 


E. FERMI: Prog. Theor. Phys., 5, 570 (1950); Phys. Rev., 81, 683 (1951). 
L. D. LaAnpavu: Izv. Akad. Nauk USSR, Ser. Piz., 17, 51 (1953). 


H. Fukupa and G. TAKEDA: Prog. Theor. Phys., 5, 957 (1950). 

H. Umezawa, Y. TakaHasHI and S. KAMEFUCHI: Phys. Rev., 85, 505 (1952). 
E. CLEMENTEL and N. DaLLAPORTA: Nuovo Cimento, 5, 234 (1948); 5, 298 (1948). 
R. GLAUBER: Phys. Rev., 84, 395 (1951). 
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Where the discs collide we have a region of very high energy density and 
extremely strong interaction and, therefore, one can assume that a large pro- 
portion of the proton energy goes into the mesonic field. The volume V in 
which the energy is released is determined by the dimensions of the nucleon 
cloud, whose radius is of the order h/mc, where m is the mass of the 
m-meson (ApAv ~ h). Since the nucleons move at high velocity the cloud will 
undergo Lorentz contraction in the direction of the nucleon’s motion. Hence, 
the volume will be of the order 


ATA dre 2G 
(1) VaVey =) 


3 \me 
where JZ is the nucleon mass and £ is the total energy of two colliding nucleons 
in the center of the mass system. The particles are created in the volume V 
at the initial moment of collision according to the laws of statistical equi- 
librium, and assuming that they no longer interact with each other escape 
from the volume in a frozen state. The states involved in the equilibrium are 
the ordinary free particle states of the mesons and nucleons not the states of 
free particles confined to the spheroid. In this way, dropping constant factors 
we have 


and since temperature is energy per particle 


T~ EIN 
so that 


He Be 


AE 


If £ is the available energy in the center of mass system, we can make the 
further assumption (Eq. (1)) 


or N~ E:. 


— 
bo 
— 
— 
N 
| 


LEWIS, OPPENHEIMER and WOUTHUYSEN (5), and others have modified 
Eq. (2) by introducing the existing field theory with the result 


(3) N & 0.9 Et (GeV), 


where £, is the energy of the primary in the laboratory system. 
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Numerous attempts have been made to crystallize existing theories into 
some workable form. The two best review papers are those by MILBURN (77) 
and LEWIS (1°). 


4. — Computations and results. 


If we assume a nucleon-nucleon collision to be responsible for the events 
No. 1 and 2, we can apply the theory of multiple meson production and try 
to calculate the probability of observing a meson of some energy Æ in the 
laboratory system. Since both observed mesons were about 9 MeV we asked 
ourselves the question, « What is the probability of observing a pion of energy 
10 MeV or less in the laboratory system when the incident nucleon energy 
is sufficient to produce one pion and two pions? » Calculations for any other 
energy pion and for any incident nucleon can be made in the same way. 

In transforming an event to the center of the mass system of the two 
nucleons it is more applicable to transform total energy U and momenta P,, 
P, and P, than velocities. By a systematic application of the Lorentz trans- 
formation it can be shown that 


U* = y(U— BeP,) 
cP* = y(cP,— PU) 


a 


DEEP: 
Dio 


U* and P* are the total energy and momentum of the particle in units 
of MeV and MeV/c respectively in the moving system, U and P being the 


TABLE III. — Data for transformation of 10 MeV pion to the laboratory system. 
(For 10 MeV pion, pe = 53.75 MeV). 


No. Ey | | E* aga 2 eee oF Proba- 
of pions|(Gev)| fe | % | 8 (om), RO Biro See ala ites 
and 6 and 7b 
Eu | i Wee TE RE 
1 | 1.38) .652 | 1.32 | 98 MeV | 600 MeV | .2775 units | .0168 | 6/100 
2 | 11.1 | .923 | 2.65 | 138 Mev | 3 MeV | 2.28 units | .00384 | 2/1000 | 
| | 


U= total energy of the particle to be transformed. | 


(7) R. H. MILBURN: Rev. Mod. Phys., 27, 1 (1955). 
(12) H. W. Lewis: Rev. Mod. Phys., 24, 241 (1952). 
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respective quantities in the laboratory system. Where y is the ratio of the total 
energy of the collision to the rest energy of the center of mass. 
On the basis of a nucleon-nucleon collision LEWIS (6) has developed an 
expression 
N(M + k) M—k dp 


(4) P(0, e, p) = exp | kE, E—7 


where P(0, e, p) is the probability of finding a pion of energy e and e+de, 
momentum p and p+dp at an angle 0 with the initial direction of motion of 
the parent nucleon. If we integrate Eq. (4) over all angles, we get 
NES): V(M — k)p 
(5) P(e, de) = exp | — - Li ‘sinh SALE M ‘de 
> 0 


where P(e, de) is the probability of finding a pion between e and e+de (or p 
and p+dp) anywhere in space. £, is the energy of the incident nucleon in 
the center of the mass system of co-ordinates, N is the multiplicity, M the 
mass of the nucleon and k the mass 


of the pion in electronic masses. ZE 
Fig. 5 and 6 are plots of Eq. le 
(5) for ineident nucleon energy of 12 
1.38 and 11.1 GeV respectively. 19 — 
To determine the probability x 
of observing a pion of 10 MeV or 0 
less emitted in any direction in  :02 q 
the laboratory system of co-ordi- ceil, ©) 1 5 5 7 
nates, we transform the pion into P(e, de)=e-1%*€ sinh (1.33 p) de; N=1, Eo(LS)-1.38 GeV 
the center of the mass system of 0925300 Mevinueleon; EolM_0=245 
Area=.278 units; Area=0.25 Z (ordinates) 


the colliding nucleon. This was 
done for both the 1.38 and Fig. 5. — Energy distribution of all pions formed 


È in a 1.38 GeV nucleon-nucleon collision. 
11.1 GeV nucleons and the results i 


are shown in Fig. 7. 

All pions of 10 MeV or less in the laboratory system will have their mo- 
mentum vector in the center of the mass system terminate in the area shown. 
If we graphically integrate these areas with the aid of Eq. (4) and compare 
them with the areas of Fig. 5 and 6, we can calculate the probability of ob- 
serving a pion of any energy produced in a nucleon-nucleon collision. It must 
be remembered that Fig. 5 and 6 give the probability of finding a meson of 
total energy e anywhere in space regardless of direction in the center of 
mass system. In our problem the meson occupies a volume given by particular 
limits of angle and momentum in the center of mass space. For this reason we 
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tone d Ave iihe 
30 Pies de) des ginh (521 p) de cannot use Fig. 5 and 6 directly 
w N=2, E(LS)=11.1 GeV in finding the probability of ob- 
D + ~ 
> E,(cM)=1.5 GeV/nucleon : Roe 
DES SS | CE ig 
5 B4|Mnct=10.95 serving a meson of total energy 
20 Area=2.28 units between the values e and e de. 
sy rdinates Eten 

eae Cp rien À The probabilities were found to: 
be 6 in 100 for nucleons of 1.38 GeV 

10 and 2 in 1000 for 11.1 GeV. 


In a total of 40000 stars 
there were observed 42 events 
e(E/M, ©) that involved heavy nuclei, en- 
12345678 9 01 121314 % % 17° dings of which 4 were thindowns, 
1 was thindown associated with 


L 


Fig. 6. — Energy distribution of all pions formed 
in an 11.1 GeV nucleon-nucleon collision. breakup (1%), 4 were breakups as- 
sociated with the emission of slow 
mesons and 33 were ordinary Draps, . {The distribution of the events was 
given in Table I. te 
It must be realized that the statistical data are not sufficient to warrant the 
drawing of definite conclusions. If we are permitted to assume that Events No. 1 
and 2 can be considered as a nucleon-nucleon interaction, the probability that. 
a meson of energy of 10 MeV or less be formed is found to be 2 in 37, or about. 


cM LS 
u3300MeV/nucleon  Ye:132 Bc=0.65 E;31.38 GeV W:1 


E, 


ci 


R°- 6725 Mevl = p*-2686Meyk 6°: ts: A R:4656 x p 7 MeV/c 
CM LS 
E,4=15 GeV/nucleon ¥ 22.65 Bc=.923 Est Gev N:2 


P = 53.75 Mevk 
26° R=46.55 MeV  P:2688 MeV 


-\ 26.88 
244 


R*=244 MeV/c R*= 2688 MeV/c 8% tg 


Fig. 7. — Transformation of all x-mesons of energy 10 MeV or less from momentum 

space in the laboratory system of co-ordinates (L.S.) to momentum space in the center of 

mass system of co-ordinates (C.M.). The transformation was made for two energies 

of the incident nucleon. The top diagram for H,=1,38 GeV, multiplicity of 1; the 
bottom diagram for E,=11.1 GeV and multiplicity of 2 


(13) A. G. Barkow, B. CHAMANY and R. E. MCDANIEL: Nuovo Cimento, 6, 125 
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was 


9 in 100, which agrees favorably with theory. The study is being continued 
and it is hoped that added events will contribute much to the data. 

To date two events have been found in a stack of plates exposed for 8 hours 
at 110000 feet at the magnetic equator (Guam). Preliminary results indicate 
nuclei of Z=7 and Z=8, with plus pions being formed at angles 55° and 61° 
with respective energies of (7.06 + 0.3) MeV and (9.3 + 0.2) MeV. 


The authors are very grateful to Professor MARCEL SCHEIN of the Uni- 
versity of Chicago for his generous assistance in regard to this problem. 


RIASSUNTO (*) 


= 


In emulsioni G-5 esposte nel Texas ad un’altezza di 104000 piedi furono trovati 
4 eventi nueleari in cui la frammentazione del primario o l’interazione risultante erano 
accompagnate dall'emissione di pioni lenti. Due degli eventi hanno potuto essere iden- 
tificati come interazioni nucleone-nucleone e le energie dei pioni lenti prodotti risul- 
tarono di (8.06 + 0.9) e (8.6 + 0.1) MeV, rispettivamente. Sulla scorta delle teorie 
esistenti sulla produzione multipla dei mesoni, si è tentato di calcolare la probabilità 
di osservare la produzione di un pione di 10 MeV o meno. I calcoli sono stati eseguiti 
per due energie del nucleone incidente. Se il nucleone incidente era di 1.38 GeV nel 
sistema del laboratorio sì produceva un mesone e la probabilità che questo avesse 
energia di 10 MeV o meno si trovò essere del 6%. Se il nucleone incidente è di 11.1 GeV 
(N= 2), la probabilità che uno dei mesoni prodotti abbia energia 10 MeV o meno fu. 
trovata del 0.2%. 


(*) Traduzione a cura della Redazione. 
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M1 Transitions in ‘Pr and !Ir. 


M. C. JosHI, B. N. SUBBA RAo and B. V. THOSAR 


Tata Institute of Fundamental Research, Bombay, India 


(ricevuto l’11 Aprile 1958) 


Summary. — /‘-conversion coefficients of the 145 keV transition in 14!Pr 
and the 129 keV transition in 19Ir have been determined to be 0.34+.035 
and 2.1+.3 respectively. These values are found to be in better agree- 
ment with the finite nuclear size corrected values rather than with the 
point nucleus values. 


1. — Introduction. 


Following the theoretical investigations of SLIV and LISTENGARTEN (+) re- 
garding the finite nuclear size effect on K-conversion coeffcients, a few expe- 
rimental investigations (*") have been made. The results of these investi- 
gations indicate that for M1 transitions, the experimental A-conversion co- 
efficients are in better agreement with those calculated by SLIV and BAND (8), 
rather than with those calculated by Rose et al. (2). More recent theoretical 
work (*!°) indicates that the conversion coefficients may depend on the finer 


(1) L. SLIv and M. A. LISTENGARTEN: Zu. Eksp. Theor. Fiz., 22, 29 (1952). 
(*) M. E. Rose et al.: Beta- and Gamma-ray Spectroscopy, Ed. K. SIEGBAHN (1955) 
and ORNL-1023. 


(*) A. H. Waspsrra and G. J. Nisan: Nucl. Phys., 1, 245 (1956). 

(4) C. NORDLING et al. Nucl. Phys., 1, 326 (1956). 

(9) F. K. McGowan and P. H. StELSON: Phys. Rev., 108, 1133 (1956). 
(5) K. O. NIELSON et al.: Nucl. Phys., 2, 476 (1957). 

(7) G. T. Ewan et al.: Phys. Rev., 108, 1308 (1957). 

(8) L. A. SLIV and I. M. Banp: AECO/Tr 2888. 

(*) E. L. CHURCH and J. WENESER: Phys. Rev., 104, 1328 (1956). 

(10) S. G. Nizssox: UCRL-3803 
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features of the nuclear structure. So, it was felt that more experimental in- 
vestigations in this field would be helpful. 

In the case of the 145 keV transition in “*Pr and the 129 keV transition 
in !9!Ir, there was considerable disagreement about the conversion coefficients 
and conversion ratios among the previous authors. So, we have re-investigated 
these transitions. 


2. — Experimental methods and results. 


By studying the internal conversion spectrum and the photoelectron spectrum 
of y-rays in the intermediate image B-ray spectrometer, we have determined 
the values of conversion coefficients and conversion ratios. The spectrometer 
has been used with 8°% transmission and 4°% resolution. In one method, the 
area under the conversion peak is compared with the area under the {}-con- 
tinuum as computed from the experimentally determined Kurie plot. When 
the decay scheme is simple and well-established this method yields satisfactory 
values of conversion coefficients. In the other method, described in detail 
elsewhere (15), from the areas under the photo-electron peaks from a radiator, 
the intensities of y-rays and X-rays are compared. 


99 Ls 141pr, 


The 145 keV transition in 1*4Pr arises due to the B--decay of 141Ce (T4= 
= 33 days) to the 145 keV excited state of !4!Pr. Though the decay scheme 
seems to be well established, the conversion coefficients determined by various 
workers (7%) indicate considerable disagreement (Table IT). 

life was obtained by irradiating a « spec-pure » sample of Ce,O, in the 
Harwell pile. 1#Ce (7 4=33 hours) which decays into !4*Pr which further 
decays with a half life of 13.7 days to stable 14°Nd, was also found with 1#Ce 
These short-lived activities were allowed to decay for about four months. 
Then, using a source of about 0.2 mg/cm? surface density, deposited on an 
aluminum backing (0.15 mg/em?), the internal conversion spectrum, shown 
in Fig. 1, was observed. The counter window had a cut-off energy at about 


(11) M. S. FREEDMAN et al.: Phys. Rev., 79, 897 (1950). 
(12) R. L. HEATH: Phys. Rev., 87, 1132 (1952). 

(13) S. A. E. Jonansson: Arkiv f. Fysik, 3, 536 (1952). 
(14) E. KONDAIAH: Arkiv f. Fysik, 4, 122 (1952). 

(15) J. T. Jones et al.: Phys. Rev., 97, 1031 (1955). 

(16) G. B. ZorzoLi: Nuovo Cimento, 5, 289 (1957). 


a 39 - Il Nuovo Cimento. 
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4keV. The results of the analysis of the f-spectrum by the Fermi-Kurie plot 
(Fig. 2) are summarized in Table I. These results are in agreement with those 


of earlier workers (1). 


+>} 


36 aCe 


( We’ —- 
PT 


Oi > 120 600 1200 80 —> 1900 2400 


Fig. 1. — Internal conversion spectrum of 14!Ce. 


— 3000 3720» 


Fig. 1 shows that A and L+M conversion lines of the 145 keV transition, 
along with a weak broad line at about 32 KeV which can be identified as the 
Auger lines of !*!Pr following the internal conversion of the 145 keV transition. 
Using the intensities of the partial 6-groups from Table I and the intensities 


TABLE I. — Analysis of B-spectra of 14!Ce. 
Energy of the 8 groups | Relative intensity 1 
fet Se | : og ft 
(keV) (%) 
8, 580 +5 | 22 + 2 
B, 440 +9 78 +2 


MI TRANSITIONS IN 141Pr AND 191Ir 603 


of conversions lines, the A-conversion coefficient, x*, was found to be 0.36.03. 
The K/(L+M) ratio was found to be 5.6 + 0.4 from which K/L was cal- 
culated, in an approximation (17), to be 6.6 + 0.8. 


(aaa 


à D nre lesene 


di 120 LOTIR VIE 060 780 — 607 


Fig. 2. — Kurie plot for 14Ce. 


F* The photoelectron spectrum of the X-rays and y-rays of 14!Pr was studied. 
For this purpose, 6-rays and internal conversion electrons were absorbed in 
350 mg/em? thick perspex absorber. In view of the low energy of Pr X-rays, 
a copper converter of about 0.25 mg/cm? thickness deposited by vacuum 
evaporation on to a perspex disc, was used. The photoelectron spectrum ob- 
served in this way, is shown in Fig. 3. In this figure X,, and X,, represent 
the Æ and L photoelectrons of K, X-rays respectively, X,, and y, represent 
the K-photoelectrons of K,—X-rays and the 145 keV y-rays, respectively. 
Adopting a method similar to our earlier one (?*), the K-conversion coefficient 


(17) D. E. ALBURGER: Handbuch der Physik, 42, 26 (1957). 
| (18) M. C. JosHI, B. N. SuBBA Rao and B. V. THosar: Proc. Ind. Acad. Sci., 45, 
| 390 (1957). 
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was calculated by comparing the areas under the photoelectron peaks of X-rays 
and y-rays after correcting for the variation in photoelectric cross-section with 
energy (1°), for the ab- 

do Ya sorption of X- and y-rays 
in the source material 
and perspex absorber, 
and for the fluorescence 
yield (2) in Pr (w,=0.88). 
Thus af was found to be 
0.32 + 0.04. Average of 
a” from these two me- 
thods gave 0.34 + .035. 
Our experimental va- 
lues are compared with 
those of earlier workers 
and with the theoretical 
ones "inc Lables Live) ihe 
theoretical values of ROSE 
were obtained from the 
Ki UL, TT, ndo Con- 
version coefficients from 
the tables of ROSE et al. (?) 
by linear interpolation 
on log « against + plot. 
From Table II, by 
comparing the experim- 
ental values of K/L with 


680 1020 1360 1700 2040 2380 2720 3060 


Fig. 3. — Photo-electron spectrum of 14Ce. 


TABLE II. — 145 keV transition in 14'Pr. 


| References Theoretical 

| | = 2 

| Quan- | Present R t al.(2)|8 Aye 

| tity | | | oil OSE et al.(?) SLIV et al. (8) 
(i) | (5) | (2) (2) (25) (23) | E 

| | | | MO STON MII 

| | | Tape ic Ne FPS 

Pei” 6A ET) — 6.5 46.35 62 066220 ST PRISE =! — 

| | | | : | | 

Nata | 0.25 | 0.46 | 0.48 | 0.33 | 0.22 | 0.37 | 0.344.035 | 0.43 | 0.42 | 0.39 | 0.37 | 


(19) C. M. Davisson: Beta and Gamma-ray Spectroscopy, Ed. K. SIEGBAHN (1955), 
p. 869. 


(2°) I. BeRGSTROM: Beta and Gamma-ray Spectroscopy, Ed. K. SIEGBAHN, (1955), p. 630. 
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the theoretical ones of Ross, it is seen that this transition should be predomi- 
nantly M1 (93+ 10%) with a small admixture of H2. AMBLER et al. (21) 
from their experiments with aligned nuclei have shown that the transition 
is mostly M1 with only about 4% #2 component. From similar experiments 
CACHO et al. (22) have found that E2/M1= 0.007. HEYDENBERG and TEM- 
MER (?°) did not observe any 145 keV y-ray in Coulomb excitation experiments 
on !4!Pr thereby indicating the almost complete absence of the Æ2 component. 
Life-time determinations led H. DE WAARD et al. (24) to conclude that the 
145 keV transition in '1!Pr is a l-forbidden M1 transition. From all this evi- 
dence, it may be inferred that the #2 admixture to this transition is negligible 
so that this should be an almost pure M1 transition. 

Our observed value of «* may be compared with theoretical values. It is 
clearly seen from Table II, that the observed value is in better agreement 
with the value of SLIV et al. rather than with that of ROSE et al. The observed 
values of «, of 0.34 is slightly less than the value for a pure M1 transition 
given by SLIV et al. (*). It is possible that the nuclear model dependence pre- 
dicted by CHURCH and WENESER (°) is operative. 


4, — Tp 


It is known that !*0s decays by B--emission with a half-life of 15 days 
to the 170 keV isomeric level (life-time = 5 s) of 19Ir. This level decays by 
the emission of a 42 keV transition followed by a 129 keV transition leading 
to the ground state. The decay scheme (*°) seems to be established. Still, 
the conversion coefficients, and conversion ratios of the 129 keV transition, 
determined by previous workers (°°*’), indicate wide differences. 

Osmium powder was irradiated in the Harwell pile and the resulting acti- 
vity consisted of !850s (74=97 days), 1*10s (T3=15 days) and 1*Os (74= 33 h). 
1930s was allowed to decay sufficiently before the present observations were 
started. Activity of the electron capture isotope, 1*°Os, was very weak. 


(21) E. AMBLER et al.: Phys. Rev., 101, 196 (1955); 97, 1212 (1955). 
(22) C. F. M. Cacuo et al.: Phil. Mag., 46, 1287 (1955). 
(23) N. P. HEYDENBERG and G. M. TEMMER: Phys. Rev., 100, 150 (1955). 
(24) H. De WaarD and T. R. GerHoLM: Nucl. Phys., 1, 281 (1956). 
(25) E. KONDAIAH: Arkiv f. Fysik, 3, 47 (1951). 
(28) J. B. Swan and R. D. Hix: Phys. Rev., 88, 831 (1952). 
(2?) R. D. Hitt and J. W. MIBELICH: Phys. Rev., 89, 323 (1953). 
(28) (a) F. K. McGowan: Phys. Rev., 93, 163 (1954); (6) F. K. McGowan and P. H. 
STELSON: Phys. Rev., 107, 1674 (1956). 
(29) E. M. BernsTEIN and H. W. Lewis: Phys. Rev., 105, 1524 (1957). 
(39) J. W. MIinELICH et al.: Phys. Rev., 96, 1450 (1954). 
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Using a source of about 0.25 mg/cm? thickness deposited on a plastic 
film (0.6 mg/cm?) coated with a thin coating of aquadag to make it con- 
ducting, the internal conversion spectrum (shown in Fig. 4) was observed. 
In this figure, L, and M, re- 
present the L and M con- 
version lines of the 41 keV 
transition, A, and L, repre- 
sent the A and L conversion 
lines of the 129 keV transi- 
tion. Fermi-Kurie plot was 
constructed with the points 
between the X and L con- 
version lines of the 129 keV 
transition. The maximum 


80 
72 
64 
56 


48 


energy was found to be 
(135 + 5) keV. This corre- 
sponds to a log ft of 5.2, in- 
dicating that this transition 
is of the allowed type. 
Using the Fermi-Kurie plot 
to reconstruct the 6-spec- 
trum, relative intensities of 
the f-spectrum and the 
conversion lines were de- 
termined. The L/M ratio 
of the 41 keV transition was 
200 200 600 600 1000 1200 1400 found to be 2.6 + 0.1. The 
K/(L + M) and K/L ratios 
of the 129 keV transitions 
were found to be 3.7-+0.15 
and 4.2 + 0.2, respectively. Using the intensity of the {-spectrum as a 
measure of the total number of disintegrations leading to the 129 keV level, 
a* was found to be 2.10 + 0.25. All the values quoted above are the mean 
values of two runs. 


40 


32 


Fig. 4. — Internal conversion spectrum of 1910s. 


Using a thin tin converter (0.15 mg/cm?) evaporated on to a perspex 
disc, the photoelectron spectrum was observed. B-rays and conversion electrons 
from the source were absorbed in perspex dises (150 mg/em?). The photo- 
electron spectrum, shown in Fig. 5, indicates the K photoelectron peaks of 
the A, X-rays, Kg X-rays and the 129keV y-ray denoted respectively by 
Xagr Xx) Y, and the L photoelectron peaks of the three radiations denoted 
by Xa,; Xg,, y,, respectively. Adopting the procedure described in the case 
of ‘Pr, the K conversion coefficient of the 129 keV transition in 191Ir was found 
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to be «,=2.1+0.35, as the mean of two runs. In the evaluation of this 
quantity, the contribution of Rhenium X-rays due to the capture isotope 
Qs to the Iridium X-rays 
due to internal conversion, was 
estimated. For this purpose 
when !90Os had spent nearly 
15 half lives, the 18508 radia- 7 
tions were studied under similar 
conditions and the required 
intensity computed. This con- 6 
tribution, however, was found 
to be very small. 

The mean value from these 
two methods is «“ =2.10-L0.30. 
Our observed values. of x" and 
K/L are compared with those 
of earlier authors and with the 
theoretical values in Table III. 3 

It is seen from Table III 
that the K/L ratio from the 
present work is higher than 2 
that reported by some authors. 
Absorption in the counter win- 
dow can affect these low energy 
conversion lines. We have used 
a window with a cut-off à 
energy at about 4keV, so " 200 
that the relative intensity Fig. 5. — Photo-electron spectrum of 1910s. 


600 900 1200 1500 


TABLE III. — Conversion coefficients and conversion ratios of the 129 keV transition in 1*"Ir, 


References Theoretical 
| Quan- Present. | Roseet. al. (2)| SLIV et al. (8) 
tity = E x n° work | 
ar OM AT (3) (89) : 
| | DAME? MINE 
ole ME PAC SÒ Pe ELIZA O a a 
{ pr | 
ont == 11.36 2 |2.07+0.1 (a)! — |2.1+0.3| 2.95 | 0.46 | 2.45 | 0.44 
2.90 £0.58 (b) 
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of these conversion lines should not be affected to any appreciable extent. 

The K/L ratio or «* alone compared with theoretical values indicates that 
the transition has a strong M1 component. In order to get a better idea of 
the transition, we may compare the K/L and «* values with theoretical values 
graphically, following the method of WAPSTRA et al. (?). 

Let Bi(i=XK, L,, Dy, Lin) and of ((=K, L,, Li, Lan) denote theoretical con- 
version coefficients of M1 and H2 transitions, respectively. Let a represent 
the fraction of the M1 component in a 
M1 + E2 transition. For such a tran- 
sition, the A and L conversion coeffi- 
cients, «“ and x" respectively, may be 
written as, 


YUE 4 (1) Etap + (1 ae, 
\ 
\ 


2 3 4 5 6 


eee 


(2) a” = ap; L(1— ax. 


Eliminating between these equations, it 
: is seen that 
Fig. 6. - cé-(K/L) plot on point- 
nucleus (A) pug finite size nucleus n (Brat — of 62) K/L 
(S) basis. (>) Mi % 3 =_= 
(Pi 2 ) (Pi x) K/L 


Fig. 6 indicates «* as a function of K/L. The shaded region being the 
experimental one. Using the conversion coefficients Bi(i—X, L) and x/(j—K, L) 
of ROSE (?), the curve R is drawn. The other curve S is drawn, using Sf and a 
of SLIV and BAND (*), and the finite nuclear size corrected values (32) of 7 
and «;. The experimental values are seen, from this figure, to be in better 
agreement with Sliv’s values than with those of ROSE. 

Using Sliv’s values of 6* and «5 and our value of «*, the M1 component 
is found to be a=0.82-+0.10, so that the transition is predominantly M1. 
This is in agreement with the value of M1/H2—5 given by SUNYAR (21), in 
his summary of lifetimes of M14E2 type of transitions. 


(81) A. W. SunYAR: Phys. Rev., 98, 653 (1955). 
(32) Privately circulated tables of M. E. Rose. 
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Our thanks are due to Mr. BHAGWAT and Mr. PANAT for their help during 
the experiments. 


RIASSUNTO (*) 


I coefficienti della transizione a 145 keV nel Pr e di quella a 129 keV nel !9Ir 
sono stati determinati, risultando, rispettivamente 0.341.035 e 2.1+.3. Si trova che 


questi valori si accordano meglio coi valori per il nucleo di dimensioni finite che con 
quelli per il nucleo puntiforme. 


(*) Traduzione a cura della Redazione» 


2431 


JL NUOVO CIMENTO Vor. IX, N. 4 16 Agosto 1958 
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(ricevuto il 2% Aprile 1958) 


Summary. — The dispersion relations for form factors are studied in 
perturbation theory. In particular, the relations are proved for a real- 
istic nucleon in the form usually assumed in literature. It is demonstrated, 
however, that in more general cases the mass spectrum in the dispersion 
relations depends on the masses of the contributing particles in a peculiar 
way, which is different from the simple relation hitherto believed. Such 
a situation arises, for example, in the case of some hyperon models and 
bound systems such as the deuteron. A formula concerning the final 
state interaction to be used in conjunction with the dispersion relations 
is also derived from the invariance of the theory under space-time reversal. 


1. — Introduction. 


In this paper we are concerned with the properties of the form factor F 
of a particle 


(LE) KP | F(p — q)|p> = <q|ji(æ)|»p> exp [— i(p—q) -x]. 


Here j(x) is a local Heisenberg operator such as the charge-current density; 
|p> and |q> are true eigenstates of a particle characterized by four-momenta p 
and 4g, and other quantum numbers which are suppressed; |p® and |g°> are 
the corresponding eigenfunctions for an interaction-free particle. 

It has been generally believed that F(p — q) = F(k) should be representable 
in a form analogous to the dispersion relations for scattering matrices. For 


(*) This work was supported by the United States Atomic Energy Commission. 
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example, if the particle has no spin, (x) is a scalar density, and F(k) — 0 
sufficiently rapidly as k? + co, we should have (1) 


foo} 


(1.2) F(k) di -- d(me). 


m 
0 


mo(— k?) is the so-called absorptive part of F(k), which should exist only above 
a threshold — k? > mì for material creation out of vacuum by the external 
probing field of frequency k, and wave number k. 

In fact, such «dispersion relations » for form factors have been proved 
under some restrictive conditions from the general properties of field theory, 
including notably the local commutativity of field operators (2). For the form 
factors of a nucleon (mass M) interacting with pion field (mass 4), for example, 
equations of the type Eq. (1.2) are valid with 


(1.3) My = 20 
provided that 
(1.3') > (V2 —1)M. 


On the other hand, the relation was proved by the author in every order 
of perturbation in an example of quantum electrodynamics where m, = 0 (8). 
The present paper will deal with an extension of the proof in perturbation 
theory so as to include the case of nucleon form factors. In particular, it will 
be shown that the relation (1.3) is valid for actual nucleons without the res- 
triction (1.3’), which is not satisfied by the real pion. 

For practical applications, it is also important to relate the absorptive part 
of F to other observable quantities, thereby making the dispersion relations, 
at least in principle, equations to be solved for #. The main idea behind this 
is very similar to the so-called final state interaction in nuclear reactions for- 
mulated by WATSON (*). 

These problems will be treated in the following sections, specifically with. 
the nucleon form factors in mind. In the course of the development, it will 
also become clear that some modifications are necessary with regard to the 


(1) G. F. CHEw, R. KarPLUS, 8. Gastorowicz and F. ZACHARIASEN: Phys. Rev. 
110 265 (1958); M. GOLDBERGER and S. TREIMAN: preprints on form factors for weak 
interactions. Also a lecture by GOLDBERGER on nucleon form factors at Stanford 
Conference, 1957. 

(2) H. J. BREMERMANN, R. OEHME and J. G. TAYLOR: Phys. Rev. 109 2178 (1958). 

(3) Y. NAMBU: Nuovo Cimento, 6, 1064 (1957). This paper will be quoted as I. 

(4) K. WATSON: Phys. Rev., 95, 228 (1954). 
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usually held statement about the lower limit m, of the dispersion integral when 
more general classes of problems are considered. The discussion of this si- 
tuation will be made in the last two sections. 


2. — Representoation of the nucleon form factors. 


We consider here the realistic case of the baryon (nucleon and hyperon) 
field interacting with 7 and K meson fields (spin zero). The interaction is 
to be renormalizable. Among the particles, the nucleon and the pion are 
characterized by the fact that they are respectively the lowest states of the 
baryon and the meson families, which will turn out to be an important point 
in the following. Other details are not essential. 

We are interested in the electromagnetic form factor of a nucleon 


(2.1) <q°|F,(k)|p®> = <alj,(@)|p) exp [— i(p — 9)-2], 


j,(v) is the total charge-current density of the participating fields. From 
Lorentz and gauge invariance it follows immediately that 


€ 


(2.2) per = oF (Vy, + oy 


BY (Ke? tO yy ky 5 
(n = proton, p = neutron). 
The above decomposition can be carried out in perturbation theory for each 


Feynman diagram. Let such a contribution from a single diagram @ be simply 
called F,. Then we want to show that 


n a) ke i em?) 


por k2 Les ! S \ FEST 
1 (4?) | 9 po I À 
2,3) . (2)? 
os" (Mm?) 
TSO) = se mi . 
puo fran 
(24)? 


The different forms of the representation for 7, and F, are due to their dif- 
ferent behavior for large %?. It is also a reflection of the fact that 7, requires 
an infinite charge renormalization whereas F,, the anomalous moment part, 
converges without a corresponding renormalization. This is well known, and 
will not be discussed here. Another point to make is that the same function 
F, can represent creation or annihilation of a nucleon pair by an externat 
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electromagnetic field. In such cases k? will be < —4M?, whereas for a form 
factor k? > 0. 


21. — According to the results of I, #, and F, should have a representation 


F,(k) age oe ACER 


ie)" 


fe (ay) 
2 -5 fed ve 


J = ap? + Bq? + yk? =— (a+ PM + yk. 


These equations establish the analyticity of F(k?) for Im (— k?)> 0. In ge- 
neral, there will be a singularity where 


J+ M=0 or —k—=(1—x—B)M?}y. 


In comparing Eq. (2.4) with Eq. (2.3), the different powers of the denominator 
will not matter since they do not affect the region of singularity. 


2°2. — In order to determine the singularities of Eq. (2.4), we must use 
Theorem (5) of I, which says that under certain conditions one can find a 
reduced (lower order) diagram G, and a corresponding F,(k?) with J,(p, q, k) = 
= oep?+fhoq?+yk? such that 


(2.5) Ja, y, 2) > F(a, Y, 2) , 


if æ, y, 2 are real and x+y+<=0. and if the «mass term» M? is always 
normalized to the same value. 

Let us consider J(p, q, k) as a function of two independent vectors 
k=p—q and l=p+q. Write 


J(P, q, k) = J(k, 1) = Jak? +2J,k41 + dal? . 


In our case k:l= p?—q?=0, and Eq. (2.5) means 


(2.6) DM TAN 200% 


op} © 


Now both &k and / are time-like if 0 << — kh? < 4M’, so that with |k| =~ /— k, 
|l|=/—EP we get 
J(k, l) = — 1 |& |? Jo], 


(2.7) 
> Ji |b)? —dy, |b? = J,(h, 1). 


N 


614 Y. NAMBU 


Suppose J, +M?>0 at the same time for 0<—k?< m,<4M?. Then 
a fortiori 


dic MB SW. 


Once this is satisfied, 7+M? and J,+ M? will have no zeros for — k? < m, 
since they are increasing functions of k?, and Eq. (2.3) will have been proven, 


2°3. — Let us classify the Feynman diagrams in the following way. Since 
the nucleon number is conserved, a diagram consists of 


a) a baryon line starting with the initial nucleon line and ending with 
the final nucleon line, carrying a mass m > M; 


b) any number of baryon closed loops; 


b') meson (x or K) lines connecting the baryon lines, carrying a mass 
m > p. 


The external field («y ray ») may act on any of the charge carrying lines. 
So classify F into F, and F, according as the y-ray acts on the baryon line a) 
or the others. Then we carry out the following operations: 


1) Reduce the masses appearing in the À, functions on the baryon line 4) 
from their original values to M. (We may leave untouched the factor y-0 —m 
of S,) . 


2) Reduce the masses appearing in the A, functions in b) and b’) to 
the pion mass u < M. We will now refer to all these lines simply as « meson » 
lines. 


Under these operations, the mass term called M? in I obviously decreases 
(or at least does not increase). 


3) If the y-acts on the baryon line (case F,), remove all the « meson » 
lines except one, so that we get a typical 
lowest order diagram Fig. la. 


4) If the y-ray acts on one of the 
meson lines (case /,), remove as many 
«meson» lines as possible. We will 
end up with another typical diagram 


i Fig. 1b (*). 
Fig. 1. — Comparison diagrams for The final diagrams thus obtained we 
the nucleon form factors. now call @,. According to Lemma (4) 


(*) In case there exists a vector meson, the situation should be different because 
of different selection rules. 
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of I, J will increase (or at least not decrease) to J° during the operation 3 or 4. 
Thus after all the operations 1 to 4, the inequality (28) of I; J/M? < J,/M, or 


for our choice of normalization, is certainly valid. 

For the diagrams la and 1b, we can easily verify the dispersion relations. 

We find (*) 
| Mo = 2M (> 2u) Tore Lig 
(2.8) 
| PAT OT. 
This is a more precise statement than is indicated in Eq. (2.3). 

In the problem of nucleon form factors, a most puzzling point is the ap- 
parent big difference of proton and neutron charge radii, which implies near 
equality of isotopic vector and scalar parts of the charge form factor. Using 
a Similar technique as above, we can prove that the scalar part has, as has 
been usually claimed, m,= 3u in contrast to m,— 2u for the vector part. 


3. A theorem about the final state interaction (°). 


The results of the previous sections show that the absorptive part of F 
exists only if & is time-like and — k? > 4u?. As was noticed before, the func- 
tion F for — k? > 0 does not correspond to a form factor, but rather to a 
matrix element for the transformation of a virtual y-ray into a nucleon pair 
if — k? >> 4M?. In this region the absorptive part can be interpreted as arising 
from the interaction of the produced pair. Below the threshold, /(k?) must 
be considered as an analytic continuation, which can still have an absorptive 
part because real mesons may be produced in the intermediate states. 

The exact mathematical statement about the final state interaction will 
be derived as follows: Let us assume that the dispersion formulas obtained 
in perturbation theory are valid for the entire form factors, and consider the 
matrix element 


—p, a" |j,(k)|O), —& >4M, 


(*) In deriving the result (2.8), it may seem that a difficulty arises from vanishing 
of the denominators corresponding to the initial and final nueleon lines in Fig. 1 when 
we set p?=q?=— M?. This can be avoided by first setting p?=q? > — M? and 
taking the limit. 

(5) This is an adaptation of a theorem formulated for the case of scattering dis- 
persion relations by S. FuBINI, Y. NAMmBU and V. WATAGHIN; to be published. 


È 
Del 
+ 
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where <— p,q | is an outgoing state of a nucleon pair. Because of the in- 
variance under time reversal and space reflection, we have 


(3.1) 04 pi Uj,0-|0>* = OU, U-1|—#, 7") = 
= <— D, qu 112 |0> — DI <— D; q | n°" <n" li, |0> > 


{— P, "| here is the space-time reversed state (p,q >P,4, c-+-— 0), and 
U the associated transformation operator. » runs over a complete set of 
states, which includes the nucleon pair state. <— p, g'"|m°") is an element 
of the S-matrix S+= S-1. By the nature of the operator j,, n will consist 
of states of nucleon number zero, angular momentum one, and odd spatial 
and charge parities. For a nucleon pair, in particular, it involves 38, and *D, 
states, reflecting the fact that there are two invariants y, and 1,,k,. 
Let us now introduce the generalized form factor , by 


(3.2) Cn |F ,|0> = <n 19,19 , 


for an arbitrary |n>. As in the case of F,, F, may be decomposed into a 
sum of different invariants 


(3.3) Fi = > Osa ier 


where O, is a hermitian operator. «, are scalar numbers, other than k?, which 

are necessary to fix the state. For a given k? there will be only a finite number 

of invariants since only a finite number of particles can be present in |»). 
The relation (3.1) then becomes in this case. 


CUI IDEA SE nt rt IF ,,|9> 
m 


(34) 7 OF 
> <n ODF (He, an) [0) = > > <n? S+|m><m? OOF ,(k?, «,|0> . 


à mar 


We may write it symbolically 

(3.5) 2 OF: DO 

Alternatively, with S=1+iR, 

(3.5') > OL FEF) I 

On the left-hand side, i(¥* — F,) is twice the absorptive part of F. Assuming 


S to be given, this expresses a set of linear homogeneous equations to be 
satisfied by the absorptive and dispersive parts of ¥;. It must be understood 
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that they are meaningful even below reaction thresholds through analytic 
continuation. 

Eq. (3.5) or (3.5') and the dispersion relation (2.3) together constitute 
coupled integral equations for the absorptive and dispersive parts of F, and F, 
as functions of k?. Since other generalized form factors also appear in Hqs. (3.5) 
and (3.5’), it may be suggested that all ¥; should have their own dispersion 
relations like 

co 


12 (m2 m2 
(3.6) ONE LAURE | Im Fi(m?, x.) dm 


x} k? + m?—ie m? ’ 


Mi 


with an appropriate choice of the parameters x,. Then we would have a set 
of integral equations for all ¥. In fact, Eq. (3.6), may be easily derived from 
the representations of general Green’s functions. This will be shown in the 
Appendix I. 

We will not attempt to solve the integral equations for #. Except for 
the single channel case where only one 7, has an absorptive part, the only 
practical way seems to be perturbation expansion starting from Re ¥ = % (0). 
But the solution may not be unique. If there is more than one solution, the 
difference of two such solutions must satisfy homogeneous dispersion relations. 


4. — Form factors in a more general case (°). 


As is clear from the proof in Sect. 2, the derivation of the mass limit 
Eq. (2.8) depends critically on the fact that nucleon and pion are respectively 
the lightest particles of baryon and meson families. In case, for example, the 
form factor of a hyperon is considered, the proof will not go through since 
a hyperon can change to a lighter particle, and destroys the necessary inequa- 
lities in the proof. 

Such a situation can be handled as follows: Since a hyperon form factor 
is not of immediate practical interest, we will adopt a simplified model, con- 
sisting of a hyperon, a nucleon, and a K-meson field, with the respective masses 
M, > M> x. Then a form factor diagram will contain: 


a) A baryon line (or baryonic charge carrying line) connecting the initial 
and final hyperon lines. 


b) A strangeness carrying line, which may overlap with a), also con- 
necting the initial and final lines. 


(8) A very similar observation to the one in this section was recently done by 
KarPLUS, SOMMERFELD and WICHMANN: preprint. 
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c) Baryonic charge loops and strangeness loops. 


These lines will overlap each other to form a connected diagram. The 
external y-ray may act on any of the lines. 

In order to facilitate the comparison with standard diagrams, we will first 
observe the following theorem. 

Theorem. In order that the inequalities (27) of I 


T (By à) M2 (Gy) Ty (2, 0) Mix)" 


J (Gi) Br)/ MB) < Jodi, Be) MB) » 


(4.1) 


be realized for two diagrams @ and @,, and for any particular sets 2, x or q, p; 
it is necessary and sufficient that the following inequality be satisfied: 


(4.2) Fes) >). 
The function 7(z,;) is defined for a diagram by 


(42) 7 (€) = Min (4 Y mr, 


r; => 0 is the length of a propagation line, m, the mass of the propagating 
particle, and the summation is over all the lines that occur in the diagram. 
The minimum is taken with respect to the co-ordinates of the internal vertex 
points. The proof of the theorem rests on the property 


I, &) M) > F(z) , 


which will be derived in the Appendix. 

Now proceeding more or less in the same way as before, and making use 
of the above theorem, we find the following comparison diagrams (with the 
corresponding J,) 


1) If the y-ray acts on the baryon line a), take Fig. 2a. 


2) If the y-ray acts on the strangeness line b), take Fig. 2b. 

3) If the y-ray acts on a (baryon or strangeness) loop, take either 
Fig. 2e or 2c’, according as it is connected to the baryon line a) or the strange- 
ness line 5). The masses on the connecting lines have all been reduced to xa 


It is not difficult to see that Fig. 2e and 2c’ may further be replaced by 
+ 6 ‘ 2 6 Mica nti = qu ; . “uc È 
Fig. 2d and 2d’ respectively, which are combinations of the vertex diagram 
2a or 2b, and a self-energy diagram. 


The main interest is in the vertex diagrams. They are essentially the same 


2440 


DISPERSION RELATIONS FOR FORM FACTORS 619 


as Fig. 1, except that we now put — p?>=— ga = Mj > M?. By direct calcu- 
lation, the lower limit m? of the dispersion relation (2.3) now turns out to be 
different for the following three cases. 


M 
Le | sf 
| M x | t 
& 
x! k M De x k M | k 
| 74 | | 
| 4 | x 
M x | 
A M 4 
p p p p 
(a) Rio) (Cc) Cer) 


(d) (d') 


Fig. 2. — Comparison diagrams for the hyperon form factors in our model. Solid line: 
baryonic charge line. Broken line: strangeness line. These are not necessarily actual 
diagrams. 


Cosel... Ti Mi =< M Ex, then 


(4.3) Mo = 2M Fig. 2a, 
= 2% Fig. 2b and 2d. 


The overall lower limit is thus m)= 2x. 


Case 2. If M?+x2< Mi <(M+x), then 


A? 1 
(4.3') Me = (200)? vin (2205) HigwZa,, 
A? 
= (2%)? — WE (> 0) Fig. 2b and 2c, 


A2 
= Min (2%, (2m)? — =) Fig. 2d’. 
x 


A = M?— M?—x’. 


The overall lower limit is m{= (2x)?— A?/M?. 


441 
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Case 3. If M,> M- x (the hyperon unstable), then 


(4:31) mi = — co (*). 


This result indicates that the usual argument for the lower limit m, as given 
in Sect. 2 does not necessarily hold for each Feynman diagram, nor for each 
order of perturbation expansion. 

In the general derivation of the dispersion relations from properties of local 
field operators, the spectrum of the total energy-momentum of intermediate 
states is specified, but not the masses of individual particles. Thus in the 
present example we would regard M,=M+% and 2x (or 2M) as the lowest 
total masses of the intermediate states. In order to have the Case 1 above, 
it is then necessary that 


M<M +», 


for fixed M, and any x< M,. Since the right-hand side takes a minimum 
value M°/2 for x= M — M,/2, that means 


(4.4) M°>2M5 cor Wie (V2 Dae 


This coincides with the condition (1.3’) obtained by BREMERMANN, OEHME 
and TAYLOR (?) for the derivability of the dispersion relations on the basis 
of local field theory. It would seem, therefore, that their condition (1.3) is, 
at least partly, caused by their inability to specify the detailed nature of the 
intermediate states, rather than their inability to take account of the uni- 
tarity condition in general. 


5. — Concluding remarks. a} 


We have examined in some generality the validity of the « dispersion re- 
lations » for form factors in perturbation theory. In particular, it has been 
shown that the nucleon form factors should satisfy the usually assumed re- 
lations. Thus the large difference of the observed proton and neutron charge 
radii still remains to be understood. On the other hand, it has become clear 
that the mass spectrum in the dispersion integral does depend on the masses 
of the individual particles in a peculiar way. 

The last point is of practical significance since any qualitative as well as 
quantitative appreciation of dispersion relations depends on the mass spectrum. 
It would be necessary to understand the mass spectrum conditions from a 


(*) The actual values of M,, M, x correspond to the case 2. 
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more general point of view. The relation between causality and the validity 
of the dispersion relations with the modified spectral conditions is not yet clear. 

Another interesting point turns up when we consider bound states. It has 
been shown recently (7) that one can describe bound systems in field theory 
in much the same way as the so-called elementary particles. One may expect 
then that the form factors of a deuteron, for example, should satisfy the same 
kind of dispersion relations as the nucleon. On the other hand, we know that 
the deuteron form factor should be determined largely by the binding energy 
according to the uncertainty principle. 

It would be absurd if in this case the mass spectrum of the dispersion 
relation started from m,= 24, as would be expected if the original argument 
about m, was correct. As a matter of fact, if we calculate the deuteron form 
factor from the virtual process (M,= M,= M assumed for simplicity) 


with local interaction, we do fall in the Case 2 of Sect. 4, with 
mì = 4M° (1 — Mi/4M?), 


~16eM , 


where e is the binding energy. Thus the radius of the deuteron should be 


falls off as exp [— VeMr]/r. (r is the proton-neutron distance, whereas the 
radius is measured from the center of gravity). 


APPENDIX I 


Representation for generalized form factors. 


A generalized form factor F, as defined in Eq. (3.2), is essentially a certain 
n-point Green’s function G,(k, Po; …, Pn) (n > 2), considered as a function of 
one momentum k, the other n —1 momenta p, being those of free particles. 
In other words 


(A.1) FU, pi) = Ly (h) IT LP) Gah, pi) , 
where the L,’s are the free-field differential operators, and pi = — mi. 


According to the results of I, each invariant ¥;, Eq. (3.3), will therefore 


(7) K. Nismigima: preprint; W. ZIMMERMANN: preprint. 
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have a typical representation lke 


a Trargi È fe ì I m(Br) A 
Fo (k, Pp) sò O, SARI DARI af Whee , 
(A2) Je pr +) pp 1 S Bip: + ph + DBs. (D; ED; + Px)? +. 
i=2 i>j i>i>k 


The summation continues until all combinations of n — 1 momenta p,: Note 
that we have deliberately avoided the explicit appearance of combinations 
like k+p,;, k+p,+p,;,.... by replacing them by the complementary sums. 
Thus 


2.n 


(A.3) T=Bh4 YD yd Ps = Bik — yum +2 DvD Ps 
id i i>j 
where f, > 0, y,; > 0 and the matrix y;; forms a non-negative quadratic form. 
Eq. (A.3) contains n+(n—1)(n—2)/2=1+n(n—1)/2 scalar products, 
whereas only n(n — 1)/2 of them should be lineraly independent since 


(A.4) k? = (— > Di) =) Pi Di. 


2,) 


So eliminate, for example, p,-p,; by the above relation, or 


(A.4') 2P2° Ps = k? — DI Pi Pi - 
Then J takes the form | 
J =ck?+ J'(a,), 


(A.5) (6) = Bi ty» 20, 
a, = Pi Pi, t 22, je, 1-07. 


Thus for fixed parameters «,, F;(k?, x.) is analytic for Im (— k?) > 0 as was 
expected. The problem remains as to the lower limit m, of the mass spectrum 
in this case. We will not try to determine it since at any rate such gene- 
ralized form factors will be of little pratical significance. 


APPENDIX II 


Proof of the Theorem in Sect. 4. 
By the results of I, the quantities [(z) and M? are 
ek x 
I(z;) = Min 5 Ja, 


(A.6) Me 
My 

Mi = y 

33° 
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a.>0 is the parameter attached to a propagation line e, in the diagram. 
The minimum is taken with respect to the variation of the internal vertex 
co-ordinates. 

Now we want to find the double minimum J(z;) of I(z;,a)M?(a,) with 
respect to both the co-ordinates and the parameters a,. By taking variation 
with respect to a,, we immediately find 


Ns / = 
(A.7) oi = VG 4 (eva) 
and 
(A.T') IM°?=(3Y mr). 
Thus ; 
(A.8) Ja) Min (ar 


the minimum now being with respect to the internal co-ordinates in the 
diagram. 
For any two diagrams, for which one can realize 


(A.9) T(z, a) M°(a) > Lo(2, do) Mo(4) ; 
for ony given set of parameters a, it follows that 
| IM? > 1, M? > Min, I, Mo = To), 


boy) | Min, IM? = F(z) > Lyle). 

Conversely, if F(z) < 7,(z), this does not necessarily mean the inequality (A.9), 
to be valid for fixed a, a, and all 2. However, for any fixed a and 2, there 
must be some set a,(a, 2) such that (A.9) is satisfied since otherwise 7(2) > F(z) 
would not hold. This means in turn that for any fixed a, and q; there is some 
sei &,.(@, gq) such that 


J(a, q)|M2(a) > Todo; /Mi(%) ; 


is satisfied. For our present purposes in Sect. 4, this should be sufficient. 


RIASSUNTO (*) 


Colla teoria delle perturbazioni si studiano le relazioni di dispersione pei fattori 
di forma. Si dimostrano in particolare tali relazioni per un nucleone reale nella forma 
usualmente assunta in letteratura. Si dimostra, tuttavia, che nei casi più generali lo 
spettro di massa comparente nelle relazioni di dispersione dipende dalle masse delle 
particelle che danno contributo in un modo particolare, differente dalla semplice rela- 
zione finora ritenuta esatta. Tale situazione si presenta, per esempio, nel caso di qualche 
modello d’iperone e di sistemi legati quali il deuterone. Dall’invarianza della teoria 
rispetto all’inversione dello spazio-tempo si deriva anche una formula riguardante 
l'interazione degli stati finali da usare in unione alle relazioni di dispersione. 


(*) Traduzione a cura della Redazione. 
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A°- and 0°- Particles Produced in Iron. 


G. ALEXANDER (*), C. BALLARIO (+), R. BIZZARRI (*), B. BRUNELLI, 
E. DI CAPUA, A. MICHELINI, G. C. MonETI and A. ZICHICHI (*) 


Istituto di Fisica dell’ Università - Roma 
Istituto Nazionale di Fisica Nucleare - Sezione di Roma 


(ricevuto il 23 Aprile 1958) 


Summary. — The results of an experiment on some properties of A. 
and 0°-particles produced in the plates of a cloud chamber are discussed. 
115 V°events are analysed by means of angle measurements and ion- 
ization estimates. The frequency of production of A°, 6°-pairs is found 
to be small. For neither type of particle is any evidence found for ani- 
sotropy in the decays. The V° production cross-section appears larger 
for a pion primary than for a nucleon, and increases with the energy 
released in the shower. 


1. — Introduction. 


For some years the masses and the decay schemes of A° and 0° particles 
have been well-known mostly from measurements in magnetic cloud cham- 
bers (1). Although some decisive results have been recently obtained for A° 
particles (*), much room for improvement remains regarding the knowledge 
of their other properties, such as spin, mean lifetimes, behaviour with respect 
to inversion operators, and the details of their processes of production. 


*) Now at Dublin Institute for Advanced Studies, Dublin. 
+) Now at CERN, Geneva. 


fe 
( 
(*) Now at «Enrico Fermi» Institute for Nuclear Studies, University of Chicago 
(!) See, e.g. the paper by C. FRANZINETTI and G. MorpurGo: Suppl. Nuovo Cimento, 


6, 469 (1957), chap. 3 and 4. 
(*) Ibid., chap. 15. F. $. CRawrorp et al.: Phys. Rev., 108, 1102 (1957); F. EisLER 


> 


et al.: Phys. Rev., 108, 1353 (1957); F. EIsLER et al.: Nuovo Cimento, 7, 222 (1958). 
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This paper gives the results of a study of some of these latter properties. 
The investigation was carried out on 107 A° and 0° particles produced in the 
iron plates of a multiplate cloud chamber exposed to cosmic radiation at an 
altitude of 3500 m above sea level. 

The dominant decay modes of A° and 0° particles are well known to be: 


(1) Ko peer M eV 
and 
(2) 0° > ++ n + 214 MeV. 


If other neutral particles exist which disintegrate into two charged bodies, 
these latter particles are rare (1). It is therefore possible to make use of decay 
events that are not completely measurable. A V° decay into two secondaries 
that is definitely not a A° decay (1) can reasonably be assumed to be a 6° decay (2) 
and viceversa. Events selected in this way are of value for statistical studies 
giving information about relative frequencies of production, angular and energy 
distribution and mean lifetimes. 

In our experiment we considered only V° decays for which an evident origin 
existed in the plane of the V° secondaries and we measured the angles 0, 
and @, which the decay products of the V° particles formed with its line of flight. 
From these angles it is possible to calculate the ionization of the two prongs 
for the possible decay schemes A° and 6°. If, as often happens, these ioni- 
zations are sufficiently different in the two cases, the nature of the 9-particle 
can be established by simple visual estimation of the ionizations. The method 
of analysis and the selection criteria, as well as the bias introduced by them 
into our sample of V°-events are discussed in Sect. 3. 

Sect. 4 presents the results obtained concerning the relative production 
frequency of A° and 0°-particles, Sect. 6 those relating to angular distributions, 
and Sect. 7 those relating to the parent showers. 

Our estimate of the lifetimes of A° and 6°-particles has already been pub- 
lished (*) and is referred to in Sect. 5. Sect. 2, which follows, describes some 
features of the cloud chamber and of the reprojection system; they have been 
described in detail elsewhere (4). 


(3) C. BALLARIO, R. BIZZARRI, B. BRUNELLI, A. De Marco, E. Di Capua, A. MI- 
CHELINI, G. C. MonETI, E. ZAVATTINI and A. ZIOHICHI: Nuovo Cimento, 6, 994 (1957). 

(4) C. BaLLarIo, B. BRUNELLI, A. DE MARCO and E. ZAVATTINI: to be published 
in Nuovo Cimento. 
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2. — Experimental apparatus and measurements. 


The useful volume of the cloud chamber is (3035 %20)cm*. As Fig. 1 
shows, there are four iron plates, each of thickness 12mm. The figure also 
shows the arrangement of the triggering counters and of the lead transition 
layer. 


CTW 


Counters 1:27 cm 


VIZIZIZIIZZIAZZIZZIZZZA 
(Z2Z7L LLL LL LLL LIL LL 22 272 


(P2722 LZ LLL LLL ZL LLL 22224 


(77777777 LL LIL LL TLL 777} 


Counters 2*35 cm? 
lo)» CSOSSCOOSOHOES. | Sie eee 


iron lead perspex glass 
0 5 10cm 20 


Fig. 1. — Cloud chamber and trigger system. 


The chamber was triggered by either of the following two coincidences: 
a) between at least 2 counters of tray A and at least 3 counters of tray B; 
b) between at least 4 counters of tray A and at least 2 counters of tray B. 


With a chamber dead time of about 3 minutes, the rate of expansions was 
15 per hour. 

The chamber was photographed by two cameras. Stereoscopic photographs 
were reprojected by a system similar to that used for photography. 

There are various possible methods for the spatial reconstruction of the 
events (1). In most cases the angles 0, and 0, were traced directly onto a screen 
on which the two stereoscopic pictures had been superimposed. A correction 
was always made for gas motion during expansion. 

For the purpose of estimating the errors in reprojection and calculation, we 
chose four typical events, and for each of them made ten calculations of the 
value of 0 (=9,--0,). We found that the half-breadth of the distribution of 


2448 


G. ALEXANDER, 0. BALLARIO, R. BIzzARRI, B. BRUNELLI, E. DI CAPUA, 
A. MICHELINI, G. C. MonETI, and A. ZICHICHI 


Typical photograph with two indipendent V°-events. 
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values was always less than 1 degree. So we generally assumed an error of Ms 
for some events in awkward positions we increased the error to 1.5°. The errors 
in the partial angles were mainly due to the indeterminacy of the origin of 
the star. 


3. — Selection and analysis. 


31. General procedure. — With the equipment described we took about 
30000 photographs, which were scanned for V°-events. 

It was decided to regard as neutral hyperon or heavy meson decays all 
those V°-events which were associated with a nuclear interaction in iron, co- 
planar within 3°, and satisfying the following conditions: 


a) 1) At least one of the two tracks showed an ionization exceeding twice 
the minimum; 07 


2) At least one of the two tracks penetrated a plate without multi- 
plying and without being deviated by more than 20°; or 


3) The angle 0 between the two tracks was bigger than 20°; and 


b) The two secondary particles displayed no perceptible multiple scattering 
in the gas; and 


c) Both branches had at least 12 mm of measurable track; 


d) The vertex was in a position where the scanning efficiency could be 
assumed to be constant and there was only a negligible chance that a fortuitous 
coincidence of two tracks might simulate a V°-event (5). 


Conditions a) and b) were designed to eliminate electron pairs; condition d) 
made it possible to calculate the probability of every V°-event being observed (5) 
With our chamber this probability varies considerably from event to event 
and this had to be allowed for. 

Concerning condition a 2) it should be noted that the thickness of each 
plate is 0.7 radiation lengths. It may be objected that this is not thick enough; 
however, in the case of 8 out of 11 V°-events admitted on this criterion (see 
Table I), one or both secondary particles traversed at least two plates. The 
critical angle for the deviation of the secondary particles in their passage through 


(5) Cf. C. BaLLARIO, R. BIZZARRI, B. BruNELLI, A. DE Marco, E. Di CAPUA 
A. MICHELINI, G. C. Moneti, E. ZAVATTINI and A. ZICHICHI: Nuovo Cimento, 6, 994 
(1957). 

(6) D. B. GayTHER and C. BUTLER: Phil. Mag.. 46, 467 (1955). 
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the plate was taken to equal twice the mean square deviation angle of a pion 
leaving the plate with double the minimum ionization. 


Taste I. — Distribution of the events admitied by virtue of the various alternatives 
of condition a). 


Alternative i Number of events 
Il 1 
2 11 
3 15 
1 and 2 2 
Inde | 14 
DENON © 50 
Pan | 22 


A total of 115 V°-events were selected according to these criteria. All the 
photographs were examined twice, and each scanning for events satisfying 
the selection conditions had a 95% efficiency. In these circumstances the pos- 
sibility that any event may have escaped detection in both scannings may 
be neglected. 

The photographs were stereoscopically projected, as described, and each 
event was geometrically analysed. We measured: 


1) The angle 0 between the two secondary tracks; 


2) The angles 0, and 0, between the line of flight of the V°-particle and 
each of the two secondary tracks. In the following we shall assume that. 
8, < A; 


‘ 


3) The path length of the secondary particles within the iron plates; 


4) When the primary particle of the star which generated the V°-particle 
was charged and visible, we measured the angle ® between the line of flight 
of the V°-particle and the track of the primary particle, the angle between 
the production plane, è.e. the plane defined by the primary track and the line 
of flight, and the decay plane. 


Finally we estimated visually the ionization of the two secondary tracks. 
The V°-events were classified according to the following method. From 
the angles 0 and 0, the ionization of the two branches of the V°-event can 
be calculated for each of the two types of decay. The visual estimate of ioni- 
zation enables us to choose between the two interpretations A° and 6°, when 


JE, NO Too Ft: 
3 ely Saad 5) re 
x Tro» Tone — (0 <3) 
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where Lx, Los Lines Inq, are the ionizations of the two tracks according 
to the two interpretations. 

Equation (3) is the quantitative expression of the assumption that the 
ratio of the ionization of one track to the known ionization of another track 
can be visually estimated only within a factor of 2. The equation is in 
a form which will facilitate the discussion of the bias introduced by our method 
of classification, and will also enable us to calculate, by statistical methods, 
how many of the unidentified particles were A° or 0°-particles. 

When equation (3) is not satisfied both interpretations are possible, unless 
information about the range or the ionization change in the passage through 
the plates enables the two types of decay to be distinguished. 

For a first rapid analysis the curves in Fig. 2, 3 and 4 are very useful. 
These curves represent 0 as a function of 0, for A°-particles, and OL G4, Lor 


TF in 
DI LISI een ORES ASA 
À 1620220020 36% 30 32 30 36038 ZO 42 44 460 G6) 50, 52 54 56 
Pigs 2 = O decay, 0 vs. 0, (i.e. the angle between the proton track and the line of 


flight of the V°-particle) for various values of the proton’s momentum p.. Angles are 
in degrees, p in 105 eV; ionization is expressed as a ratio to minimum ionization. 


00-particles, the momentum (and hence the ionization) of one secondary par- 


ticle being assumed constant. 
In a more detailed graphical analysis the experimental point defined by 
the angles 0 and 0,, for each V'-event was marked in the plane (x, €) (7). 


(7) J. PODOLANSKI and R. ARMENTEROS: Phil. Mag. 45, 13 (1954). 
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Fig. 4.— 0° decay. 0 vs. 0, for various values of the momenta P and p, of the secondaries. 
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This makes it possible to calculate rapidly, for each interpretation, the mean 
values and the corresponding maximum errors of the momentum p ot the 
V°-particle, of momenta p, and p, of the secondary particles, and of the cosine 
of the emission angle 0* of the positive secondary particle in the centre of mass 
system (*). 

The following values were assigned to the constants of the two decay 
schemes: 


(A) a* e* 
“AG 37 MeV 0.693 0.180 
0° 214 MeV 0.000 0.33 


By the method described we obtained the following classification of the 
115 V°-events: 


- === = = 
Number of events for which 
ra | Total number ‘ A : ts to I 
Identification ; range information was needed 
of events SAC E RAR 
for identification 
PAU 27 & 
| 
0° 42 2 
A° or 0° 38 ae 
Neither A° nor 0° 8 = 


32. Apparently «anomalous » events. — In most cloud chamber experiments 
V°events have been observed that are incompatible with the decay schemes (1) 
and (2). It is now well established that 3-body decays of mesons occur and 
many «anomalous» V'-events are probably of this type. In our experiment 
the size of the chamber and the selection requirement that an origin for the 
V°-particle should lie in the plane defined by two charged secondary particles 
makes it improbable that our «neither A® nor 0°) events should be due to 
3-body decay processes. Unfortunately, we have insufficient measurements 
of the individual events to enable us to draw any conclusions regarding their 
nature. 

There remains the possibility of a statistical explanation of the presence 
of such events without recourse to the hypothesis of anomalous decay. There 
may be a number of reasons: 


1) The shower in which the V°-particle is generated may not be the one 
with which the analysis associates the event. The origin of the latter shower 


(*) For a 0°-decay it is actually not possible to distinguish the positive secon- 
dary track from the negative one; hence we can only get the absolute value of cos OF 
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may have been in the decay plane by accident, while the former was not visible 
in the chamber. 


2) The margin of error allowed for ionization estimates and angle measu- 
rements may be insufficient, because in exceptional cases the error exceeds 
the maximum values found from a number of standard events. 


A simple calculation based on the frequency of the showers found in the 
photographs which show a V°-event leads to the conclusion that about one 
spurious event is expected in our sample due to the first cause. 

The second cause seems even more plausible. All the events classified as 
«either A° nor 0°) have some geometrical peculiarities, generally tracks with 
a steep inclination towards the window of the chamber, which raise special 
difficulties both for ionization estimates and for stereoscopic projection. It 
seems reasonable, therefore, to assume that at least some of the « neither A° 
nor 0°» events are produced by abnormally large errors of measurement. 

Other causes, such as two-pronged stars produced by neutral primary 
particles are less plausible. 


3°3. Bias in selection and analysis. — At this point it will be well to examine 
and discuss the bias introduced by our criteria of selection and analysis. 

Given the decay scheme, each selection condition defines a region of the 
plane (cos 0*, p) which comprises all the points, and only the points, corres- 
ponding to the events meeting the condition. 

Moreover, if we neglect the few events that could not be identified without 
range information, equation (3) defines a region of the plane (cos 0*, p) cor- 
responding to the identifiable events. 

Fig. 5 and 6 show the curves delimiting these regions. We have drawn 
curves corresponding to two different values of ionization, range and total 
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angle, so as to see to what extent these criteria are critical. Note that for the 
A°-particle the curves I,=2 and I,=2.5 nearly coincide with the curves 


(3) Se = oral A 


(The curve (I, ,0/1,g0)(I,g0/Izx0) = 4 is of no practical interest). 


5 10 15 20 25 30 35 40 


Fig. 6. — (cos 0*, p) diagram for 0° decay (see text). 


An examination of Fig. 5, shows at once that for the A°-particles the se- 
lection criteria introduce no bias with respect either to the momenta or to 
cos 0*, while the impossibility of distinguishing 6°-particles cuts off the spec- 
trum between 0.7 and 1.0 GeV/e and introduces a slight bias against small 
values of 0*. 

Fig. 6 shows that for the 0°-particles the selection condition 0 > 20° intro- 
duces a fairly definite bias with respect to the momentum, by excluding mo- 
menta exceeding 2.4 GeV/c. However, this bias applies only to events the 
position of which in the chamber is such that it is not possible to establish 
whether the secondary particles do or do not penetrate through the plates. 
On the other hand, the identification equation (3) introduces a bias against 
small values of 6* and high momenta. 

As regards the scanning, we can say that, looking at the few V°-events 
first detected in the second scanning, the scanning does not seem to introduce 
any noticeable bias against any particular category of events. 

Finally, as regards the number of events analysed on the assumption of 
a line of flight mistakenly associated with the event by reason of chance co- 
planarity, we estimate that this number is about unity. 
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Before going on to present the results of our experiment, we must attempt 
at least a statistical determination of the nature of the V°-events classified 
as « A° or 0°». 

The curve Z,\,L,g/L,gL,ns — 2, Which we shall call the «identification cur- 
ve», never has zero ordinate (see Fig. 6). This fact enables us to estimate the 
percentages of A° and 6°-particles in the category in question. If we introduce 
the hypothesis of isotropy of the 0°-decay in the centre of mass system, the 
ordinate of the identification curve for any given momentum represents the 
probability that our method should lead to identification, without range infor- 
mation, of a 0°-particle with that momentum. 


ms 
cos @ 


BRE 
/ | (a a Momentum in 10° ev/c 
0 


5 10 15 20 30 40 50 


Fig. 6-bis. — Distribution of 6° events and of « A° or 0°» events, analysed as 60, in the 

plane (p, cos 0*). The represented errors are due to the errors of the angle measure- 

ments. All the events that, within the errors, are represented by points above the 
curve, are classified as « A° or 0°» events. 


Let n(p)Ap be the number of V°-events identified as 0° without range 
information, with a momentum in the interval Ap around p (*) and h(p) be the 
corresponding (average) ordinate of the identification curve. The corrected 
number of 0° in the interval Ap is n'(p)Ap= (n(p)/h*(p)) Ap, where h*(p) = 


= h(p) +4 A cos 0*, A cos 0* being the average value of the errors with respect 
to cos 9* for 0°-events represented in the plane (p, cos 0*) by points lying on 
the identification curve in the interval Ap (+). 


(*) The 6° events identified by means of range information are in this context 
included with the « A° or 6°» V°-events. 

(+) Since the errors with respect to cos 0* are due to almost constant errors of 
angle measurement, the former errors are fully determined by the co-ordinates of the 
decay in the plane (p, cos 6*). 
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The number of 6° among the events which were classified as «A° or 6°» 
would then be 


N(0)=Y 


h*(p,) —n(p.) ane 


and the number of A° would be 
N(A°) = N(A or 6°) — N(6°). 


Setting apart the two 0° identified by range information, we find that our 
38 «A° or 0% include only 4 to 8 0°, all with a momentum exceeding 7-10* eV/e 
and with cos 0* between 0.8 and 1.0. This means that if we regard all 
« A° or 60» V°-events as A° and add them to the 27 certain A°-events, we should 
expect the resulting 65 A°-events to have something like 8 per cent 6° conta- 
mination. This is insignificant for the shape of the histograms which we shall 
presently examine. 

Note that this result rests in the assumption of isotropy of the 0°-decays. 
We shall, in Sect. 6, discuss to what extent this hypothesis is consistent with 
the experimental results so far obtained; here, it may merely be stated that 
unless the deviations from isotropy of our 0°decays (see Fig. 11) are due 
merely to statistical fluctuations, they would tend to diminish the percentage 
of 0° in the V°-events classified as «A° or 6°». 


4. — Relative production frequencies. 


The frequencies discussed so far are the frequencies of observation of JIN 
and 6° events with our apparatus. As we have stated at the outset, we limited 
ourselves to V°-decays generated by coplanar showers and we subdivided those 
events into decays of the type (1) and (2). 

Now, while it is highly probable that these are the only two-body decays 
with charged secondaries, it is also known that A° and 0° particles can disin- 
tegrate in a different manner, that is into two neutral bodies, two charged 
bodies and one neutral one, and probably into three neutral ones. Moreover, 
the §°-particle almost certainly has no determinate mean lifetime but is a 
superimposition with equal amplitude of 0° and 03. Since only 0} particles 
could disintegrate into two charged or neutral pions and 0} would have a 
mean life at least 500 times as long as 0} (**) and almost all of them would 


(8) F. EisLer, R. PLANO, N. Samros, M. SCHWARTZ and J. STEINBERGER: Nuovo 
Cimento, 5, 1700 (1957). 

(9) M. Borpon, W. CHINOWSKY, M. Fucus, K. LANDÉ, L. M. LEDERMAN and J. 
TIinLor: Nevis 61. 
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thus escape from the chamber, we can assume that our 0° particles are all 0}. 
The probability that contamination may simulate a two-body decay com- 
patible with the decay schemes (1) and (2) is negligible in the light of our 
discussion in Sect. 3. 

To calculate the probability of a A° decay (1) and of a 0° decay (2), we 
can use the experimental results of EISLER et al. (8). We obtain 


(4) PAD; ro 0,032. 0.058 
and 
(5) Poo(x*, x) = 0.86 + 0.06. 


The probability of a 0° decay (1) is found by halving this last value, to 
take account of the probability of a 0° particle being 0°: 


(6) Part, 7) = 0.43 + 0.03. 


Now, it should be noted that owing to the size of our chamber and the 
spectrum of momenta of the V°-events observed, the probability of obser- 
vation varies considerably as between various events. To arrive at the cor- 
rected ratios of production frequencies, each event has to be weighted with 
a statistical factor w equal to the inverse of its observation probability (°). 

For the 27 «certain A°» we get: 


D'un 32", > wi = 524 
and for 42 «certain 0°»: 
37) 09 > w= 274. 


The weighted numbers of V°-events are 


These figures still have to be corrected to take account of the A° and 
0° decays included in the V°-events classified as «A° or 6°». To this end we 
can repeat the calculation of the preceding section, applying it this time to 
the weighted spectrum of momenta of our « certain 0°» decays. We find that 
the total weight corresponding to the number of 0° among the « A° or 6°» 
has the value 20 + 6. 

Now we can attribute to the approximately 32 A° included in the « A° or 6°» 
group the average weight of the 38 «A° or 0% events analysed as A°. This 
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weight can obviously not be appreciably affected by the small contamination. 
In this manner we get for the weighted number of A° among the « n) 
93 -|-22. The error allows also for the (not statistical) errors in momentum 
measurement, which are altogether negligible for the « certain N° » and « cer- 
tain 0°» categories. 

Let us now consider the distribution of the angles © (see Fig. 7). To allow 
for the events which remained unobserved because ® = 90°, we must increase 
the weighted number of A° by about 20% and that of 60 by about 10%. 
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Fig. 7. — Distribution of the angle ® for A° and 9° events. The shaded area represents 
the contribution of 19 unclassified events. The errors shown are calculated as es wî) È, 


è 


Another bias tending to reduce the frequency of V°-events is due to our 
selection criterion a) 3), which requires a total angle of more than 20° between 
the two tracks of a Vt-event, in the absence of secondary tracks seen to pe- 
netrate the plates or having I > 21 min. Consider the spectra of momenta 
of the 65 A° and the 42 0° events and compare them with the curves in Fig. 5 
and 6 which delimit, in the respective planes (cos 0*, p), the region in which 
the points satisfying the above-mentioned selection criterion are located. It 
will be seen that the bias has practically no effect at all for 0°, and a negligible 
one, in relation to the other bias effects previously considered, for FASI 

To determine the true production frequencies, those corrected figures finally 
have to be weighted with the inverse of the probabilities (4) and (6) of charged 


decays; we find 
Ni = 309 3279, 


Ny = 305 +67. 
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Hence the ratio of production frequencies is 


ay NO Orde 
ne 


This figure should be compared with various cosmic ray and machine results, 
which have been fully discussed by COOPER et GIO 


45 


numbers 
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Momentum in 10 ev/c 
Fig. 8. — Distribution of momenta for 65 A°-events. The shaded area represents the 
contribution of 36 « A° or 6°» V®-events. Errors as in Fig. 7. 


The cases of associated V° production which we observed had best be 


treated separately. There were three showers associated with a pair of V°- 


events that satisfied our selection conditions. In only one of the three cases 


we were able to identify both V°-events, one as A° and one as 0°. In 
each of the other two cases only one V-°event was identified, and it was 0°. 


(10) W. A. Cooper, H. FILTHUTH, L. Montanet, J. A. NEWTH, G. PETRUCCI, 
R. A. SALMERON and A. ZicHicHi: to be published in the Nuovo Cimento. 
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We can proceed as for the identification of the events classified as « A° or 0°», 
and conclude that the two incompletely identified pairs are almost certainly 
made up of one A° and one 0°. Assuming then that they are (A°, 0°) pairs, 
we get 

NA) 


40 


30 


weighted numbers 


20 


erre i e O e 16 18 20 22 26 26 28 30 32 34 36 38 40 


Fig. 9. — Distribution of momenta for 42 0°-events. Errors as in Fig. 7. 


This figure has to be weighted and corrected as above leading to the result 


N'(A°, 0°) = 79 +65. 


In speaking of A° production, we have so far included A° particles produced 
directly or by X° decay. The ratio between the frequency of A° particles pro- 
duced without being paired with a 0° (A° sing) and of those produced together 
with a 0° particle (A° ass) is 

N(A° sing) 


4.5 
SII OO 


N(A° ass) 


Tf we assume that all the observed A° particles are produced by x or n and 
that each final state has the same frequency as its charge symmetric state, 
the value of the ratio should be unity. 
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The deviation may of course be due simply to statistical fluctuation. How- 
ever, it is interesting to consider its significance if it were real. One possible 
cause of charge asymmetry in the final states may be asymmetry in the in- 
cident flux of the 7-mesons and nucleons. It will readily be appreciated, 
however, that a marked effect of this kind is hardly probable; another pos- 
sible cause is the presence of =~ and K-particles in the incident flux. 

It should be noted, finally, that in the calculation of the ratio under dis- 
cussion the effect of 0° absorption in the iron plates has been neglected. If 
we assume that the absorption length coincides with the geometrical one, the 
value of the ratio would be diminished by only about 10%, that is, by much 
less than would be required to reduce it to unity. 

If the observed charge asymmetry were to be attributed entirely to ab- 
sorption, then we must assume that the absorption length is about + the geo- 
metrical one, but the available data on the interaction of K and K particles (11) 
give cross-sections never exceeding the geometrical. 


5. — Mean lifetimes (72). 


We estimated the mean lifetimes of 25 of the 27 «certain A°» and of 38 
of the 42 «certain 0°) particles. The particles excluded from the estimation 
sample were analysed on the assumption that one of the secondary particles 
stops in the iron plate. This kind of information is not always reliable, be- 
cause of possible interactions in iron with emission of neutral products only. 

A discussion of these estimates, and particularly of the errors involved, 
has been published (*). Here we merely repeat the results, namely 


T A0 = (2.143 0:40) -10-10 Ss ; 


Tes == (0.92 p28)" 1051 Bs 


6. — Angular distributions. 


The angular distributions I(p) of the angle gm between production plane 
and decay plane, and /(9*) of the angle 0* formed by the tracks of the decay 
products and the line of flight of the V°-particle in the latter’s center of mass 


(41) W. K. H. Panorsxy, V. L. Fircx, R. W. Morzex and W. G. CHESNUT: 
BNL 3410. 

(12) For a summary of earlier results and literature see: C. FRANZINETTI and G. 
MorPurGo: Suppl. Nuovo Cimento, 6, 469 (1957); more recent estimates are given by 
CooPER et al. and by C. Graves and D. GLASER, circulating preprints. 
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system, have been studied by various authors in connection with the pro- 
perties of V°-particles (1°). 

The angle g is described by the vector k(z-) x k(V°), perpendicular to the 
decay plane, when, by counter-clockwise rotation from the direction k(V°), it 
is moved to the position occupied by the vector k(prim) x k(V°) where, k(prim), 
k(V°), k(r-) are the momenta of the primary particle, the V°-particle and the 
m secondary, respectively. 

In this manner we get 0 < y < 2x for the A° particles, and 0< gy <2 for 
0° particles, for which rt cannot be distinguished from 7. 

Similarly, we get 0 < 0*< x for A° particles (taking the proton trajectory 
as positive direction) and 0 < 0* < 2/2 for 0°. 

The distribution of other angles would also be of interest, but im our expe- 
riment we were unable to measure them. 

The presence of anisotropy in one or both distributions indicates first of 
all that the particles emerge at least partly polarized from the production event. 
Any asymmetries of Z(g) with respect to 90° and 180°, and of I (0*) with respect 
to 90° would be very interesting, in that they would indicate that one of the 
following assumptions is not correct: 1) that the particles have a determinate 
parity; 2) that there are no parity doublets; 3) that parity is conserved in 
the decay. Other forms of anisotropy can be studied better in distributions 
of g folded over with respect to 180° and 90° and of 0* folded over with 
respect to 90°. 

In the conditions of our experiment a number of circumstances might have 
attenuated or destroyed such anisotropy: 


1) The V°-events may have been produced either by pions or by nucleons 
and the primary energies may have varied over a wide range. 


2) V°-particles may have been generated by a secondary collision inside 
the nucleus and hence the primary track was not relevant to the elementary 
reaction producing the V° particle. 


3) We could not calculate the V° emission angle in the center of mass 
system of the colliding particle which produce the V°-event. It is not known 
how the polarization depends on this angle. 


4) Before leaving the nucleus in which it is produced, a V°-particle may 
experience nuclear collisions (4) and change its polarization. 


(13) See, e.g. G. MoRPURGO: Nuovo Cimento, 3, 1069 and 4, 1222 (1956); 5, 1785 (1957); 
S. B. TREIMAN and H. W. WyLp: Phys. Rev., 100, 879 (1955); S. B. TREIMAN, G. T. 
ReywnoLps and A. L. Hopson: Phys. Rev., 97, 244 (1955); R. K. ADAIR: Phys. Rev., 
100, 1540 (1955). 

(4) G. D. James and R. A. SALMERON: Phil. Mag., 46, 571 (1955). 
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In our conditions, observation of anisotropy is therefore rather difficult, 
but, provided it were significant, it would all the same yield useful infor- 
mation (15). 

Let us consider the distributions /(9) and (cos 0*) for A° and 0° in Fig. 10, 
11, 12 and 13. The category A° here again includes not only those individually 
identified, but also the « A° or 0° » events analysed as if they were all A° decays. 


6 
5 
3 
2 5 
i 4 
oe 60° 120° 180° 240° 300° 360° 3 
Fig. 10. — Distribution of the angle = 
for 33 A°%-events. The shaded area cor- 
responds to « A° or 0°» V°-events. 0° 60° 120° 180° 9 


Fig. 12. — Distribution of the angle @ 
for 21 0°-events. 


VEVENT 


= OMMI (6 =.2 cos 0"*.2 iG 1.0 0 ia .L cos 0°.6 8 1.0 


Fig. 11. — Distribution of cos 0* for 65 
A°-events. The shaded area corresponds 


Fig. 13. — Distribution of cos 0* for 42 
to « A° or 0°») V°O-events. 


0°-events. 


(5) J. BALLAM et al.: Phys. Rev., 97, 245 (1955); D. B. GAYTHER Phil. Mag., 45, 
570 (1954); 46, 467, 1362 (1955); F. ErsLer, R. PLano, A. PropeLL, N. Samios, M. 
SCHWARTZ, J. STEINBERGER, P. Bassi, V. BorELLI, G. PuPPI, H. TANAKA, P. Wa- 
LOSCHECK, V. ZoBoLI, M. Conversi, P. FRANZINI, I. MANNELLI, R. SANTANGELO, V. 
SILVESTRINI, G. L. BROWN, D. A. GLASER and C. Graves: Nuovo Cimento, 7, 222 
(1958); C. Graves and D. GLASER: circulating preprint. 
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1) I(g) for A°: this distribution suggests some slight anisotropy; we find 


N(g<180°) 21 È 
Nip 180°) A lago. 


N(g>180°) 12 


The folded distribution (0 < g < 90°) shows no other anisotropies. 


2) I(cos 0*) for A°: this distribution shows a crowding of events around 
9* = 90°. However, it is not a true crowding, since it is due: 


(i) to the presence of a considerable number of V°-events with very 
indeterminate cos 0* and hence always comprising the central part of the 
distribution, and 


(ii) to the contamination of 0°-events present among the « AS Ors G8). 
since one can see that, when analysed as A°, they give a contribution to small 
positive values of cos 0*. We may draw the conclusion that I(cos 0*) is con- 
sistent with isotropy. 

It may be of interest to consider separately the low momentum A° which 
should be emitted backwards in the center of mass system of the production 
reaction (1°). We have selected 11 A°events, with momentum less than 
500 MeV/c; 8 of them have cos 0*< 0, 3 have cos 0*> 0. If we add these 
numbers to the previous results given by the CERN Jungfraujoch group (1°) 
we get a total backward/forward ratio = 69/39. The probability that such 
a distribution is due to statistical fluctuation is 3%. 


3) I(g) for 0°: this histogram shows no indication of anisotropy. 


4) I(cos 0*) for 0°: this histogram apparently deviates from isotropy. 
If, however, we keep in mind that approximately 6 6° events have to be 
added in the interval 0.4 < cos 0* < 1.0 in order to allow for doubtful iden- 
tification, it will be seen that the deviation from isotropy is quite insignificant. 


7. — Analysis of the parent interactions. 


While scanning the photographs we have collected information on nuclear 
interactions with a view to learning something about the production of V°- 
particles. 

In the scanning we considered as useful nuclear interactions all those asso- 
ciated with V°-events and those without V° which satisfied the condition that 
at least two particles were capable of traversing a plate and that it was pos- 
sible to establish whether the primary particle was charged or neutral (the 
uncertain cases were less than 8%. Showers without V° and with less than 
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two penetrating prongs were excluded because, in comparing them with NE 
showers, we must somehow allow for the fact that there can be no V°-par- 
ticles in showers with an energy below the production threshold. We also 
excluded showers with more than 8 penetrating particles, because there is 
a great likelihood, in these cases, of not seeing any V°-particles which may 
be present. 

Since the error in our results is mainly due to the small number of V°- 
events available, we did not think that there was any point in extending the 
laborious scanning for non-V° showers to all the 30000 photographs. We 
limited ourselves to a group of photographs chosen at random. The number 
N_ of showers generated by pions was calculated from the number N, of charged 
primaries (n° and protons) and from the number N, of neutral primaries 
(neutrons), taking into account: (i) that equal numbers of protons and neutrons 
impinge on the shower-producing lead above the chamber; (ii) that conse- 
quently we have to expect below the lead an excess of neutrons over protons 
roughly equal to the excess in lead-nuclei, (iii) that the triggering system 
introduces bias against neutrons entering the chamber; this bias has been 
evaluated looking at the interactions with charged primaries and gives a loss 
of ~ 15% for neutral primaries. i 

The following results of this scanning may be of interest: 


1) The ratio between the number of A°showers with a primary n° and 
the number of all showers, whether or not accompanied by V°-events, with 
a primary «+, was found to be 


ee Se ae 


similarly for 0° we have N%Z/NZ = 0.06 (1 + 0.4). 


2) The ratio between the number of A°-showers with nucleon primary 
and the number of all showers with nucleon primary, was found to be 


Tn 
Ao x 
N° = 0.03 (1 Ss 0.6) , 
similarly for 0° we have 
ADs tog (1+0.6 
nue 0 LE DO) 


These production frequencies are corrected for unobservable decays and geo- 
metrical biases. 
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These results. indicate once more (2) that the V°-particles belonging to 
showers triggering our chamber (see the histograms of Fig. 14, 15), are more 
often than not generated by z+. 


i 


1000 


ber 


n 
© 
=} 
Tf 


Extrapolated num 


ap 


Fig. 14. — Distribution of showers without V°-events vs. multiplicity m (1355 showers 
scanned). 


3) The two histograms of Fig. 14, 15: The histogram of Fig. 14 re- 
presents the number of showers of multiplicity m, defined as the number of 
particles penetrating at least one plate, that of Fig. 15 represents the number 


30 


Weighted number 


20 


0 1 2 9 4 5 6 if TSE ia 


Fig. 15. — Distribution of showers with V°-events vs. multiplicity m (number of showers 
considered: 48). 


(16) R. B. LEIGHTON, S. D. Wantass and C. D. ANDERSON: Phys. Rev., 89, 148 
(1953); W. B. FRETTER, M. M. May and M. P. NAKADA: Phys. Rev., 89, 168 (1953); 
H. 8. BRIDGE, C. Peyrou, B. Rossi and S. SAFFORD : Phys. Rev., 91, 362 (1953). 
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of V°-particles produced in showers of multiplicity m. These two histograms 
are normalized to the same number of pictures. Since many V°-particles be- 
longing to zero-multiplicity showers and having a neutral primary could not 
be included in this sample, the ordinate for m = 0 is open to doubt. We have 
not distinguished between A° and 9° because from the analysis of the asso- 
ciated interactions no significant differences have arisen. 


A comparison of the two histograms agrees with the observation of other 


authors (17) namely that V°-production increases with the energy released in 
the shower. 


We are very grateful to Dr. J. P. AstBURY, Dr. J. A. NEWTH and the 
theoretical group in Rome for helpful discussions; to Dr. A. DE Marco for 
help in running the experimental apparatus and in the analysis of the events; 
to Prof. E. AMALDI for his interest and encouragement. One of us (G.A.) 
expresses his gratitude to Prof. E. AMALDI for the facilities extended to him 
in the Institute and to Ing. A. MAYER for a financial grant. 


(7) H. S. BRIDGE et al.: loc. cit. 


RIASSUNTO (*) 


Si discutono i risultati di un’esperienza su aleune proprietà delle particelle A° e 6° 
prodotte nei tetti di una camera a nebbia. Per mezzo di misure angolari e stime della 
ionizzazione si analizzano 115 eventi V®. Si trova che la frequenza di produzione di 
coppie A°, 00 è bassa. Per nessun tipo di particella si trova alcuna prova di anisotropia 
nei decadimenti. La sezione d’urto per la produzione di VO risulta maggiore per un 
pione primario che per un nucleone e aumenta con l’energia liberata nello sciame. 


(*) Traduzione a cura della Redazione. 
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Sommaire. — Il existe en relativité générale deux méthodes d’approxi- 
mation pour déduire les équations du mouvement d’un système de corps 
à partir des équations du champ. L’une utilise une représentation de la 
matière par les singularités du champ extérieur; l’autre fait jouer un role 
au tenseur d’impulsion-énergie. La liaison entre ces deux points de vue 
peut ètre établie en observant que les discontinuités du tenseur d’Ein- 
stein sont représentées par le tenseur d’impulsion-énergie. Utilisant alors 
le théorème de Stokes, nous montrons quo les deux méthodes sont équi- 
valentes. 


I. - Rappels géométriques. 


1. — La variété espace-temps. 


Dans la théorie de la relativité générale, l'élément primitif est constitué 
par la variété espace-temps. C’est une rariété différentiable V, à quatre 
dimensions, de classe de différentiabilité (e?, ct par morceaux) sur laquelle est 
définie une métrique riemannienne partout de type hyperbolique normal, à 
un carré positif et trois carré négatifs : 


ds? = 9,,,(x") da” da? (x, B et tout indice grec = 0, 1, 2, 3). 


Les dix coefficient g, a(t") sont appelés potentiels de gravitation pour le systeme 
de coordonnées locales (2°). 

L’équation ds? = 0 définit en chaque point x de V, un còne réel C,, le cone 
élémentaire en x. Une direction est dite orientée dans le temps ou orientée 


a 
q 
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dans l’espace selon qu’elle est intérieure ou extérieure au cône C,. Pour qu’une 
hypersurface à trois dimensions, définie localement par l'équation LIU; 
soit orientée dans l’espace, il faut et il suffit que: 


© 
A,f “i g°4 of Og 105 (2, | = 


~ da 


Un système de coordonnées sera dit physiquement admissible ou inter- 
prétable dans l’espace et dans le temps si: 


1) les hypersurfaces W, définies par 4° = const. sont orientées dans l’es- 
pace (ce qui entraine 9° > 0); 


2) les lignes Z sur lesquelles la variable x° varie seule sont orientées 
dans le temps (alors go > 0). 


Dans un tel système de coordonnées nous dirons que æ° est une variable 
temporelle, «1, a, x définissant les variables spatiales correspondantes. Les 
varietés W, seront appelées sections d’espace et les lignes L, lignes de temps. 

Nous supposerons dans la suite que l’espace-temps Vi, a été rapporté a 
un système de coordonnées physiquement admissible, caractérisé par l’en- 
semble des deux propriétés go) > 0 et y°°> 0, et nous adopterons sur les sec- 
tions d’espace W; la métrique ds*? intrinsèquement définie, c’est-à-dire inva- 
riante par tout changement de systeme de sections d’espace: 


o Jo Jos ; : a) ese F 
RE (a -- malte) a dx’ (i, j et tout indice latin = 1, 2, 3). 


\ 00 


2. — Le système des équations d’Einstein. 


Les propriétés géométriques de l’espace-temps introduisent uniquement les 
phénomènes de gravitation. Le tenseur métrique g, 8 est astreint à vérifier le 
système des dix équations d’Einstein: 


Shap = jae, 


qui généralisent les équations de Laplace-Poisson qui sont A la base de la 
théorie du potentiel. S,, et 7; sont deux tenseurs symétriques, x désigne 
un facteur constant lié à la constante d'attraction universelle. 

Le tenseur 7°; dont l'interprétation est purement physique doit décrire au 
mieux, en chaque point de l’espace-temps V,, l’état de la distribution éner- 
gétique (cas intérieur). Dans les régions de V, non balayées par l’énergie T,, 
doit étre identiquement nul (cas extérieur) 
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Le tenseur S,,, d'origine géométrique, qui ne dépend que des I. et de 
leurs dérivées des deux premiers ordres, est linéaire par rapport aux dérivées 
du second ordre et satisfait aux identités de conservation: 


VA Sa (V; opérateur de dérivation covariante). 
Ce tenseur s'exprime en fonction du tenseur de Ricci R,, de la variété rieman- 
nienne par la formule: 


Sip = Rap — 4918 (R + 22) (A, constante cosmologique). 


II. - La déduction des équations du mouvement 
à partir des équations du champ. 


3. — Généralités et hypothèses. 


La théorie de la relativité générale est une théorie du champ pur en ce 
seus que l’on peut déduire les équations du mouvement à partir des équations 
du champ seules. Cette déduction est basée sur l’approximation des équations 
du champ. Ayant calculé une solution approchée de celles-ci, on définit les 
équations du mouvement comme des conditions que cette solution doit satis- 
faire en raison même des équations du champ. Nous verrons dans les deux 
paragraphes suivants quelles conditions il conviendra de prendre suivant que 
l'on veut utiliser les équations d’Einstein du cas extérieur ou celles du cas 
intérieur. Il existe en effet deux techniques différentes pour traiter par ap- 
proximations le problème du mouvement de N masses gravitantes. L’une 
utilise une représentation des masses par de pures singularités du champ ex- 
térieur, l’autre fait jouer un rôle au tenseur d’impulsion-éenergie 7° ;. Elles 
sont connues respectivement sous les noms de méthode des singularités (12%?) 
et de méthode du tenseur d’impulsion-énergie (45:78), 


(a 
Li 
209 ( 
(s 


) A. EINSTEIN, L. INFELD et B. Horrmann: Ann. Math., 39, 65 (1938). 

) A. EINSTEIN et L. INFELD: Ann. Math., 41, 455 (1940); Can Journ. Math., 1, 
1949). 

) L. InreLp: Acta Phys. Pol., 16, 177 (1957). 

(4) V. A. Fock: Journ. Phys. Acad. Sci. URSS, 1, 81 (1939). 

(5) M.me F. HENNEQUIN: Etude easy bai des approximations en relativité géné- 
rale et en théorie unitaire de Jordan-Thiry, « Bull. scient. Com. Trav. hist. et scient.» , 
1 (1956), 2% partie: Mathématiques (Paris, 1957), p. 73-154 (Thèse Se. math. Paris, 1956). 

(6) A. LicaneROWICZ: Les théories relativistes de la gravitation et de Vélectromagne- 
tisme (Paris, 1955). 

(7) A. Papaprerrou: I: Proc. Phys. Soc., A 64, 57 (1951); II: Proc. Phys. Soc., 
A 64, 302 (1951). , | 

(8) N. M. Petrova: Akad. Nauk SSSR, Zu. Eksper. Teoret. Fiz., 19, 989 (1949) 
(en russe). 
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La méthode des singularités peut étre présentée d’une manière autonome 
car il n’est pas nécessaire, du moins théoriquement, de préciser la nature des 
singularités du champ. Mais une comparaison entre les deux méthodes n’a 
évidemment de sens que si les divers champs intérieurs et le champ exterieur 
sont créés par le même ensemble de distributions énergétiques dans l’espace- 
temps. On doit done envisager des métriques globalement réguliéres; la notion 
de modèle d’univers (6) que nous allons rappeler nous en donnera une idée 
précise. 

Un modèle d’univers est une variété V, munie d’une métrique partout ré- 
gulière satisfaisant aux conditions: 


1) Dans les domaines de V, balayés par une distribution énergétique et 
limités par des hypersurfaces frontières S, la métrique est régulière et satisfait 
aux équations d’Hinstein du cas intérieur. 


2) Dans les domaines de V, qui ne sont pas balayés par aucune 
distribution énergétique, la métrique est régulière et satisfait aux équation 
d’Einstein du cas extérieur. 


3) A la traversée d’une hypersurface S, les potentiels et leurs dérivées. 
premières sont continus conformément aux conditions de raccordement. 


Quand il est possible de construire un tel modèle d’univers, le champ ex- 
térieur peut être considéré comme effectivement produit par les différentes 
masses ou distributions énergétiques en mouvement et c’est le raccordement 
des champs intérieurs des différentes distributions avec un même champ qui 
assure Vinterdépendence des mouvements (5). 

Les hypothèses que nous faisons dans le problème du mouvement sont les 
suivantes : 


a) Nous considérons un modèle d’univers à N tubes massiques repré- 
sentés dans chaque section d’expace W, par N domaines à trois dimensions 


s 
B* (s=1,..., N) connexes et bornés, sans point commun, qui correspondent à 
des corps. 


b) La métrique est quasi-euclidienne et à comportement asymptotique 
euclidien. Les potentiels admettent des développements limités selon les puis- 
sances de c-? (c, vitesse de la lumière), dont les premiers termes sont les po- 
tentiels euclidiens (+1, —1, —1, —1) de Minkowski. 


c) Nous supposons également que les vitesses des corps sont faibles par 
rapport à celle de la lumière. La dérivée partielle 0, est alors d’un ordre infi- 
nitésimal supérieur aux dérivées partielles @,: 


d = 0(1/e) . 


n 
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4. — Identités de conservation et équations du mouvement (*). 


41. Identités de conservation. — La représentation de la matière peut s’ef- 
fectuer en relativité générale par des schémas du type hydrodynamique. Un 
fluide relativiste est alors considéré comme un milieu matériel continu, défor- 
mable, doué des propriétés suivantes: 


1) il posséde une densité propre 0; 
2) il est possible de définir un vecteur vitesse unitaire w pour chaque 
point du milieu. 


Le tenseur d’impulsion-énergie 7, est destiné à assurer une représen- 
tation aussi complète que possible d’une distribution énergétique déterminée. 
Ce tenseur contiendra ainsi différents termes correspondant aux différentes 
sortes d'énergie de la distribution et à leur interaction. La description de ces 
différents termes, leur importance relative, sont fournies par l’expérience. 

Soit T; le tenseur d'énergie du schéma matériel considéré. Le tenseur #4 
étant conservatif, il en est de même du tenseur T,, qui est astreint, d'après 
les équations d’Einstein, à n’en différer que par un facteur constant: 


(4.1) V, 0% =0. 


Ces relations sont appelées identités de conservation du schéma considéré. 

On déduit de (4.1) un certain nombre de propriétés de «conservation » pout 
les éléments physiques qui figurent dans l’expression du tenseur d’énergie 
considéré (masses, pressions, vitesses, champ électromagnétique, champ des 
températures propres). Examinons les relations (4.1) pour les tenseurs as- 
sociés aux schémas usuels. 


1) Schéma matière pure. — Le terme le plus important dans 7°; lorsqu'il 
existe correspond toujours à l'énergie pondérale. Si l’on néglige toutes les 
autres formes d'énergie, on aura alors le schéma matière pure avec: 

D) == f 
(4.2) Tg = QU,Ug- 


Les identités de conservation s’écrivent: 


(4.3) Ve re — V, (ouf) “w+ ouf Vu, #0 


(*) Cf. M.me F. HENNEQUIN (°) 
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Formons la combinaison linéaire w'(V, TS) et tenant compte de wu,=1, il 
vient: 

(4.4) UV, T5) = V,(0 u*) = 0. 

C'est l'équation de continuité du milieu. Les équations (4.3) se réduisent 
alors a: 


5) WON OS 


ra 
DI 


Le système (4.5) exprime que les lignes de courant, définies comme trajectoires 
du champ de vecteurs u, sont géodésiques de la variété riemannienne. 

Ce qu'on nomme le principe des géodésiques apparait comme une consé- 
quence de (4.1) et des conditions de raccordement de Schwarzschild. Si l’on 
considère une masse d’épreuve dans un champ de gravitation extérieur donné, 
la nécessité d’un raccordement entre champ intérieur et champ extérieur fait 
que la masse d’épreuve doit décrire une géodésique du ds? extérieur. 


2) Schéma fluide parfait. — On appelle ainsi un état de la matière dont 
la description peut être faite à l’aide des grandeurs u, 0, p, les deux grandeurs 


scalaires étant supposées liées par une relation fonctionnelle dite équation 
d'état o — g(p). Le tenseur d’énergie associé est de la forme: 


p D 
(4.6) T8 = (0 — | Uy, Up — È CAE 


Ecrivons les identités de conservation pour ce tenseur: 


(4.7) We = [ 0% Ae + (0 - a) UNU VS I 
i 10 LL e x e? pa % 2 
Alors 
(4.8) (Vs Th = V,|(0 +E) 0 
$ 0 e? x C2 È) 
et les conditions de conservation (4.7) deviennent 
A (PIE) ‘5 
4°9) PIAVE ip, > AI RA rio 
|: eee) ee ae 


Ces équations définissent les géodésiques d’une métrique conforme à la métrique 


+ 
~ 
N 


COMPARAISON DES DEUX MÉTHODES D'OBTENTION DES EQUATIONS ETC. 653 


riemannienne ds? = F? ds? où la fonction /, appelée indice du fluide, est définie 
par 


42. Equations spatio-temporelles du mouvement. — Nous voyons par ce qui 
précède que les identités de conservation ont un contenu physique important. 
En particulier elles fournissent le système différentiel aux lignes de courant. 

Nous donnerons done la définition suivante: 

Les équations spatio-temporelles du mouvement d'une masse Wépreuve placée 
dans un champ sont fournies par les identités de conservation du tenseur d'impul- 
sion-énergie représentant les distributions énergétiques du champ. 

Il revient évidemment au même de considérer des combinaison linéaires 
indépendantes de ces relations. Les quatres identités de conservation du 
tenseur d’impulsion-énergie correspondent à trois équations du mouvement 
au sens classique Sol 0) et a une équation de continuite u*(V T°) = 0. 


1) Equations du mouvement au sens classique. — L’introduction des den- 
sités tensorielles permet de mettre les identités de conservation sous une forme 
intéressante, notamment du point de vue des calculs d’approximation. Nous 
serons amenés alors à multiplier ces identités par VW —ÿ}, et nous écrirons les 


trois équations spatio-temporelles, généralisant les équations classiques du 
mouvement sous la forme: 


(4.10) V=gVy Tî = OT V—g) — HT V— 9) 0: Gap = 0. 


Les équations spatio-temporelles du mouvement d’un corps seront obtenues en in- 
tégrant les équations spatio-temporelles d'une masse d’épreuve sur le volume d'espace 
occupé par le corps. En mettant en évidence l'élément de volume de W, 
(V— g* da) dæ?/\daæ®) les équations du mouvement d'un corps B? sont donc: 


> 


(E.1) [VE viv or aor dat des = 0. 


B3 
di . / 4 
Nous poserons dans la suite w= Vg/9*. 


2) Equation de continuité. — Multiplions l'équation (4.4) par V—g et 
écrivons le premier membre de l’équation obtenue sous la forme d’une diver- 
gence ordinaire. Posons d’autre part 
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ou v' = da'/dt est la vitesse d’espace et w°— 1. Il vient 
(4.11) Ve gV (ou) = 0, (ow? ADE gv’) = () 3 


La comparaison avec l’écriture classique des conditions de continuité nous 
conduit & choisir comme densité de masse la quantité 


(4.12) m= ou V— g . 


Intégrons l’équations w(V, T°)V—g = 0 sur un domaine D, à quatre di- 
mensions limité par un tube de courant 7? et deux sections d’expace B? (x° = ct) 
eu be (4° — ct’), LI vient 


(4.13) | UV, 02) V— g da? À da A da? A da = 0 : 


. 


Di 


En transformant (4.13) en une intégrale étendue à la frontière de D,, on déduit 
les propriétés suivantes relatives aux schémas usuels: 

Dans le cas d’un schéma matière pure, la masse M = in dat dx? dx? 
est constante au cours du temps. a 

Dans le cas d’un schéma fluide parfait, la variation de la masse M au cours 
du temps est uniquement liée aux termes de pression p. Le choix de m nous 
donne donc la masse la plus conservative associée au schéma fluide parfait (5). 


43. Rôle des conditions d’isothermie. — Nous avons à rechercher une solution 
approchée des équations du champ et à la substituer dans les équations du 
mouvement. Les dix équations du champ ne sont pas indépendantes, leurs 
premiers membres satisfont a quatre identités de conservation. On peut done 
fixer localement quatre potentiels, ou plus généralement se donner quatre 
relations de coordonnées. Pour avoir des résultats qui puissent étre interpré- 
tables physiquement, il nous faut choisir parmi l’ensemble des systèmes de 
coordonnées un système qui soit adapté le mieux possible à la réalité physique. 
A ce point de vue, les coordonnées qui s'imposent sont les coordonnées iso- 
thermes (G. DARMOIS). 

Dans un système isotherme, SŸ se réduit à une expression que nous dé- 
signons par $$. Dans chaque équation S'5,= ~T*’ ne figurent que les dérivées 
secondes de g*"/—y. Cette propriété nous permet de faire usage d’un pro- 
cédé en cascade d’approximations. Quand nous aurons à notre disposition 
une solution approchée, nous appellerons équations approchées du mouvement 
ce que l’on obtient en substituant la solution approchée dans les équations 
du mouvement. 


~ 
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Qt 


Cependant, notre but initial était la résolution de: 
SR PTE 


En choisissant un système de coordonnées isothermes, nous nous sommes 
ramenés au système plus simple 


+B B 
SOs, = ve >) 


pour lesquelles nous avons écrit les équations du mouvement sous la forme 
habituelle V;7%= 0; on peut se demander si les conditions d’isothermie sont 
vérifiées pour la solution trouvée, c’est-à-dire si le problème premier a bien 
été traité dans sa généralité. 

La réponse à cette question a été apportée par Mme F. HENNEQUIN. L’auteur 
a étudié le rôle des conditions d’isothermie et démontré que les équations du 
mouvement approchée d’ordre p entrainent les conditions d’isotermie appro- 
chées d'ordre p (°). 

Il y a ainsi équivalence des conditions d’isothermie et des équations du 
mouvement relatives à une solution approchée. Certains auteurs (N.M. PE- 
TROVA, V. A. Fock) ont utilisé effectivement les conditions d’isothermie pour 
former les équations du mouvement; A. PAPAPETROU et Mme F. HENNEQUIN 
ont utilisé les conditions de conservation (*). 


5. — Flux et équations du mouvement. 


La déduction des équations du mouvement pour un systéme de corps a 
partir des équations du champ extérieur: 


(5.1) see. 


repose essentiellement sur les propriétés que possède le flux du vecteur 
C= wS,,. de composantes wSî, (x, indice fixé quelconque) défini sur une 
section d'espace W,. 

= gig* définit un scalaire sur W,. Ce facteur est mis dans Sy = VS, 
simplement pour faciliter la comparaison des équations (E.1) (p. 653) et (E.2) 


(*) La dépendance entre les équations du mouvement et les conditions d’isothermie 

subsiste pour le champ extérieur, d’après la continuité des g,; et 029xp (les conditions 
5: ‘ a Bp SRI : 

d’isothermie sont définies par — gm 0) et d’après les résultats obtenus dans la com- 


paraison des deux méthodes (voir plus loin). 
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(p. 659). Son emploi n’est pas indispensable dans l’étude des équations du mouve- 
ment par la méthode des singularités (*). Naturellement, on aurait pu aussi définir 
les équations (E.1) sans le facteur w et utiliser alors le flux du vecteur $,, 

Considérons dans W, un champ d’intégration à trois dimensions contenant 
un corps B® déterminé. Soit 0° ce champ d’intégration et soit CC® sa frontière 
Désignons par ©,, la forme différentielle extérieure d’ordre 2 associée au 
vecteur S,,,: 


(5.2) Qe = EN in (SE) dat AE , 


7*,, est le tenseur complètement antisymétrique attaché à la forme élément 
de volume de W,. Par définition le flux du vecteur $,, à travers la frontière 
de 03 est l’intégrale de la 2-forme Q,, sur le 2-champ @C3: 

~ ro < a 

(5.3) flux 5e So = | 2 5 


lo] lol 3 


1) Soit dQ,,, la différentielle extérieure de la forme ©, ,: 
(5.4) dQ = teind(S? V—g) Ada Ada = d(StV— g) dt A de? A das. 
Etudions la condition dQ,,,=0 à l’aide des identités de conservation 


il 


CARS TAO) 
(5.5) VS = 
MSI 


9;(8° V— 9) De 480, au — 0. 


Ces identités peuvent s’écrire encore, après multiplication par V— g: 
Re k 0 oA 
(5.6) Or SE V— g) I (Se V— q) nd? LV— 9 S "Oy Gay =0. 


Par conséquent la condition 0,(S'%/—g)=0, c’est-à-dire di. = 0, ne pews 
pas être satisfaite identiquement. Cependant la dérivée partielle è, est d’ordre 
1/e, de même que les Jan (*). On en déduit que la quantité dQ. d'ordre 1 


est nulle si les équations du champ d'ordre (J — 1) sont vérifiées, ce que nous 
traduirons par 


(5.7) GOR iim 


(*) T. H. PHam: La méthode des singularités pour les équations du mouvement en 
relativité générale et en théorie du champ unifié (Thèse Sc. math., Paris, 1957) (multi- 
graphiée); Compt. Rend. Ac. Sci. Paris, 246, 61 (1958); voir aussi dans: Réunion 
d’études et de mises au point tenue sous la présidence de L. pp BROGLIE. Paris.1958: 
Problèmes actuels en théorie de la relativité (à paraître). N 


(*) Avec les développements habituels des potentiels en coordonnées isothermes 
les 0,9, Son d’ordre 1/c? au moins. 


7 
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la notation | 
approximations antérieures. 

Considérons deux champs d’intégration C” et OC” dont l’un est intérieur 
à Pautre. La forme ©, étant régulière dans le domaine (0° — 0") compris 
entre les frontières 00” et QC”, il vient par application du théorème de Stokes: 


- signifiant que l'égalité a lieu à chaque ordre / en vertu. des 


= = . 


(5.8) Ja 2x sa fa, ari Ja. 


. 
dc 03) c'3- 0"3 


> 

a 
a 

©) 
è 


Les relations (5.7) et (5.8) entrainent que le flux du vecteur S.,, ne dépend pas 
du choix de C*. En particulier on peut évaluer ce flux en déterminant sa valeur 


sur la frontière du corps contenu dans 03: 


(5.9) 10e = [2x C 


J x 
00° OB 


2) Examinons maintenant la manière dont les potentiels des divers 
ordres interviennent dans l'intégrale (5.3). 

Pour calculer cette intégrale nous utiliserons un espace euclidien de repré- 
sentation €; à trois dimensions. Désignons par A le domaine représentatif 
de 03, et par 0A sa frontière que nous rapportons à des paramètres (t, ¢?). 
Sur 003 la forme £Q,,, admet une forme induite Q,,, qui peut s'exprimer à l’aide. 
des variables (¢!, 1?) par la formule 


(5.10) 0 =) dn da 
ou 

n DERE, 405) 
~ fs 5 QE 7! 
(5.11) eR = EN iin (Ul Se Dit, €) c 


L'intégrale de la forme 2,,, sur le champ 00° n’est autre que l’intégrale de la 
fonction f,, sur le champ 04. On sait qu’elle peut se mettre sous la forme 


que TT A 
(>. 12) 12 = Jos | di A dé?, 
00% 1 
où: (v,) = normale unitaire à 005, 
gî = déterminant de la forme quadratique fondamentale de Cl? rap- 


porté aux variables (tt, #2) (*). 


(*) Voir A. LiCHNEROWICZ: Algèbre et analyse linéaires (Paris, 1947), p. 192-194. 
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Explicitons maintenant l'expression S* qui figure sous le signe d’integration. 
Posons: 


(5.13) PS fy Ho Ot (8 0 Vg) 


Par un caleul classique on obtient la formule 


(5.14) S xp} == Ca (g'# anes —s g°° OL Lee gl 0, E° DE ge Oo F°) Ae Hae 
2V— 9 


HŸ désignant une fonction des potentiels et de leurs dérivées premières. De 
(5.14) on tire: 


(5.15) SESSI 
ou 
(5.16) DIO SIG Ce 00) 


est l’ensemble des termes linéaires qui ne comportent aucune dérivation par rap- 
port à la variable 2°; A" comprend les autres termes linéaires ainsi que les termes 
non linéaires de S*. 

Nous avons séparé les terms de S* en deux groupes. La raison en est la 
suivante: si l’on considère l’expression S° d’ordre 1 quelconque, seul ®° contient 
les potentiels d'ordre 1 tandis que A* ne fait intervenir que les potentiels d'ordre 
au plus égal à ([—1). Ceci est une conséquence de l’hypothèse c). 

On constate que @* est de la forme 0,F", où F% est antisymétrique par 
rapport aux indices k,/. On peut alors montrer que dans l’espace euclidien 
de répresentation €; rapporté a des coordonnées rectilignes, la partie de l’in- 
tégrale (5.12) qui porte sur @* peut se ramener à l'intégrale sur 04 d’une 
2-forme homologue a zèro. Le champ d'intégration ¢4 ayant une frontière 
nulle, une telle intégrale est toujours identiquement nulle (+). Ainsi, seul A” 
contribue à l'intégrale (5.12), ce qui montre que le flux de se d'ordre l ne 
dépend effectivement que des potentiels d'ordre (1— 1). Il résulte que le calcul 
du flux d'ordre / n’est lié qu'à la solution des équations du champ d’ordre 
({— 1). Cette propriété importante signifie que si l’on substitue aux g° une 
solution des équations N,,— 0 le flux de SP ne sera pas identiquement nul. 


te U : . Se 72 È 8 E 

(*) Sur l'emploi de la densité tensorielle g comme variable de champ, cf. (?) 
2° référence. 

(*) Cf. (F2) (cas de la variété V,) et (°) (cas d’une variété V, à n dimensions). 


. 
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Or les équations du champ S,,— 0 impliquent nécessairement 


éc3 


En raison des propriétés du flux de ST. ces conditions ne sont pas vides. 
Elles constituent des relations que doit véritier toute solution approchée des 
équations du champ, relations qui ne dependent que de la frontière du corps 
contenu dans 08. Nous pouvons done les prendre comme définition des équa- 
tion spatio-temporelles du mouvement. 

De facon plus précise, les équations du mouvement au sens classique s’ob- 
tiennent à partir de 


(E.2) fes = [20 0! 


003 êB3 


La condition (5.17) pour «=0, ou plutôt une combinaison analogue à (4.13) 


(5.18 Q=0, Q= u® S2) da” Ada" A dx”, 
à) 


31 Mauro 
OD, 

renseignera sur la variation de la masse (*). L’intégrale (5.18) est étendue a 

la frontière d’un domaine D, à quatre dimensions limité par un tube de cournat 

et deux sections d’espace, 17,,,, désigne le tenseur complètement antisymétrique 

attaché à la forme élément de volume V4. 


TIT. - Comparaison des deux méthodes d’ obtention 
des équations du mouvement. 


6. — Généralisation de la méthode utilisant le flux. 


Nous nous proposons dans cette section d'établir le rapport entre la méthode 
du tenseur d’impulsion-énergie et la méthode des singularités. Nous utili- 

(*) Dans les applications, la méthode des singularités a toujours été envisagée 
comme un procédé autonome. La condition : 


joue de ce fait un rôle tout autre. Elle intervient alors dans l’approximation des équa- 
tions du champ; associée avec les conditions de coordonnées elle permet de déterminer 
les fonctions harmoniques qui s'introduisent dans la solution des équations du champ 


extérieur. 
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serons principalement les relations entre le tenseur d’impulsion-énergie et les 
discontinuités du tenseur d’Einstein. Ces discontinuités se produisent pendant 
la traversée d’une hypersurface frontière S limitant une domaine de Va. 
balayé par une distribution énergétique. D'autre part, pour la clarté de l'exposé 
nous considérerons des équations intermédiaires qui se prêtent bien à la com- 
paraison avec les équations du mouvement (E1) et (E2). Nous obtiendrons 
ces équations en généralisant la méthode utilisant le flux au cas des équations 
du champ intérieur. 

Ecrivons en effet les équations d’Einstein du cas intérieur sous la forme: 


(6.1) Se ag AT ap = 0. 


= 


Introduisons le vecteur d’espace &,, qui admet les composantes contra- 


variantes Wes et soit 
(6.2) Qu MW ae del 


la forme différentielle extérieure d’ordre 2 associée à ce vecteur. Considérons 
les équations suivantes 


(6.3) Jo 0 


qui expriment que le flux du vecteur s,,, à travers la frontière du corps B* 
est nul pour toute solution des équations (6.1). 

Nous raisonnerons sur les équations (6.3) en vue de montrer l’équivalence 
complète des deux méthodes. Les résultats obtenus s’appliqueront a fortiori. 
aux équations 


(E.3) Jo AK 


oB* 


qui correspondent aux équations du mouvement au sens classique. 


7. — Relation entre les équations (E.1) et (E.3). 


La métrique satisfaisant aux équations (6.1) est regulière dans B®. Nous 
pouvons appliquer le theoreme de Stokes a Vintégrale de surface (6.3) et la 
transformer en une intégrale de volume étendue au domaine B?: 


(7.1) | On | do, , 


BS B3 
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où do, dèsigne la différentielle extérieure de la forme @,,): 
CR) do. = 36% dive g) A dx \ dai = dx(stV— g) dat A dr? \ da. 
D'autre part, le tenseur s,, qui est symétrique vérifie la formule 


/— gVysb = OV — 9) —$ V—9 8" 0,9 95 ) 


— 
Can] 
C9 


que Von peut écrire encore 

At pag Babich 3 ,B A 07 A1 300 1 

(7.4) Ast V— g) = V—9Vpsh — [dsx V— 9) — 4 V— 98° Jos] - 
Il vient done 


(7.5) Jor = Re gVs sf det \ da? \ da? — 


ep? 33 
= | [ao(s° V— g) — + V— 98° 0, Foal dat À dr? À dat. 
B 
Examinons l’intégrale (7.5) d’ordre /. Le deuxième terme du second membre 
qui ne fait intervenir que les $,; d'ordre (J —1) est nul parce que les équa- 
tions du champ d’ordre (/ — 1) sont vérifiées. Il en résulte que cette intégrale 
se réduit au premier terme du second membre qui est égal a 


ee) | Dea” ne gV(Sî — xT%) dat \ da? \ dar? = 
eB? B3 | 
Ta HC VEZ IV, Te aa /\ da? A dx = 


Bs 


puisque le tenseur d'Einstein est conservatif. Rappelons que !_ signifie: éga- 
lité à chaque ordre / en vertu des approximations antérieures. 

Ainsi, on voit que les équations du mouvement (E.1) peuvent être remplacées 
par les équations (E.3) qui en constituent une autre formulation. 

Nous noterons que le tenseur S,, a disparu des équations (7.6) gràce aux 
identités de conservation. Dans les équations du mouvement ne figure done 
plus que le tenseur T,,. De là vient la simplicité relative des caleuls d’ap- 


proximation auxquels conduit l'emploi du champ intérieur. 


8. — Relation entre les équations (E.2) et (E.3). 


La comparaison de ces équations est immédiate. Affectons respective- 
ment du signe + et du signe — les expressions de Sp à l’intérieur et à l’exté- 
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rieur de B3 et désignons par 


[Se] == Sy TA Srp 


È È se it fe DEAD, I TUR tes 
la discontinuité subie par le tenseur d’Einstein à la traversée de cB*. Sur cB 


on a l’égalité 


— 
(0 0) 
Ha 
De 


AN GE IEEE 


En portant (8.1) dans l’expression de ©,, on voit que l’intégrale (6.3) se met. 


sous la forme 


(8.2) fo = [2 a frate — T°) de’ A da’. 


… 


ass aps 0B3 


Mais la discontinuité du tenseur d’ Einstein est précisément représentée par y7,,:. 


iS = 4 Loe ) 


et la deuxiéme integrale de (8.2) est nulle. Done 


lo = fe ) 


ape aps 


ce qui montre en particulier que les équations (E.3) sont équivalentes aux équa- 


tions du mouvement (E.2). 


9. — Equivalence des deux méthodes. 


En raisonnant sur les équations (6.3) nous obtenons non seulement l’équi-. 


valence de (E.1) et (E.2) mais aussi celle de (4.13) et (5.18). Par un calcul ana- 
logue et en tenant compte de l’ordre de grandeur des 0,u*, nous pouvons d’ail- 
leurs démontrer directement que les équations (4.13) et (5.18) se ramenent 
l’une à l’autre, Ainsi se trouve établie l’équivalence entre la méthode du 
tenseur d’impulsion-énergie et la méthode des singularités. 

Nous pouvons énoncer: 


THÉORÈME. — Dans la théorie de la relativité générale les équations spatio- 


temporelles du mouvement de N masses gravitantes peuvent étre définies par les: 
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4N conditions (6.3): 


al 
Voi nix ws, dat \dx', 


suivant que Von considère le champ intérieur ou le champ extérieur. 
Les conditions 
fo = 0 


ôB* 


correspondent aux équations du mouvement au sens classique tandis que la com- 
binaison linéaire 


È DZ œ 2 
| wo=0, D = 37 Mayo Suda” \ da? A dx”, 


OD, 


joue le rôle d’une équation de continuité généralisée. 


RIASS TN TOM) 


Esistono in relatività generale due metodi d’approssimazione per dedurre le equa- 
zioni del moto d’un sistema di corpi partendo dalle equazioni del campo. Uno utilizza 
una rappresentazione della materia per mezzo delle singolarità del campo esterno; 
l’altro fa intervenire il tensore energia-quantità di moto. Si può stabilire il legame tra 
questi due punti di vista osservando che le discontinuità del tensore di Einstein sono 
rappresentate dal tensore energia-quantità di moto. Allora, utilizzando il teorema di 
Stokes, si dimostra che i due metodi sono equivalenti. 


(*) Traduzione a cura della Redazione. 
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On the Mechanism of the (n, 2n) Reaction 
at 14 MeV: Neutron Energy. 


E. Remy (*) and K. WINTER (*) 


Laboratoire de Physique Atomique et Moléculaire du College de France - Paris 


(ricevuto il 28 Aprile 1958) 


Summary. — Neutron spectra emitted by °Be, Mn, Fe, 1271 and ÆITa 
under bombardment of 14 MeV neutrons have been measured by using 
time-of-flight techniques. Continua of the first neutrons are separated 
from low energy groups due to the second neutrons. Successive emission 
is predominant for medium weight and heavy nuclei. A possible con- 
tribution of simultaneous emission by light nuclei cannot be decided. 
The conclusion is used to interprete inelastic scattering data from the 
interaction of 14 MeV neutrons. 


1. — Introduction. 


There now exist some investigations of the process of the (n, 2n) reaction 
on ‘Be (17). However, no experiments on the mechanism of this reaction are 
known for heavier elements. The knowledge of this mechanism and of the 
Subsequent energy spectrum of the emitted neutrons is important for the 
interpretation of inelastic scattering data from the interaction of 14 MeV 
neutrons in the aspects of the statistical model. 

HUBER and WAGNER (') investigated the energy spectrum of neutrons 
emitted by ‘Be under bombardment of 3.7 MeV neutrons. Using photographic 
emulsion techniques, they found discrete neutron groups at energies corres- 


) Now at Max-Planck-Institut für Physik, Gottingen, Germany. 
) Now at CERN, SC-Division, Geneva 23, Switzerland. | 
(1) P. HuBER and R. WAGNER: Helv. Phys. Acta, 30, 257 (1957). 
) N. S. ANSEROVA et al.: Rec. Trav. Inst. Rad. V. G. Chlopin, 7, 114 (1956). 
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ponding to known levels in *Be. They concluded that the two neutrons of 
the (n, 2n) reactions are emitted in succession. 

ANSEROVA et al. (2) investigated the respective directions of emission of 
the neutrons in the same reaction by the interaction of 14 MeV neutrons. 
By using Be-loaded photographic emulsions they determined the total linear 
impulsion of the two neutrons emitted, by measuring the energies and the 
relative angle of emission of the «-particle fork from the disintegration of *Be. 
The impulsion will be optimal for neutrons emitted in the same direction and 
with the same energy. This process is suggested by observation of pick-up 
reactions (n,d) and (p,d) and by the possible existence of a stable « di- 
neutron ». After examination of 84 forks they concluded that the two neutrons 
are not emitted in the same direction and with the same energy, but in dif- 
ferent directions and with different energies. 

A third possible mechanism of the (n, 2n) reaction would consist in simul- 
taneous emission of the two neutrons by a compound nucleus, è.e. a three-body 
breakup. This process would be favored by simple nuclear structures such as 
magic neutron cores with one or two extra neutrons. 

The relative importance of all three mechanisms considered can be clas- 
sified by the relative probabilities of the existence of two intermediate quantum 
states, in the compound nucleus and in the residual nucleus (neutrons emitted 
in succession), of one intermediate quantum state, in the compound nucleus 
(simultaneous emission by the compound nucleus) and of no intermediate 
quantum state (pick-up reaction). This classification is related to the inter- 
pretation of inelastic scattering 
data from the interaction of 
14 MeV neutrons in the aspects 
of the statistical model. GRAVES 
and ROSEN (°), STELSON and 
GOODMAN (*), and WHITEMORE 
and DENNIS (5), by using pho- 
tographic emulsion techniques, 
measured the distribution in 
energy of the inelastically scat- 
tered neutrons from the inter- 
action of 14 MeV neutrons with SO ito To 300 
various elements. They assigned 


Fig. 1. — Nuclear temperature versus mass 


a parameter analogous to a tem- number of the target nucleus. Data taken 
perature, 7, at an excitation from Ref. (3). 


(3) E. Graves and L. Rosen: Phys. Rev., 89, 343 (1953). 
(4) P. H. Stetson and C. GooDMAN: Phys. Rev., 82, 69 (1950). 
(5) B. G. WHITEMORE and G. E. Dennis: Phys. Rev., 84, 296 (1951). 
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energy E, = (10--13) MeV of the residual nucleus, which is defined as 
(1) T=—[dno(B— £,)/dE,]. 


Fig. 1 shows a plot of their results versus A, the mass number of the residual 
nucleus; the temperature seems to be essentially independent of A. Whereas 
in the aspects of a statistical theory, based on the hypothesis of immediate 
formation of a compound nucleus and equal reduced widths for all channels, 
T would be proportional to A7*. 

The temperatures quoted are deduced from spectra composed of inelastically 
scattered neutrons and neutrons emitted in (n, 2n) reactions, if such reactions 
are energetically possible, the effect of the second contribution being not re- 
solved in these experiments. Assuming successive emission of neutrons in 
(n, 2n) reactions, the correction of their effect on the spectra would increase 
by (20--30)% for the heavier nuclei the temperature due to emission of the 
first neutron, increasing still more the deviations from the prediction of the 
statistical model. GRAVES and ROSEN (3) concluded, that their results would 
not support the hypothesis of simple sharing of the energy introduced into 
the nucleus by all nucleons. 

This situation can be interpreted as well by a considerable contribution 
of simultaneously emitted neutrons, which would be favored by a partial 
sharing of the excitation energy. The maxima of their energy distributions 
would shift to higher energies with increasing A, due to the decreasing binding 
energy of neutrons. In this way the decrease of the temperature with in- 
creasing A would be compensated. 

To investigate the principal mechanism of the (n, 2n) reaction from the 
interaction of 14 MeV neutrons and to decide the above mentioned alternative 
we used time-of-flight techniques to measure energy spectra of neutrons emitted 
by various elements. 


2. — Experimental method. 


21. Time-of-flight technique. — Two measurements on the neutrons must 
be made to establish the zero as well as the arrival time. The time zero has 
usually been established by using a pulsed neutron source (*). This method 
is not recommended for high energy neutrons, because of the difficult back- 
ground problems. Neutron pulses are delivered in 47, giving rise to time- 
correlated background coincidences due to inelastic scattering in the surround- 
ing objects. 


(°) W. WEBER, C. W. JOHNSTONE and L. CRANBERG: Rev. Sci. Inst., 26, 526 (1956). 
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For the experiment described here, a time-of-flight technique in conjunction 
with the (d, t) reaction was chosen (7). A 100 keV deuteron beam bombarded 
a tritium-titan target. The reaction T(d, n)‘He yields a 14 MeV neutron as- 
sociated with a 3 MeV recoil «particle emitted at 180° with respect to the 
direction of the neutron. A plastic scintillator 0.1 mm thick and 5mm in 
diameter within the target chamber of the accelerator detected recoil x-par- 
ticles. The photomultiplier viewing the scintillator is protected against light 
coming from the target by a 0.5 um thick aluminium foil. The «-signals thus 
supplied a time zero for associated 14.1 MeV neutrons in a conical beam, the 
width of which is defined by the solid angle of the a-particle detector. By 
establishing the time zero, a beam of monoenergetic neutrons having a known 
direction is obtained without the use of absorbers and collimators. The neutron 
beam did not meet, for a path of 10 meters, any material except the target 
chamber walls and the studied scatterer. Within this cone, at 20 cm of the 
neutron source, a scatterer, 5 em long and 3 cm in diameter, was placed. It 
was faced by a neutron counter, containing a plastic scintillator, 5 em thick 
and 5 em in diameter, placed outside the time-correlated neutron beam at 
a distance of 50 cm from the scatterer, this distance representing the flight 


n(14MeV) 
RCA 5819 


200MHz FORE 200 MHz 

eee fast |-_—----- SS PETE 
< coincidence 
2% h-p 2*h-p 


timesorter 


= eee ae 


inci hannel 
slow coincidence 100 ¢ 


Fig. 2. — Block diagram of the electronic circuits. 


(7) G. K. O'NEILL: Phys. Rev., 95, 1235 (1954). 
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distance of neutrons emitted by the scatterer. The neutron signal thus supplied 
the arrival time. 

Fig. 2 shows a block diagram of the electronic circuits. RCA 5819 photo- 
multipliers with convex cathode having relatively small transit time fluc- 
tuations (8) were used. They were operated at 1250 V, with equal dynode- 
to-dynode voltages and doubled cathode-to-first-dynode voltage and sub- 
sequent distributed amplification. | 

The electronic circuits consist of a fast-slow coincidence circuit connected 
with a time-to-pulse-height converter. A complete description has been given 

elsewhere (?). The discrimi- 
500r 50 nators in the slow lines are 
adjusted in such a way that 
the slow coincidence output 
signal gates open the am- 
plitude analyser only for 
signals analysed in the fast 
circuit which satisfied the 
required cut-off pulse height. 
The discriminator in the 
slow neutron line served as 
well to establish a well-de- 
fined minimum pulse height 
accepted by the electronics, 
a quantity needed in the 
calculation of the energy 
sensitivity of the neutron 
detector. 

The timesorter has been 
a ania calibrated, using ?2Na an- 


Slope 4.18:10 s/v 
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100} 


Pulse height 
ee 


1 LÌ 


100 75 50 25 0 


nihilation radiation or (n-«) 
Fig. 3. — Calibration of the timesorter using 22Na coincidences, by inserting 
annihilation radiation. The response is linear from known delays in one of the 

8 to 43 ns. One of the resolution .......... fast lines. The response of 

the time-sorter for a given 

delay was defined by the pulse height at the maximum of the corresponding 
resolution curve. Fig. 3 shows a calibration curve as well as one of the 
resolution curves. The response as a function of the delay was found to be 
linear in the region from 8 to 43 ns, representing flight times over 50 cm 
of neutrons of (10--0.8) MeV. In this region the channel width of the am- 


(5) R. V. Suit: Westhinghouse Res. Lab., 60-94511-6-P6 (1956). 
(*) E. Remy and K. WixntER: Journ. Phys. et Rad., 18, 112A (1957). 
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plitude analyser was 4-10-! s/V. Absolute calibration of the time scale with 
respect to the reaction to be studied was obtained by the known flight time 
of y-rays accompanying (n, 2n) events as well as inelastic scattering in the 
studied scatterer. 


2°2. Background. — The number of real time-of-flight events in our, ex- 
periments was about 


(2) 6,= 810-5 Ne 


where N, means the x counting rate per second. For the experiments described 
here the background was due to random coincidences and to time-correlated 
events due to inelastic scattering in the target chamber walls. The number 
of random coincidences is given by 


(3) CON ANA TEN 1 


where N, means the counting rate of the neutron detector and f is a factor 
containing solid angles, sensitivity and possible attenuation of direct neutrons 
by a shield. The ratio of time-of-flight events to direct neutrons on the neutron 
detector being of the order of 10-?, N, can be taken as the counting rate of 
direct neutrons. t, the effective resolution of the circuit, has been measured 
to be 1.01-10-7s, by counting singles and random coincidences of two independent 
y-ray sources. In measuring it by (n-x)-singles and coincidences, without 
the scattering, we found t = 2.00-10-7s, indicating time-correlated coinciden- 
ces due to inelastic scattering events in the target chamber walls. By inserting 
a lead shield between neutron source and neutron detector, the direct 
neutron flux was attenuated by a factor of 10, and 7 was found to be 1.02 -107° s, 
indicating the elimination of time-correlated coincidences. The ratio of real 
to random coincidences was found to be 


Cr 


LA 
GE 


(4) 


Admitting a ratio of 4, an « counting rate of N= 6.7-10° e/s could be 
used yielding 0.5 real events per second. 

The background spectrum was found to be uniform in the linear region 
of the timesorter (random coincidences) and showed a peak in the upper 
channels due to their increasing width in time. In order to establish the sub- 
traction of this peak from a measured time-of-flight spectrum with the same 
statistic as that of the uniform part, we measured its form with very good 
statistics by using two independent sources of y-rays. This spectrum has been 
normalized to the total background rate found in the experiments for a given 
neutron source strength and a given flux. 
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2°3. Neutron detector sensitivity. = The neutron detector energy sensitivity 
is defined by 


, xf BERE 
(5) s-0-1,( = ) 


E, 


Y, meaning the interaction probability of neutrons of energy #, in the 
volume of the scintillator and H™" the minimum recoil energy accepted by 
the slow neutron line. Y; 

5 has been calculated up to 
neutron energies of 4.5 MeV, 

admitting only elastic neu- 

tron-proton interactions for 

detection, VC recoils not 

being accepted by the slow 

neutron line. £" has been 

determined by comparing the 

respective pulse heights of 

proton recoils of 14.1 MeV 

and E?" with the known 

response curve of plastic 

cimev)|  Scintillators. E°" has been 

10 75 20 7.5 3.0 35 zo (500 + 50) keV in all exper- 

Fig. 4. — Neutron detector sensitivity calculated for  iments. The calculated sen- 
a minimum recoil energy accepted by the slow neu-  Sitivity passed through a 
tron line E?" = (500 + 50) keV. maximum at 2 MeV neutron 

energy, its total variation 

between 1 MeV and 4.5 MeV is 25%; the relative sensitivity is shown in Fig. 4. 


24. Time and energy resolution. — The total finite resolution of the time-of- 
flight measurement is determined by two groups of uncertainties; one which 
gives rise to constant relative errors AE/E = 2At/t, another which gives rise 
to energy dependent relative errors. In the first group there are: uncertainties 
due to: 


1) the finite dimension of the scatterer in the direction of the flight path, 


2) the finite dimension of the neutron detector; 
in the second group there are all time uncertainties: 


3) the finite resolution of the timesorter, mostly due to transit time 
fluctuations in the photomultipliers (~ 3 ns), 
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4) the transit time of 14 MeV neutrons through the scatterer, i.e. the 
uncertainty of the moment of the reaction (1 ns), 


5) the uncertainty of the time zero, due to the target thickness and the 
x path geometry (0.5 ns). 


Time-of-flight data have been taken at 90° with respect to incident neutrons, 
so uncertainties due to 1) and 4) can be considered to be approximately inde- 
pendent. Admitting a Gauss function for the different uncertainties, the time 
resolution for y-rays has been computed to 3.2 ns in good agreement with 
the width at half height of observed y lines. The total resolution for neutrons 
over 50cm flight path varies from 10 to 20% for neutrons of (1—5) MeV, 
i.e. neutron groups spaced by 200 keV in the region of 1 MeV can be sepa- 
rated. 


3. — Results and discussion. 


Time-of-flight spectra have been measured for 25-10° «-particles, or 2:10 
neutrons in 47, the total statistics varying from 2500 to 4000. Spectra ob- 
tained with ‘Be, Mn, Fe, 271 and “Ta are shown in Fig. 5-9, the time 
scale is linear, the corresponding neutron energies are indicated. Background 
has already been subtracted, but the spectra have not been corrected for the 


‘Be (n, 2n) 


20 25 30 4.0 50 


Fig. 5. — Neutron time-of-flight spectrum obtained with ‘Be: 
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energy sensitivity of the neutron detector. The peak marked y at 1.7 ns 
of the time scale is due to y-rays from the decay of the residual nucleus. The 
differential cross-section for elastic neutron scattering is small at 90° for heavy 


200 


100 


15 


i 
th (mus) 30 


Fig. 6. — Neutron time-of-flight spectrum obtained with 55Mn. 


nuclei, aS can be seen in the spectra. The spectra after correction for the 
energy sensitivity of the neutron detector are indicated in Fig. 10. They 
show neutron groups at low energies which are attributed to the subsequent 
emission of the second neutron in (n, 2n) reactions. 

We will discuss here the case of ?*I with the help of Fig. 11, which shows 
a diagram of reactions in competition with the (n, 2n) reaction. The com- 
pound nucleus can decay by emission of a neutron or a proton, the latter branch 
leading to the (n, pn) reaction. Table I gives the threshold energies for (n, 2n) 
and (n, pn) reactions, values of the branching ratio, G,/G,, are given in co- 
lumn m, calculated for proton and neutron energies which permit the sub- 
sequent emission of a neutron of 1.3 MeV in contribution to the above men- 
tioned low energy group. In the case of 7I this ratio turns out to be 4-10-?, 
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TABLE I. 
Element Q(n, 2n) Q(n, pn) GG (14 MeV —Q,) 
I II III IV 
55Mn — 11.0 — 7.2 | 1.4-10-4 1.5 | 
54H'6 —1 4.1 — 10.0 = 0 
Fe a2. RSA = | 1.0 
AL — 8.2 — 7.2 49270079 2.9 
181Ma, — 7.7 — 6.2 Aiko INQ" 3.2 


favoring the branch leading to the (n, 2n) reaction. We conclude that the 
low energy neutron group is not due to neutrons emitted in (n, pn) reactions. 


{et 
660° 685 
200 ne 
100 
Eg (n.2n) 
10 7315 20 25 30 40 5.0, 80 1460 En (MeV) 
40 Èp(mys) 30 0 10 vœu 


Fig. 7. — Neutron time-of-flight spectrun obtained with Fe. 


Its interpretation as inelastic scattering of the first neutron would lead 
to a particular variation of the level density with excitation energy. Taking 
f(8)/S, -E proportional to the level density, the latter would turn out to change 
suddenly its slope, a behavior which is not supported by other investigations (7°). 


(10) P. C. Guertor: Phys. Rev., 98, 425 (1954). 
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Furthermore, the shape of the spectrum is in qualitative agreement with 
calculations for successive emission on the basis of the statistical model, achieved 
by WEISSKOPF (1). 


200 


E,(n,2n) 


| 


100 


——— 
ee 
Fmps 30 20 10 y 0 
Fig. 8. — Neutron time-of-flight spectrum obtained with 1271. 


Simultaneous emission of the two neutrons of the (n, 2n) reaction has to 
be considered as a three-body breakup which would imply a continuous distri- 
bution of energy from zero to the maximum energy available. The maximum 
intensity would occur at 3(14 MeV —Q,), where @, means the binding energy 
of a neutron in the target nucleus. Values are listed in column Iv of Table I. 

In the case of !*"I no indication for simultaneous emission has been found 
(Fig. 8 and 10). The listed energy value for maximum intensity occurs on the 
slope of the continuum due to the emission of the first neutron. Simultaneous 
emission by the compound nucleus will be less probable for heavier nuclei 


(11) V. F. WEISsKOPF et al.: NYO, 636, 153 (1951). 
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(compare !8!Ta in Fig. 9 and 10), the reduced width for this process being 
smaller in the presence of a still greater number of nucleons in the compound 
nucleus. 


200 


100 


£,(n,2n) 


Fig. 9. — Neutron time-of-flight spectrum obtained with 1!81Ta. 


The neutron spectra observed for **Mn and Fe finally provide a quite sen- 
sitive test for simultaneous emission. The neutron configuration (magic +2) 
should favor the pick-up as well as simultaneous emission by the compound 
nucleus. Again, there is, in the limit of statistics, no indication for a three- 
body breakup. We conclude that this mechanism is very improbable for me- 
dium weight nuclei (and for incident neutrons of 14 MeV) in comparison with 
successive emission, even if the probability of a simultaneous distribution of 
energy on two neutrons in the presence of other nucleons is favored as in the 
cases of **Mn and *Fe. 

It is concluded that the principal mechanism of the (n, 2n) reaction on 
medium weight and heavy nuclei induced by 14 MeV neutrons is successive 
emission. The consequences for the interpretation of inelastic scattering data 
will be discussed in Sect. 4. 
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The pick-up mechanism in (n, 2n) reactions on *Be is favored by the simple 
structure of the target nucleus. HAEFNER (1?) has shown that *Be can be con- 
sidered as a system of a weakly bound neutron in the average potential of 


30 20 


Fig. 10. — Neutron time-of-flight spectra corrected for the energy sensitivity of the 


neutron detector. 


two preformed «-particles. This simple structure would favor as well a three- 


body breakup of the compound nucleus. 


The energy spectrum observed for 


‘Be (Fig. 5) is simple because inelastic scattering is impossible or very impro- 


target 


Fig. 11.— Competition diagram for the interaction 


of 14 MeV neutrons with 


bable since the ground state 
energy of ‘Be is lower than the 
first excited state in *Be. So 
only elastic scattering and (n, 
2n) reactions have to be con- 
sidered. This simplification is 
supported by the small number 
of y-rays observed, which are 
emitted in competition with 
the second neutron. We have 
listed neutron energies expected 
for successive emission to the 
levels of °Be and *Be in Table IT, 


(12) P. R. HAEFNER: Rev. Mod. Phys., 28, 228 (1951). 
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TABLE II. 

n,(*Be) ny (Be) 
14 MeV ces | = = sa | == 
12.2 = ee = = =i 
11.5 08 ft dr = | = 
| 10.9 1.4 | = = = = 
| 9.2 3.1 | — — _ | — 
10) 5.1 | 2.2 0.9 | = = 
| 6.1 | 6.2 | 3.3 | 2.0 0.8 - 
2.7 | 9.6 | 6.7 | 5.4 1 | 1.95 


calculated for all levels reported by AJZENBERG (!*). Neutron groups at 1.3 MeV 
and 2 MeV are resolved, agreeing with calculated transition energies to levels 
of ‘Be. Neutron energies between 5 and 12 MeV are not resolved, containing 
either a great number of successive transitions (see Table II) or the maximum 
intensity for simultaneous emission or a mixture of both. In the spectrum 
observed there is no indication of the transition to the 11.3 MeV level in °Be 
which would yield 2.7 MeV neutrons. We cannot decide whether this level 
does not exist or whether this transition does not occur because the time 
between the emissions of the first and the second neutron at this excitation 
energy is very short, e.g. of the order of the transit time of nucleons over 
nuclear dimensions, so that there is no defined excited state between the two 
emissions. 

A contribution of the pick-up mechanism could be identified by its angular 
correlation. Experiments on this question are in progress. 


4. — Interpretation of inelastic scattering data. 


The principal mechanism of the (n, 2n) reaction in medium weight and 
heavy nuclei is successive emission of the two neutrons. The compound 
nucleus and the residual nucleus have lifetimes which are long compared with 
the transit time of nucleons over nuclear dimensions: ¢.9. there are well de- 
fined intermediate quantum states between the emissions of the first and the 
second neutron. A partial distribution of the excitation energy, on a small 
fraction of the nucleons, would favor simultaneous emission. 

The variation of the nuclear temperature versus A (Fig. 1) has consequently 
to be interpreted in the concepts of the compound nucleus. The correction 
for the effect of the subsequent emission of a second neutron, if energetically 


(13) F. AJZENBERG and T. LAURITSEN: Rev. Mod. Phys., 27, 77 (1955). 
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possible, increases the temperature due to the emission of the first neutron 
by (20--30)% for the heavier nuclei, increasing still more the deviations from 
the predicted A-* behavior. The only hypothesis involved in this prediction 
which could have been violated is that of equal reduced widths for all levels. 
There is some evidence for a corrective transition factor, which favors certain 
transitions and forbids others. The importance of this factor would increase 
with increasing mass number. This suggestion is supported by the form of 
spectra observed in inelastic scattering of protons (1). LANE and WILKIN- 
SON (#4) have interpreted this «selection rule» by the concept of fractional 
parentage. The energy spectrum of inelastically scattered neutrons would 
consequently not provide a direct measure of the level density of the residual 
nucleus. 


We would like to express our gratitude to Professor FRANCIS PERRIN for 
his hospitality and for his continued encouragement. 
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(14) LANE and WILKINSON: Phys. Rev., 97, 1199 (1955). 


RIASSUNTO (*} 


Si sono misurati gli spettri dei neutroni emessi da ®Be, 5°Mn, Fe, 1271 e 18!Ta bom- 
bardati con neutroni di 14 MeV servendosi di tecniche basate sul tempo di volo. Si 
separano gli spettri continui dei primi dai gruppi di bassa energia dovuti ai neutroni 
incidenti. Nei nuclei di peso medio e pesanti è predominante l'emissione successiva. 
Non si può decidere sull’esistenza di un possibile contributo da parte di un’emissione 


simultanea dai nuclei leggeri. Si usano i risultati per interpretare i dati dello scattering 


anelastico derivante dall’interazione coi neutroni di 14 MeV. 


(*) Traduzione a cura della Redazione. 
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A New Possible Theory of the ;-Meson. 


I. ULEHLA (*) 


Joint Institute of Nuclear Research - Moscow 


(ricevuto il 30 Aprile 1958) 


Summary. — The paper presents a theory of the y-meson as a particle 
with spin + described by anomalous equations differing from the Dirac 
equations. Calculation is carried out of the intrinsic magnetic moment 
of the u-meson, the cross-section for Coulomb and Compton scattering, 
the self-energy and the electrodynamic correction to the magnetic mo- 
ment of the u-meson. An explanation is given of the similarity of the 
u-meson and electron at low energies and their dissimilarity at very high 
energies in electromagnetic interaction. The difference in the masses of 
the electron and u-meson is determined by means of the same electro- 
magnetic cut-off for both particles. 


Introduction. 
It is known from experimental material that for electromagnetic inter- 
action in the non-relativistic region the y-meson negligibly differs from the 


electron. For example, its magnetic moment, measured by LEDERMAN (!), is 


m, = 1.0024 +0.0008 , 


while the magnetic moment of the electron has the value 


m, = 1.001167 + 0.000005 . 


(*) Now: Institute of Nuclear Physics, Academy of Sciences, Prague. 
(1) L. M. LEDERMAN: Proc. of the Internat. Conf. on Mesons and Recently Discovered 
Particles, Padua-Venice (1957). 
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Measurements of the cross-sections for y-meson scattering in the relativistic 
region (*) and other such measurements point to the fact that at high energies 
the u-meson does not behave like an electron. If these measurements are correct, 
the difference in the behaviour of the u-meson and the electron in the rela- 
tivistic region will have to be clarified. 

Up to now the relatively large difference in the masses of the p-meson 
and the electron have not been satisfactorily explained. As long as it was 
not definitely known that the meson has spin 4 it was thought that this dif- 
ference could be explained by means of the u-meson as a particle with spin $ (*). 

As will be seen, the effects mentioned can be explained by describing the 
u-meson by anomalous equations which were found in paper (*). The existence 
of such equations for particles with spin & seems at first debatable if we 
remember the proofs of WILD (5) and GEL'FAND and YAGLOM (6). These proofs 
are not however general, as is shown in paper (*). 

Of all permissible anomalous equations we chose the simplest because it 
is the nearest to the Dirac equations from the structural point of view. These 
anomalous equations are characterised by the fact that they contain one 
free parameter called the parameter of anomality g. They are complicated 
and the calculations carried out were therefore not simple. In order to avoid 
mathematical difficulties we confined ourselves to the lowest approximations 
of the perturbation theory. 

When working with anomalous equations the general method elaborated 
in quantum electrodynamics can be used without change except that the Dirac 
matrices y, must be replaced by new matrices and the changes incurred by 
this substitution taken into consideration. 

We shall first deal with the equations for a free y-meson and will then 
show that the anomalous equations are compatible with «switching-on » of 
the electromagnetic field. We shall investigate the electromagnetic interaction 
of y-mesons in detail and in conclusion will briefly mention other types of 
interaction. 


1. — Free u-meson. 
The motion of a free u-meson is described by the equation 


(1) (BY 0, —ix)p = 0, 


J. Luoyp, A. W. WOLDENDALE: Proc. Phys. Soc., A 68, 1045 (1955). 
J. MATHEWS: Phys. Rev., 102, 270 (1956). 

I. ULEHLA: Zu. Éksper. Teor. Fiz., 33, 472 (1957). 

E. Witp: Proc. Roy. Soc., A 191, 253 (1947). 
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where the wave function g has twenty components. The matrices f, are 
twenty-row irreducible matrices which satisfy the relations 


(2) Dai Ca vay, di 5200) =0 ’ 


where SE denotes summation over all permutations of the indices. The re- 
lations (2) are not the only ones satisfied by these matrices. There also exist 
certain relations of the fourth order such as 


BoP1 Ps Po 2 Bo Ps Pr Bo = 0 , 


which are an analogue of the equations 


Bobi Ps =) o) 


valid in the Duffin-Kemmer algebra for five-row representation, and which 
point to the fact that we are dealing with a special representation of some 
more general algebra. Relations of fifth and higher order in f,, not following 
from (2), also hold for the matrices poe We shall deal with the algebra of 
the matrices f, in the appendix. The matrices 5, contain one free parameter 
— the parameter of anomality g — whose value lies in the intervals (— co, 0) 
and (0,4). For g=0 wave equation (1) becomes the Dirac equation. 

If we apply the operator Ba, three times to Eq. (1) we obtain the Schro- 
dinger-Gordon equation 


(O+#)p9 = 0, 


from which it follows that the particle has a single mass x. 
Since there exists a matrix A= A* satisfying the condition 


AB, == 024) 5} 
we can derive wave equation (1) from the Lagrangian 
(3) L= i(gp" 0,9 — ixpp) , 


where p= g*A. From this function we can derive the general momentum- 
energy tensor which is symmetrical and Hermitian 


na a E 
(4) Ty À PBr0uP + 92 0°{9 [Polat By Ih ait L,61P} ; 


where 1,, are generators of infinitesimally small rotations. We can also derive 


44 - Il Nuovo Cimento. 
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from (3) the current density 
(5) j, = — PB ,@- 


We shall now show that wave equations (1) contain 16 subsidiary con- 
ditions by which the superfluous components are eliminated and that equa- 
tions (1) describe particles with spin }. The subsidiary conditions can be 
derived the quickest in the spin-tensor notation of wave equations (1), in 
which these equations have the form 


(6a) (204% y Ory Hy) 42 PB", 
(65) 28, 0,BY + ky" 0,4? = ing, 
(6e) 38, 0, B2 a k,y" OS = tay ? 
where 
DRE D k 
tb Sg yey Gao 


= Sor k, + k: = À, kg; k= 0: 


1=l = el for, ga sh = 0, and A= ye, or 0) 
The twenty components of the wave function y are formed by two Dirac spinors 
XP and y and 12 independent components of the spin-vector B“, for which 
the following condition is valid 


(7) ye. ai 0 A 


Condition (7) also follows directly from Eq. (6a) if it is multiplied by the 
operator y,,. 
From Eqs. (60) and (6c) we easily obtain the equation 


(3) ye ds (ri ge To”) — Tx & yD + 3° n°) 


Vi 
by means of which (64) can be rearranged to 


(9a) 40u(r, 40 + yy) — ing (; 0 + È a — inBe, 
“il 2 


In Eqs. (9a) we have 12 subsidiary conditions (for u —1, 2,3) not containing 
time derivatives. 
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The four remaining subsidiary conditions are obtained by applying the 
operator 0, to Eq. (6a) and the operator y’d, to Eq. (8). The second deriv- 
atives are eliminated by combining these two equations. By means of Eqs. (6b, c) 
we eliminate the expressions 70,47 and y*o 4 from the resultant equation 
and finally obtain 


Alla 3 181 | 3 7282 jp 
Bin (Tx Ted se 2 E Paty Die 


1 “a 


Since the expression in the brackets at 0,B" is zero, the last four subsidiary 
conditions read 


(9b) Hye aes — OF 


12 | no 


Tf we consider the spinor 4° independent, we obtain for it from Eq. (8) 
by means of (9b) 


(10) (vasti “0-0, 


which is the Dirac equation. If we know its solution — 4.6. the spmor UE 
we determine the spinor © from conditions (9b) and the spinvector B, from 
conditions (9a) and from Eq. (7). 

That we obtained the Dirac equation for the spinor 7 was only to be 
expected. This result shows that we are actually dealing with a particle having 
spin 3. It cannot be inferred from this however that the y-meson described 
by anomalous equations is identical with an electron, for the energy-momentum 
tensor (4) and current density (5) cannot be reduced by means of the subsidiary 
conditions to the energy-momentum tensor or current density for an electron. 


9. — Electromagnetic interaction of the u—meson. 


If the problem can be formulated so that the wave equations for the 
u-meson can be derived from the Lagrangian, then the « switching-on » of the 
electromagnetic field leads to no difficulties. The interaction term can be de- 
fined in the «classical » way as 


(11) Dés Cr ila 7 


i.e. the electromagnetic interaction of the p-meson can be assumed normal. 
Although the algebra of matrices B permits the construction of a number 
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of different and independent vectors such as 


By PBB BiB > ME 


we choose for the interaction term only that vector the divergence of which 
is zero — the vector of the current density (5). 

For the case of electromagnetic interaction the number of subsidiary con- 
ditions does not change but their form differs from (9). After « Switching-on » 
the electromagnetic field we write D,—¢,—ieA, instead of 0, in Eqs. (1) 
and (6). Proceeding in the same way as when deriving the subsidiary conditions 
in part I, we derive the following equations instead of Eqs. (9a) 


n GP, Ve : 
(124) ADM rg + 12710) — ig [TE y + xo) = inte, 


VI 2 
and instead of (95) 


È ‘9 : il 
(120) x (7 hoy ae = kx) = deg I’ (3) Polti 4 de 42) 7 
vil 2 2 


where 1(4)= }(y“y" — yy"). When deriving Eqs. (125) we used the relations 


[D,D,) = — ieF yy DD, = — iel'(4)F,,4 DD 


uv I uv 


Equations (12) again contain 16 subsidiary conditions and always have 
a solution; the introduction of electromagnetic interaction does not therefore 
lead to any incompatibility. 

The equations for the y-meson have a certain relation to the Pauli equations 
for a particle with anomalous electromagnetic interaction. If we introduce 
a new spinor y defined by 


y = Te + vor ; 
then by means of conditions (125) we can easily derive for it 


e 


av IRE .% 
(13) pon + gl' BEE =: 


A 


which is known as the Pauli equation. At first glance it would seem that the 
remaining components of the wave function » can be forgotten and only the 
Pauli equation used. However the transition from equations for the y-meson 
to Eq. (13) in the sense implied does not mean only a transition from one 
field to another, as we have already mentioned, but also a transition from a 
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particle with normal interaction to a particle with anomalous interaction. 
A particle described by Eq. (1) (with DEAD) cannot therefore be considered 
equivalent to a particle described by Eq. (13). As will be seen below, the 
results of the calculations differ in both cases. 


91. Intrinsic magnetic moment of the p-meson. — In a non-relativistic approx- 


imation the intrinsic magnetic moment of an anomalous particle — of the 
u-meson — is easily determined. Its value, measured in Bohr magnetons, is 
(14) m =—1--29, 


and it is seen to be anomalous. If g>0 then 


L< m, <1 : 


while for g< 90 
lime Oe 


The magnetic moment of the u-meson was measured very accurately by 
LEDERMAN (!), who obtained the value 


m, = 1.0024 40.0008 . 


This can serve as a basis for estimating the magnitude of the parameter of 
anomality g. The parameter of anomality cannot be too large since it can 
hardly be expected that the electrodynamie correction to the value of the 
magnetic moment of a y-meson would be much bigger than corrections in 
the electron theory. Thus 9 will have a value of the order 


(15) \g| = 10°. 


9°2. Coulomb scattering. — When calculating Coulomb scattering of a u-meson 
on a fixed centre we shall confine ourselves to the first approximation of the 
perturbation theory. The vertex of Feynman’s diagram is represented by the 
operator ep, and the amplitude of the transition is therefore given simply 
by the expression 


UsPoUr - 


The cross-section is determined, as is known, from the expression 


eZ\? n2E? Meo 
(16) dg = | o | pid — cos OY 2 Sp Pz BoPi Bodo , 
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where H=v/p?+w and pi, are projection operators which, with respect to 
relations (2), have the form 


a: a | sali PSE 
x 


The scalar product flp, is denoted by (BP). 
If we calculate the trace in (16) we obtain the cross-section for the scat- 
tering of the u-meson 


x 


(17) do = (=) fel ee AS ie L 1+ cos 0 + 4g (1 +17) (1 — cos o). 
p? (1— cos d)? | p? di 

If we put g=0 in this expression we obtain the cross-section for scattering 
of the electron. The do obtained differs from the cross-section for Coulomb 
scattering of a particle described by the Pauli equation (13) as is seen in 
paper (?). 

In the non-relativistic approximation do is the same for the electron and 
u-meson since g is small. In the extreme relativistic case, however, do for a 
u-meson is independent of its energ 


sin Ÿ dd 
1— cos 3’ 


pi 


— 
ft 
(oe) 


e2Z\2 
do | 
\ x 


and its angular dependence differs from the angular dependence for an electron. 
Both deviations obtained agree with those discovered by J. LLOYD and 
A. W. WOLFENDALE (?) when investigating the scattering of high energy 
u-mesons in cosmic radiation. The region of energies >10° MeV, in which 
the deviations were observed, is also in agreement with expression (17) if 
we put |g|< 10-°. 


2°3. Compton scattering. — In calculating the scattering of a y quantum on 
u-mesons we again confine ourselves to the first approximation for the sake 
of greater mathematical simplicity. For Compton scattering the transition 
amplitude will have the known form 


M = us(Bex){(P, Pr + ky) — x {Be} + (Be, )L(B, Pi — ke) — x (Pes), , 


where p, is the momentum of the u-meson, k, the momentum of the + quantum 
and e, the polarization. 
The inverse operator must be known in order to calculate the traces. This 


(7) H. C. CorBEN and J. SCHWINGER: Phys. Rev., 58, 953 (1940). 
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is easily found if we start from Eqs. (2) 


1(BD=—x1 > = = ae tas + (BQ) 
qq) — % 


We determined the exact values of the cross-sections do, and do, for light 

I Al È 
polarized after scattering parallel and perpendicular to the plane of polarization 
of the incident light: 


1 ek; i k, — ka 
do, = 4B [aca 4 cost 941-891 hn + 16041) haw, 
(19) Vj 12 
Lek | Ge cope lena Jes AA 
| do, — 12 ta — 2g)(1— 69)7 + 169?(1 — 29) - ay + 16g’ ai do , 


where 7 = k;/k:+ ka]k1 — 2. 
In the general case the cross-section is given by the sum 


do = do, + do, . 


For g=0 both cross-sections become the corresponding expressions known 
from the theory of the electron and for g~0 differ from the cross-sections 
for a particle described by the Pauli equation (13), as can be seen from paper (*). 

By integrating expressions (19) we obtain the total cross-section. Here, 
however, we only give their values in the non-relativistic and extreme rela- 
tivistie approximations. For an energy of the incident quantum k, &%* up to 
to the linear term in hk, 


È 82 e! | 2k, 
(20) O non-rel. — ee 2 x ’ 


from which it follows that Orma. for an electron and p-meson are identical. 
For k,>> however 


4 
(21) c3=l4agi SAL 


In this case the cross-section increases linearly with the energy. The origin 
of the linear growth lies very high (> 10° MeV) due to the factor g in ex- 
pression (21). 


(8) S. B. BATDORF and R. Tuomas: Phys. Rev., 59, 621 (1941). 
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2°4. Self-energy of the u-meson. — The self-energy of the p-meson is determined 
by starting from the integral 


e? 
Ari 


i il BALE, p— 8) — 414) Be dE | ue 


Calculation is carried out by the method of Pauli and Willars (°), by which 
we use one auxiliary mass M = Ax only. 
We get 


il 7 PA 
+ (3 — 6g + 5g?)ln À = (1 29g 4 7. 2] uu (*). 
For g=0 we have 


62% 1 
(2) 
Ee. = > ( In A "al À 


a 


which is an expression well known from the theory of the electron. 

It is seen from expression (22) that the self-energy of the p-meson is qua- 
dratically divergent. If we consider the parameter A as the parameter of cut-off 
we can explain the difference between the mass of the u-meson and that of 
the electron. If we choose the value 4 ~ 8:10? for |g| = 10-* and the value of 
the parameter 1 ~ 2.5-10* for |g|=10-* we get the ratio of the masses of 
the u-meson and electron 

VA 


LES SOND) 5 
He 


2°5. Electromagnetic correction to the magnetic moment. — In calculating 
corrections to the value of the magnetic moment we limit ourselves for mathe- 
matical reasons to the correction of the second order. This means that we 
must determine the correction to the vertex part of the Feynman diagram. 
As is known, this is given by the integral 


e? 


CHE i dp — OI *Bul(By 4 + k) — x} B,[(Bq) — 2180 aq, 


(*) A numerical error is contained in the expression for the self-energy in the 
pre-print called Electromagnetic Interaction of Anomalous Particles with Spin 4. 
(9) W. Pauri and F. WILLARS: Rev. Mod. Phys., 21, 434 (1949). 
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where the meaning of the individual quantities is seen from the diagram 


p DE p+K 


The calculation of integral (23) is carried out by the same method as was 
used in calculating the self-energy. The correction (23), in the same way as 
the self-energy, contains logarithmic and quadratic divergence but does not 
contain terms with divergences of higher order. In the non-relativistic approx- 
imation, which we obtain by multiplying expression (23) from the right and 
left by the projection operator 


IRA + Bo); 


all terms with quadratic divergence vanish for y=1, 2,3. In this approx- 
imation, which is the only interesting one, we can therefore limit ourselves 
to terms of the zero-th and first order in g and first order in k. 

Tf we then determine the value of the magnetic moment from the gene- 


ralized interaction term 
pp, + B,)A"p 1 


we obtain, after subtracting expressions with infra-red catastrophe, 


my = 1— 29 + - [1 + g(a + bn À)], 


where a and } are constants ~1; thus the magnetic moment of the u-meson 
will be 


mr 1 29 5 


if we consider the values of À from Sect. 2°4. The result we obtain confirms 
the correctness of the estimate (15) of the magnitude of the parameter of 


anomality g. 
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3. — Conclusion. 


This new interpretation of the properties of the u-meson is based on one 
assumption. We assume the u-meson as a particle described by the simplest 
anomalous equations, which contain one parameter of anomality g. The small 
value of this parameter explains why the p-meson, interacting with the electro- 
magnetic field, does not differ from the electron at normal energies and shows 
that differences between the two particles can only be expected at high energies 
(> 10° MeV). It also permits a suitable explanation of the difference in mass 
of the u-meson and the electron. 

So far we have not dealt with questions of other types of meson inter- 
actions than electromagnetie interaction. This problem is complicated by the 
fact that several linearly independent Lagrangians can be constructed for the 
given type of interaction. For example, for scalar interaction we do not have 
only one scalar operator I, as in the Dirac algebra, but also other linearly 
independent scalars 


(BB), (880, B,BBB - 


The same is the case with vector interaction, tensor interaction ete. In general 
not all these operators are permissible, since some can lead to disruption in 
‘the number of subsidiary conditions and must therefore be eliminated. But 
even then we still have more coupling constants than in the Dirac theory. 

Experimental data on y-mesons are not yet so extensive that they could 
not be explained by a suitable choice of constants. We did not therefore think 
it of any great use to deal with other types of interaction except the decay 
Ta Ve 

If we describe the u-meson by anomalous equations we can regard the 
electron as a limiting case of the u-meson (9g +0). As can be seen from 
expressions (17) and (19) the transition probability for an arbitrary process 
with a y-meson will always be a polynomial in g. 


(24) Wim F Ag + dog? +... + Ang”, 


where the coefficients a; depend on the momentum. If we put g = 0 in this 
expression we get the probability of transition for the electron 


W, = Ao è 


If we now consider the 7 >u+y decay as a direct process and choose the 
Lagrangian for it so that the coefficient a, in expression (24) is zero — this 
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can be done for arbitrary coupling for example by introducing the operator 
(BB)-4 into the interaction term — we get a finite probability W, for the 
m—>u-+v decay and the probability W,= 0 for the decay x —e+v. For the 


processes of the type Up = DIEEVAWO must choose such Lagrangian (with 
S, T, V coupling only), which leads to non-zero coefficient a, in (24). 


In conclusion the author would like to thank N. N. BOGOLJUBOV and 
A. SALAM for encouraging him in this work which is full of mathematical 
difficulties. 


APPENDIX 


We give here the algebraic relations necessary for the calculations in this 
paper. An important factor which greatly facilitates calculation of the traces 
is the possibility of writing the matrices 7 in the form of a direct product 
(A.1) Bey pn oes Bu = yh X oy - (no summation). 
The matrices y, are the well-known Dirac matrices, and the matrices «, form 


a certain algebra A with which we shall deal in detail later. 
If the matrix % is known the matrices «; (t=1, 2, 3) can be derived suc- 


cessively by. means of the unitary transformations 
a, = (1 — w)e,(1 — w) , 
a = (1— 4) (1 — U) ; 
dr ARMATA 


The elements wu, v, w are Hermitian operators satisfying the relations 


ua lu, v= 20, w= 2W, 
(A.24) uvu=U; MU ig OUD Ie ROE, 
vw =wv=0. 


The algebra A, which is formed by the elements x, WU, 0, w 18 determined on 
the one hand by the relations (A.2a), and on the other by the minimum 
condition for % 


(A.2b) ad = Aa , where 14 EME 
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and the relations between the element « and the elements u, v, w: 


ko = Uy = HY = VKH = 0, 


WU 0 
(A.2c) À 
Wow = gw, 


XW + Lp WG = ÀxgWx, + YX - 


From relations (A.2) it follows that the algebra 2{ has 26 linearly inde- 
pendent elements of the basis and two elements of the centre. It thus has 
two non-equivalent representations, one one-row, which is trivial, and one 
five-row: 


toe vas 
ry V2 2 3 Vi 
mal. à Vir, —Vir, | 4=| vi Vi 


3 Vs, ki 


So a | Vis, . ko 


v = diag|2, 0, 0, 0, 0|, w = diag |0, 2, 0, 0, O|, 


where 
2 k 2 k 
"8, =~ ki —*— eg SE: 
oO Do ano 
Kki+k: = À, Ta = E1S1 » Va = E82 » g=4kka, it, SOR [fe So se 
ene 1 10M Aa and 6 — 6== — _L for 1 — | 


Of all the traces of linearly independent elements of the basis of the 
algebra % the only non-zero ones are 


Spee = AG Spiu =Sp? Sp, 
Spo = Spuv = Spuw =1, Spagw=g. 
The generators of the infinitesimal rotations are given by the expressions 
In =yoyrX(¢—), Ia=%y:X(3-%, Ia= Yea X (3 — 0). 


The remaining I,, are determined from the above by means of the known 
relations for these operators. By means of a direct product we can also express 
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the matrix of spatial inversion as 
J Sel; 
and the matrix, which we need for constructing ¢: 


A = vo X diag |1, ily eros 


RIASSUNTO (*) 


Il lavoro presenta una teoria del mesone p come particella di spin 4 descritta da 
equazioni anomale differenti da quelle di Dirac. Si esegue il calcolo del momento ma- 
gnetico intrinseco del mesone y, della sezione d’urto per lo scattering coulombiano € 
Compton, della self-energia e della correzione elettrodinamica del momento magnetico 
del mesone y. Si dà una spiegazione della similarità del mesone y con l’elettrone alle 
basse energie e della loro dissimilarità alle altissime energie nell’interazione elettro- 
magnetica. SÌ determina la differenza fra le masse dell’elettrone e del mesone p. per 
mezzo di un cut-off elettromagnetico uguale per le due particelle. 


una ci 


(*) Traduzione a cura della Redazione. 
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On the Determination of the Space-Group 
by the Statistical Method. 


M. S. FARAG 


National Research Centre - Dokki, Cairo 


(ricevuto il 6 Maggio 1958) 


Summary. — Following A. J. C. WiLson and his collaborators, it is pos- 
sible, by a study of the statistical distribution of X-ray intensities of Bragg 
reflections, to differentiate between centro- and non-centro-symmetric stru- 
ctures. This method was applied to a centred and a non-centred projection 
in the crystal structure of 1.3.5-triphenylbenzene. The results in both zones 
showed big discrepancies which could lead to an incorrect conclusion. 
These discrepancies are discussed in terms of the established structure 
which is of a layer type. The molecules are projected edgewise in the 
centred projection, and the electron density distribution appears in the 
other one, with some approximations, as if centred. 


1. — Introduction. 


One of the important early steps in X-ray crystal structure analysis is the 
determination of the space group symbol. From the systematically absent 
Bragg reflections, it is possible to deduce first the lattice type and then any 
screw axes or glide planes with their appropriate translations. The parts of 
the space group symbol which are then left to be found by methods other 
than X-rays are: rotation axes and mirror planes (1); the apparent presence 
of which in X-ray diffraction photographs is due to Friedel’s law. 


(*) N. F. M. Henry and K. LONSDALE: International Tables for X-ray Crystallo- 
graphy, Vol. 1 (1952). 
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9. — The statistical method. 


Methods usually employed to find the complete space group symbol include 
studies of: crystal morphology, etch figures, optical activity, pyro- and piezo- 
electricity. However, a new method using X-ray data has been recently worked 
out by WILSON (?) and developed further by HOWELLS, Pres and ROGERS (°). 
According to these authors, if systematic absences and reflections having 
sin 0 < A/a, (where a is the shortest axis in the unit cell) are ignored, then 
the ratio <|F|>2/<I> has the value 0.785 for a non-centro-symmetrie and 0.637 
for a centro-symmetric structure, where the brackets < > indicate the average 
value for a certain sin Ü group. Accordingly if each intensity is expressed as 
a fraction Z, (Z=I[KI) of the local average intensity) then the fractions N(Z) 
of the reflections whose intensities are less than or equal to Z will be: 


N(Z)=1— XP [—e] for non-centro-symmetric structures designated by (1), 
and 


1 


N(Z)=erf (32) » centro-symmetrie structures designated by Gis 
where erf is the error function. 
The assumptions underlying the validity of this theory are briefly : 


(i) that there should be a large number of atoms all occupying general 
positions and such that none of them dominate the intensity distribution. 


(ii) that a sufficiently large number of intensities should be used in 
calculating the averages. 


3. — The case of 1.3.5-triphenylbenzene (CH) (CH 5)s « 


The statistical method has been successfully used in many cases (3), but 
the space group of 1.3.5-triphenylbenzene presents an interesting deviation 
from this method. The crystal was first examined by X-rays by LONSDALE 
and ORELKIN (*), who found the unit cell dimensions and space group. The 
latter could be then Pna 2,, or Pnam, but since the crystal showed a positive 
pyro-electric effect, i.e. the crystal is polar, it was then unequivocally settled 
as the non-centrosymmetric one i.e. Pna2,. The crystal structure was fully 


(2) A. J. C. WILSON: Acta Oryst., 2, 318 (1949). 
(3) E. R. HOWELLS, D. ©. Prices and D. ROGERS: Acta Oryst., 3, 210 (1950). 
(4) K. LONSDALE and B. ORELKIN: Proc. Roy. Soc., A 144, 630 (1934). 
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elucidated by the author (5). The unit cell which is orthorhombic has the 
following dimensions: 


a =7.47+0.0249°, b=19.66+0.054°, ¢ =11.19+0.03A°, 


with four molecules per unit cell. The c-projection is the only one with a 
centre of symmetry. The number of independent reflections observed and ex- 
pected theoretically for each zone are given in the following table: 


(hkl) Number of reflections | Number of reflections 
observed expected 

okl 114 155 

hol | 49 61 

hko | 130 | 195 


The probability distribution of X-ray intensities was studied for the (okl) 
and (hko) reflections; the (hol) type were ignored on account of their limited 
number. The intensities used were corrected for the Lorentz and polarization 
factors. Overlapping sin 0 
groups were found useful. 

The experimental results 
for the 2 zones are shown in 
Fig. 1 together with the 
theoretical statistical distri- 
butions for centro- and 
non-centro-symmetrie struc- 
tures. Fig. 1 shows the 
following details: 


1) The (hko) reflect- 
ions: The intensity distri- 
bution for these reflections 
appears as a hyper-centro- 
symmetric one, i.e. as if the 
molecule had a centre of 
Symmetry in a centro-sym- 
metric projection. Such examples have been shown by Lipson and WooLr- 
SON (°) in the case of pyrene and 1:1:6:6 tetra-phenyl-hexapentaene. In the 


80+ N(2 %o) 


(5) M. S. Farag: Acta Cryst., 7, 117, 380 (1954). 
(°) H. Lipson and M. M. WooLrson: Acta Oryst., 5, 680 (1952). 
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present case the deviation may be due to the fact that the separate parts 
(phenyl groups) of the molecule are centro-symmetrie, besides a violation of 
the first condition of the requirements necessary for the application of the 
theory. As has been shown by the author (5), the structure is of a layer 
type where the molecules are seen edgewise in the c- projection. 


2) The (okl) reflections: The intensity distribution for these reflections 
shows that the a-projection is centro-symmetric rather than non-centro-Sym- 
metric; a conclusion which is simply incorrect. The discrepancy between the 
curves (a-) and (1) in Fig. 1 may have the following explanation. 

The a-projection, Fig. 2, has the two dimensional space group symbol em 1, 
and the elements of symmetry are shown in the key diagram of the same figure. 
Although em 1 is not a centro-symmetric space group yet the molecules are 


= A [ANY 
AID TERA 
Caos DERE 
SAG) ET 3 
GC? SO e See 


Ra SNES 


S 


= CSA 
BASE a IS 


arranged in such a manner that, to a first approximation, the benzene rings 
A, B. C, ... ete. are apparently related to A’; B’, C',=. etc. by à centre of sym- 
metry. This reason, of course, may be considered as a violation of the first 
condition for the requirements of the application of the theory, but unfor- 
tunately such an explanation could only be given at a late stage in the analysis 
of the structure. 


45 - Il Nuovo Cimento. 
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4. — Conclusion. 


From this example, and in the opinion of the author, a centre of symmetry 
may safely be detected by the statistical method if hyper-centro-symmetry is 
observed, otherwise physical tests are to be used. Such tests usually are powerful 
enough to place the crystal either in the hemihedral or in the holohedral class. 


The author wishes to record his deepest thanks to Prof. K. LONSDALE, 
F.R.S., D.B.E. under whose supervision the determination of the crystal struc- 
ture was carried out. 


RIASSUNTO (*) 


Secondo A. J. C. WILSON e coll. è possibile, per mezzo di uno studio della distri- 
buzione statistica delle intensità dei raggi X delle riflessioni di Bragg, distinguere le 
strutture centro-simmetriche dalle non-centro-simmetriche. Questo metodo è stato ap- 
plicato a una proiezione centrata e a una non centrata della struttura cristallina del- 
1°1.3.5-trifenilbenzene. I risultati mostrarono nelle due zone gravi discrepanze che po- 
trebbero portare a conclusioni errate. Si discutono tali discrepanze in termini della 
struttura accertata che è del tipo a strati. Le molecole si proiettano di taglio nella proie- 
zione centrata e la distribuzione della densità elettronica appare nell’altra proiezione, 
con qualche approssimazione, come centrata. 


(*) Traduzione a cura della Redazione, 
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Summary. — Conversion coefficient of 800 keV y-rays from the decay 
of Polonium-210 has been measured using scintillation method to be 
0.039 + 0.01. The value of conversion coefficient has also been obtained 
by the measurement of conversion electrons with an end window type 


G.M. counter. The absorption curve for electrons is given and its shape 
has been interpreted in the light of the results of ALPHONSE LAGASSE 
and JACQUELINE DOYEB (Compt. Rend., 239, (11), 670 (1954)). The value 
of conversion coefficient obtained in this way is 0.014 + 0.005. The 
discrepancy between the two values has been explained as due to K and 
L shell ionization by «-particles. No coincidences were found between the 
different parts of y-ray spectrum. 


1. — Introduction. 


Many workers (GRACE et al. (2), ZAJAC et al. (?), SIEGBAHN and SLATIS (8), ete.) 
have investigated the radiations from 2°Po. The following facts have been 
established: 

1) Ground level of 2!°Po decays predominantly through an x of 5.3 MeV 
to the ground state of 206PD. 

2) A small fraction (~10- %) of the ground state of 2!°Po decays through 
an a of 4.5 MeV to the excited level of 206Ph followed by a y-ray of 800 keV 
(GRACE et al. (1), DE BENEDETTI and MINTON (*)). 

3) A soft radiation of about 80 keV has also been observed. 


(*) Research Scholar of Ministry of Education, Government of India. 
(1) M. A. GRACE et al.: Proc. Phys. Soc., A 64, 493 (1951). 

(2) B. ZAJAC et al.: Proc. Phys. Soc., 60, 501 (1948). 

(3) K. SIEGBAHN and H. SLATIS: Nature, 159, 471 (1947). 

(4) S. DE BENEDETTI and G. H. Minton: Phys. Rev., 85, 944 (1952). 
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ZAJAC et al. (2) concluded by the critical absorption method that the soft 
radiation was of nuclear origin. GRACE et al. (*), on the other hand, consi- 
dered the soft radiation to be X-rays of 2°*Pb, by the same critical absorption 
method. These X-rays were interpreted by them as due to the conversion 
of 800 keV y-rays. Conversion coefficient obtained from their measurements 
of X-rays and conversion electrons was found to be 6.7%. GRACE et al. (1) 
suggested that 800 keV y-rays should be due to M, transition. DE BENEDETTI 
and Minton (*) measured the angular distribution between 800 keV y-rays 
and 4.5 MeV «-rays from which they concluded that 800 keV y-rays were 
due to £, transition, for which the Æ conversion coefficient should be 0.009. 
This latter conclusion also fits better with the theoretical considerations of 
the decay of 2?°Po (GOLDHABER and SUNYAR (*)). BARBER and HELM (5), who 
measured the total conversion coefficient to be 14%, suggested that a fraction 
of X-rays may be due to the A-shell ionization of «’s passing through the 
source. PRINGLE et al. (7) estimated the conversion coefficient to be 20 to 30%. 
The difference in the values of conversion coefficient with GRACE et al. (4) was 
explained as due to the A-shell ionization due to «’s. ALBURGER and FRIED- 
LANDER (5) measured the conversion electrons with a 6-ray Spectrograph and 
found the conversion coefficient to be 1 to 5%. 

In view of these controversial results about 80 keV radiation, a further 
study of the radiations from ?!*Po was undertaken. 


2. — Measurements by scintillation spectrometer. 


210Po, procured from Atomic Energy Commission Canada, was prepared 
by bombarding Bi,O; with slow neutrons from the reaction ?°*Bi(n, y)?!°Bi + 
> 20Po. Its activity at the time of measurement was 0.05 micro curie/milli- 
gram. A source of 0.184 gm/cm? was prepared on a perspex disc. The spectrum 
of the radiations was taken with a scintillation spectrometer using NaI(TL) 
crystal by placing the source right over the crystal. The spectrometer was 
calibrated for different energies with y-rays of 17Cs, 2*Hg and ‘*Ru. Enough 
absorber (perspex and aluminium) was placed between the source and crystal 
to absorb the æs and conversion electrons. The NaI(TL) crystal was one 
inch in diameter and half inch in thickness. The differential spectrum thus 
obtained is given in Fig. 1 (after subtracting the background) which clearly 
indicates the presence of two radiations, one at 76 keV and the other at 


5 


) M. GOLDHABER and A. W. Sunyar: Phys. Rev., 88, 906 (1951). 
) W. C. BARBER and R. H. Helm: Phys. Rev., 86, 275 (1952). 

) R. PRINGLE et al.: Phys. Rev., 87, 384 (1952). 

‘) D. ALBURGER and G. FRIEDLANDER: Phys. Rev., 84, 523 (1951). 
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800 keV. No other y-ray was noted between 20 to 150 keV. The small humps 
at about 0.27 and 0.40 MeV are due to back scattering, which was confirmed 


by their relative reduction for very thin source. The measurements of relative 


intensities of 76 keV and 800 keV 
radiations were carried out as 
The detection efficiency 
of the erystal for the two ra- 
diations was taken to be (1— 
— exp [— pd}) where y is the total 
absorption coefficient for the par- 
ticular radiation and «d» is the 
effective path length of the ra- 
diation through the crystal. The 
value of 


follows. 


«d» was obtained by 
actually drawing various paths 
in the crystal and taking their 
averages. The value of u was 
obtained from the graph of ab- 
sorption coefficient of the y-rays 
(K. SIEGBAHN (°)). 

For determining the total 
counts of 800 keV y-rays, their 
Compton distribution was drawn 
as shown in Fig. 1 by the dotted 
line. Knowing the resolution of 
the spectrometer at 300 keV, it 
is possible to get the resolution 
at 76 keV. From this the base of 
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Fig. 1. — The differential pulse spectrum of the 
y-rays and X-rays from Polonium-210. 


76 keV y-ray was determined. Compton distribution of 800 keV y-ray can then 


be drawn by smoothly joining the base to the Compton edge. 


Total counts for 76 keV can, of course, be determined directly by the area 


under the peak. 


Finally the absorption of these radiations in the aluminium covering and 


the source were also taken into account. 


Escape peak correction was not 


required because in both cases their contribution has been added to the 


total counts. 


Apart from statistics, various sources of error are: (a) uncertainty due to 
i ’ à 


the numerical integration used in determining «d ». 


estimated was about 10%. 


The maximum error 


(b) The uncertainty in drawing the total ditter- 
ential spectrum for the two radiations. After 


‘areful consideration the max- 


(9) K. SIEGBAHN: Beta-Gamma Ray Spectroscopy (1955). 
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imum error in conversion coefficient due to this uncertainty was estimated 
to be 10%, the final value of conversion coefficient, thus obtained, was 


0.039 + 0.01. 


3. Measurement of conversion electrons. 


A source of 6 mg/em? was prepared on 8 mg/cm? cellulose tape. The source 
was put right over the window of an end window type G.M. counter, with 
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Fig. 2. — Curve A, shows the coincidences between 
the A and Z X-rays and electrons for various thick- 
nesses of absorber between the source and the 
electron detector as reported by LaGasse and 
Doyen (1954). Curve B, shows the number of elec- 
trons for various thicknesses of absorber between 
the source and the detector as observed by the 
authors. 


the tape between the source 
and the window. The thick- 
ness of the window was 
about 3 mg/cm?. In this 
way 11 mg/cm? absorber is 
inserted between the source 
and the 
enough to check «-particles. 
were taken 
for various thicknesses of 


counter which is 
Counts, then, 


aluminium inserted between 
the tape and the window. 
An absorption curve was 
drawn which is shown in 
Fig. 2, curve B. The thick- 
nesses shown in the curve 
include half the thickness of 
the source, the thickness of 
the tape and window and 
the thickness of aluminium. 

Apparently, the graph 
can be divided into two 
parts; the first one extending 
up to about 100 mg/em? and 
the other extending beyond 
that. The first portion can 
be interpreted as due to the 


electrons from A, L, shell ionization by «°s and the second portion due to 


the genuine A and L conversion electrons. 
the following reasoning. 


This interpretation is based on 


After Migdal’s theoretical work on the A, L, shell ionization by æs 
(MIGDAL (!°)), it was considered probable that «°s of ?!°Po contribute to the 


(1°) A. MiGpAL: Journ. Phys. (U.S.S.R.), 4, 449 (1941). 
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X-rays and electrons through K and L ionization (BARBER and HELM (5 
PRINGLE et al. (7), LAGASSE and Doyen (!!), Roy and GOES (!°)). LAGASSE 
and DovEN (1!) measured the coincidences between the K X-rays and the 
electrons and also between the K and L X-rays and electrons, for various 
thicknesses of aluminium for selecting various energies of electrons. Their 
curves for the number of coincidences versus energy of electrons, give a con- 
tinuous spectrum (Fig. 2, curve A), resembling very much the first portion 
of our curve (Fig. 2, curve B). Like our curve they also have a steep portion 
near about 80 keV followed by a lesser steep curve extending up to 300 keV. 
These coincidences are not due to conversion electrons because, in that case, 
coincidences should follow the curve of MARSHALL and WARD (!). These 
coincidences have, therefore, been interpreted by them as between the K, 
L X-rays and electrons emitted by the process of K and L ionization due 
to «’s. Similar reasoning can also be applied to the first portion of our curve. 
The second portion of the curve, when extended back, is similar to the one 
expected of monoenergetic electrons of 712 keV. The number of counts at 
zero thickness from this extrapolated curve has been, therefore, taken as due 
to the net conversion electrons. In the measurement of conversion electrons, 
the back scattering effect for electrons, which is only due to the finite thick- 
ness of the source, was considered very carefully. It was estimated to be 
much less than one percent in our case. This number includes A as well as 
L conversion electrons. No detectable counts were obtained after intro- 
ducing 400 mg/cm? of aluminium between the counter and the source. This 
ensured that there was negligible contribution to the counts due to y-rays. 

The over-all efficiency of the G.M. counter was determined by the cali- 
brated 8 sources of ***Pa and 2%T1 obtained from Atomic Instrument Co., 
U.S.A. It was further checked by calculating the solid angle subtended by 
the source on the effective volume of the counter (BURTT (é)) the efficiency 
obtained from both these methods was 0.083 +0.009. The total number of 
--rays of this source were obtained by the same scintillation spectrometer 
as used previously. The source was placed right over the erystal and the pulses 
of all sizes were counted. It was ensured, however, that no æs or B's were 
counted. The geometrical and detection efficiency of the crystal was cal- 
culated as before. The contribution due to 76 keV radiation to the total counts 
and the corrections due to absorption were also taken into account. 

The chief errors in the measurements of the conversion coefficient with the 


(11) A. Lagasse and J. DovEN: Compt. Rend., 239 (11), 670 (1953). 

(2) R. R. Roy and M. L. Goes: Compt. Rend., 237, 1515 (1953). 

(13) J. MARSHALL and A. G. Warp: Can. Journ. Research, A 15, 39 (1937). 
(14) B. P. BURTT: Nucleonics, 5, 29 (1949). 
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above method were due to the uncertainty in the estimation of the over-all 
efficiency of the % counter and in the extrapolation in Fig. 2, curve B: 

The final value of the total conversion coefficient thus obtained was 
0.014 + 0.005. 


4. — Coincidence measurements. 


It was also attempted to detect coincidences between the various portions 
of the y-rays spectrum. For this two scintillation spectrometers in coincidence 
were used. The resolving time of the coincidence circuit was 2:10-?s. No 
genuine coincidences were observed between 800 keV and 76 keV radiations. 

7 


Also there were no coincidences between 76 keV and the rest of the spectrum. 
Similarly no self coincidences were observed for 800 keV y-rays. 


5. — Discussion. 


Within the experimental error the value of 0.014-0.005 for the total con- 
version coefficient obtained by the measurements of conversion electrons is 
same as expected for E, transition (Rose tables (*)). This is in agreement 
with the angular distribution results of DE BENEDETTI and MINTON (*). Com- 
paring this result with the value of 0.067 obtained by the conversion electron 
measurements by GRACE et al. (1), the discrepancy seems to be due to the 
fact that they also measured the electrons due to A and Z shell ionization 
by as as discussed in the measurements of conversion electrons. 

On the other hand our value of 0.039-+0.01 for the conversion coefficient 
obtained from X-ray measurements is nearly the same (within the experi- 
mental error), as that of GRACE et al. (1). The reason for this higher value 
compared to 0.014-0.005 seems to be due to the contribution of X-rays from 
the ionization by æs, which comes out to be 60% of the total intensity of 
X-rays. This is in qualitative agreement with the calculations of Roy and 
GOES (12) based on Migdal’s theory. 

The above results help to remove the ambiguity about the exact value of 
the conversion coefficient of 800 keV y-rays as well as about the reason for 
the discrepancy in the measurements as reported in the literature. 
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Col metodo della scintillazione si è misurato il coefficiente di conversione dei raggi + 
da 800keV del decadimento del 2!°Po, ottenendo 0.039 -L 0.01. N valore del coeñi- 
ciente di conversione è stato anche ottenuto misurando gli elettroni di conversione 
con un contatore G.M. a finestra terminale. Si dà la curva di assorbimento degli elet- 
troni e la sua forma é stata interpretata al lume dei risultati di ALHPONSE LAGASSE © 
JACQUELINE Dove (Compt. Rend., 239, (11), 670 (1954)). Il valore del coefficiente 
di conversione così ottenuto è 0.014 + 0.005. La diversità tra i due valori è stata 
attribuita alla ionizzazione degli strati A ed L dovuta alle particelle x. Non si sono tro- 
vate coincidenze tra le varie parti dello spettro dei y. 


(*) Traduzione a cura della Redazione. 
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Summary. — The results of measurements of ionization made in Ilford G-5 
and K-5 photoemulsions by means of a semiautomatic apparatus are 
reported. The experimental fluctuations are in good agreement with 
values derived by simple formulae of practical use. While all the models 
proposed as yet for the mechanism of track formation appear to be over- 
simplified, an assumption is proposed to explain the results. The different 
methods of measurement are compared and some suggestion is given. 


1. — Introduction. 


The measurement of ionization in nuclear emulsions is a very important 
tool for the interpretation of the events, as it 
cities of the involved particles. 


allows us to determine the velo- 
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A number of papers have been published on the subject (127) in recent 
years and several of them (132!) have described some semiautomatic apparatus 
intended to speed up the measurements. 

The attempts to outline a model for the formation of the tracks (2727) may 
appear to be less important. However, as we shall see, the solution of this 
problem is of relevance both for the adoption of a correct method of measu- 
rement. and for the calculation of fluctuations. 

In the present paper we give further experimental results on the measu- 
rement of ionization (Sect. 2) and on the fluctuations, which fit rather well 
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(2) A. H. Morris: Phil. Mag., 43, 533 (1952); Phys. Rev., 41, 423 (1953). 
) 
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(5) G. Kayas: Compt. Rend., 238, 2153 (1954); Journ. Phys. et Rad., 15, 34 (1954). 
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with the fluctuation formulae easily calculated by means of the Blatt method. 
(Sect. 3). A short survey is made of the models of track formation proposed. 


by different authors: all of them appear to be oversimplified, but a simple 
assumption on the growing of grains seems to fit well with the experimental 
results (Sect. 4). Finally, we compare the most important proposed methods 


of measurement and give some suggestion on these (Sect. 5). The most im- 


portant results are summarized in Sect. 6. 


2. — Experimental results: measurement of ionization. 


We will use the term «r-gaps » for indicating a gap of length > r. The 
experimental quantities are: 


(1) H(r) = number of r-gaps per unit length, 

(2) B =4H(0), 

(3) lr) = total length of r-gaps per unit length, 
(4) w(r) = Ur)/H(r) 


As we shall see at once, according to our results, H(r) is an exponential 
decreasing function for r larger than a value 7, of the order of (0.2 —0.4) um. 
Therefore, we can add to the already defined quantities this other: 


(5) g = absolute value of the slope of the straight line which gives 
log [H(r)] as a function of r, for r> 1. 


The bulk of the measurements were made in a Ilford G-5 stack, exposed 
to 140 MeV protons. Some measurements will also be reported made on slow 
protons in K-5 Hford emulsions and on one heavy track in a G-5 plate exposed 
to cosmic radiation. All the measurements were performed by means of a 
semiautomatic apparatus which is described elsewhere (2°). It gives, at the 
same time B, H(r) for 5 different values of >; and {(0); it allows to deduce (0) 
from the ratio ((0)/B as well as g, from the distribution of H(r). 

The «quantum» of measured lengths is 0.04 um. The cell length is 
t—40um The velocity of work is 1.2 mm per hour, included dead times. 
The stability of measurements was tested in various ways (2°). Measurements 


were made by 3 different observers. While the generale results were comple- 


tely coherent, some detectable, though rather small, differences between their 
results were observed; in order to avoid the use of 3 different calibration curves, 


at last only the measurements made by the most skilled observer were con- 


sidered. 
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The general results are exhibited in Fig. 1. 

The curves show the variation with the proton residual range R of (0), 
B, H(0.2), H(0.4) and W(0) (all figures in um). They support the general 
results already obtained by other authors Ce): 


100 
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Fig. 1. - Measurement of ionization parameters in underdeveloped G-5 emulsions. 


Abscissae: proton residual range. Ordinates: number of gaps or gap length in pm, on 
the same scale. All values reported for a cell {— 100 um. 


The known flat maximum of B begins at about 2 mm while in the FOWLER 
and PERKINS (27) work, it begins at about 4 mm and in the Rome work (42) 
it starts from — 10 mm. The Rome plates were very heavily developed, while 
those of the present work were rather underde- 
veloped, in order to discriminate between slow 301 60° 
deuterons and protons. 

Fig. 2 gives a selection of the measured 
H(r) curves for 5 different values of the 
proton range È, that is of the ionizing power. 
Fig. 3 gives similar results for K-5 emulsions, 


Fig. 2. — Integral gap length distribution for protons 

at various residual ranges. The values reported are 

for a cell {=40 um. Measurements on underveloped 
G-5 emulsions. 


9531 


710 


G. CORTINI, 


G. 


LUZZATTO, G. TOMASINI and A. MANFREDINI 


while Fig. 4 refers to a heavy track (assumed to be a boron or carbon mini- 
mum track) in G-5 emulsions. 
As it is seen in the G-5 emulsion, all H(r) curves turn out to be exponential 
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Fig. 3. — As Fis. 
urement on K-5 emulsion. 


LAIGH (#8). 


range À —3.2 mm. The results of (ii) and (iii) are reported in Fig. 5, together 


0.1 


2. Meas- 


only for r>7,=(0.3+0.4) um. For r<% the 
curves are always bent toward the r axis and are 
similar to each other. 

As this result does not agree with those by 
other authors (27), we tested it in order to ascertain 
that it could not be due to systematic errors intro- 
duced by our measuring apparatus. Three different: 
methods were used: 

(i) Direct measurement by means of a 
standard eyepiece scale. The results seemed to be 
very subjective, and dependent on the microscope 
magnification. 

(ii) Manual measurement by means of a filar 
micrometer. 

(iii) Manual measurement by means of a filar 
micrometer arranged with a set of spider webs spaced 
in geometrical progression as proposed by O'CEAL- 


All these measurements were made on proton tracks with residual 


tone 


Fig. 4. — Integral gap length 

distribution for a heavy track. 

Measurement on G-5 emulsion 
strongly developed. 


(3) Seo at, EL 


r) 


Fig. 5. — Gap length measured with 
different methods: + filar miero- 
meter displaced by hand; © O’Ceal- 
laigh method; @ our apparatus. 


Our grateful thanks are due to Dr. O’CEALLAIGH for having sup- 


plied to us a filar micrometric eyepiece. 
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with results obtained by means of our apparatus. The agreement between 
these measurements makes us confident that they are reliable. The different 
experimental methods can probably explain the difference between our results 
and F.P. ones. 

An interesting connection has been found between two of the measured 
parameters; for every value of À, we found the equation: 


(6) B= H(0) = g exp [— 279], 


fitted well, as already found by various authors (see (34) eq. (2) and (29) 
eq. (3.1 a)). 

Here 2y is the diameter of the developed grain. Table I gives the 2y values 
deduced from (6), by measurements at different proton ranges in G-5 and 
in K-5 emulsion. The values obtained by direct measurement and by the- 
heavy track results are also reported. 


Tage I. — Value of 2Y (um). 


G-5 K-5 


| Protons 

| From eq. (6) 

Residual range mm | | 
0.16 0.605 a | 


0.48 0.615 — | 
0.84 0.618 — 
Dadi 0.597 0.524 
She 0.618 — | 
| 4.7 0.634 0.543 | 
| 6.7 0.639 DE | 
9.2 0.633 — 
13.4 — 0.583 
Mean | 0.620 | 0.553 
| Direct measurement | 0.63 | 0.52 | 


Heavy track | | 
From eq. (6) | 0.605 | = | 


Direct measurement | 0.65 | — 


Note that eq. (6) would be an obvious consequence of a pure exponential 
H(r) distribution (see below eq. (18A)) while it cannot be deduced in any 
simple way from our bent H(r) curves. 
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3. — Experimental results: fluctuations. 


The important problem of the fluctuations of the measured parameters 
was recently investigated by LOVERA (2°) and by O'BRIEN (*). According to 
O’BRIEN’s definitions, one must consider separately subjective (S) and emul- 
sion processing (Æ) fluctuations, in order to determine the true track (D) 
fluctuations, that is the pure effects of the statistical process of track formation. 
E and S effects largely depend upon experimental conditions. If these are 
good, T effects are the most important ones. 

We have calculated the fluctuations as follows. In each track a number 
of groups of cells has been selected. The number m of cells in each group was 
so small (m=4--10) that within it the ionization was practically constant. 
For each group we measured, besides other, the number n= Bt of gaps per 
cell and the integral gap length «= lt, and we evaluated 


I 
M 

Sl 

| 


: RO =, 
(S) Hees a DS 
È 1 I È 
= = dh IA 
de zii ) 


By putting m —1 instead of m in the denominators, we take account of 
the fact that the deviations were taken with respect to the mean values (7) 
calculated for each particular group. Finally, the values (8) obtained from the 
particular tracks were combined with each other by calculating the quadra- 
tical mean, in order to obtain our best values of 4} and uw? at any given value 
of the range R. 

With this procedure, we compare measurements made on the same track 


and in immediate succession so that we minimize the Æ and $ effects. We 


think that the obtained values of the Blatt fluctuation parameters (8) cor- 
respond correctly to the pure 7 effects. If one calculates them directly from 
the general mean values, uw, and y, are increased by some 10%: 

In his work (8) O°BRIEN compares the experimental values of Hz and pu, 
with the Blatt predictions based on the model by Herz and Davis, and finds 
a good fit. As we shall see, it is reasonable to use, instead, the Blatt formulae 
derived from the O’Ceallaigh model. These give about the same results (with 


(2°) G. Lovera: Nuovo Cimento, 4, 1476 (1956). 
(50) B. J. O'BRIEN: Nuovo Cimento, 7, 314 (1958). 
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differences of a few units per cent—see also Appendix—) and are much less 
involved, so that they can be of practical use. In the particular but important 
case r=0 (the one investigated by O’BRIEN and by us) they are: 


(9) fig =1—4yB 
2 | 1 

Bas ae. 

(10) pi A (2y al 


The grain diameter 2y, as well as the total blob number B are easily mea- 
sured, while g (if not actually 
measured) can be derived from 
eq. (6). From (9) and (10) one 07 sl 

sb 


immediately obtains the standard 0. 
fluctuations $B and dl: 05 
04 

B I | 
(11) oe 2 ri 5 ob = int 0.3 
VN l VX 02 


01 02 O4 06081 2 4 6810 20 40 60 80 


where N is the total number of Fig. 6. — Experimental fluctuations of B 
counted blobs and Y the total compared with equation (8). 
measured gap length. 

The Fig. 6 and 7 show the experimental values of parameters (8) in com- 
parison with eqs. (9) and (10). For 4, there is good agreement, while for y, 
we note a small but almost constant difference. 

We shall see in Sect. 5 a pos- 
sible reason for this difference. 

In any case we think it useful 
to stress that (9) and (10) give 
fluctuations considerably smaller 
R(mm) than pure Poisson distributions, so 
that their use is definitely advan- 


pî (micron) 


Fig. 7. — Experimental fluctuations of 1 (0) tageous for the correct evaluation 
compared with equation (10). of experimental errors. 


4. — Model of track formation. 


Referring to the paper by BLATT (2°), we put: 


Q(t) — «survival probability », that is the probability that the distance 
between the centre points of adjacent activated AgBr crystals along 
a track exceeds 1; 


G(y)dy = the distribution of the final grain radii, after focussing. 


46 - Il Nuovo Cimento. 
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Two extreme models have been investigated in detail by various authors. 
In both of them G(y) is a 6-function; that is the radius y is kept constant. 
For Q(t) we have: 


Model « À »: 


(12) Q(t) = exp [— gt], (g = constant for a given ioniz- 
ation power of the particle). 


Model « B »: 


| Q(t) 0<t< fp; (6 =constant crystal distance 
(13) | =1-—p, DD; along the track; p= constant 
[ =(1—p)?, 28<t<38. for a given ionization power). 


The model « B », proposed by HERTZ and DAVIS, is derived from the as- 
sumption that the crystals are arrayed along the path of the particle in a 
regular lattice of definite lattice spacing. This assumption is exaggerated, 
because in the emulsion there cannot be any long-range order. On the other 
hand, model « A » completely ignores the « graininess » of the emulsion before 
the processing. A more realistic attempt to take into account the real situation 
is a modification of model « A », that is: 


Model «Cy»: 


UTI Ot =, 
(14) 


| 


exp [— g(t — @)] PE CL 


The models described include most of the assumptions made by different 
workers. The work by O’CEALLAIGH (21) can be reduced to model « C'», with 
minor modifications. The same is true for the work by FOWLER and PERKINS 
(see the integral curve b in Fig. 5 of their paper) who deduce the exponential 
character of Q(t) by considering in a very detailed way the mechanism of track 
formation. 

The model by HAPP, HULL and Morris (21) is similar to model « A » but 
these authors assume a bell distribution function for G(y) . The model by 
DELLA CORTE et al. (?) is similar to model « B », with the particular assump- 
tion 2y = and the same is true for the fluctuation results by LOVERA. 
et al. (2°). 

By means of his general method, BLATT gets from the different assumptions 
on Q(t) and G(y) the corresponding distribution formulae for the relevant ex- 
perimental quantities (1) to (3). 

To avoid confusion, it should be recalled that the symbols used by BLATT 
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in his paper are: 


t = cell length 


(15) RAY = tr), 
(16) x(t) = tl(r) — rn(t) , 
(17) WT) = Lin — wr) re. 


In eq. (17) the index B in w” has been put in order to distinguish the dif- 
ferent meaning of the symbol w in the BLATT paper and in eq. (4). The para- 
meter x means the total length per cell of r-gaps, computing for each gap only 
the length in excess of r. For r=0, æ—tt. The symbols n and x have al- 
ready been used by us in eqs. (7) and (8), with the meaning of n(0) and x(0). 

In the limit t>2y+r; t>1/9, that is for suitable cell length, the Blatt 
results are as follows: 

Model « À »: 


(184) H(r) = gexp|—g{(2 7) | = B exp |—9rl, 


1 
(19 A) w(r) = n +r, that is AAG) = 


The ionizing power is related to the value of g. 
Model « B»: 
1 
(18B) HU pp) 


(19B) WC) (1 + 
with / integer part of (2y+-7)/6. 
The ionizing power is related to the value of p. 


Model « C »: 
(180) H(r) = exp [ga][1 + gu}-*g exp [— g(2y + r)], 


C a n 
(190) IG ot th 


All the results (184) to (190) can be directly deduced from elementary 
considerations. The Blatt method allows one to derive the same parameters 
under more complicated assumptions, and to deduce in a simple way the 
fluctuation formulae (see Sect. 5). 


716 G. CORTINI; G. LUZZATTO, G. TOMASINI and A. MANFREDINI 


An extensive comparison of the models « A » and « B» with experimental 
results was performed by CASTAGNOLI, CORTINI and MANFREDINI (!*) and by 
FOWLER and PERKINS (2°). 

CASTAGNOLI et al. studied in the first place H(0) and w(0). 

Their results could not be fitted with eqs. (184) and (194). A good fit 
was possible with (18B) and (19B), by adjusting the two parameters f and y 
which are relevant in model B (while in model « A » the only parameter at 
disposal is y). 

On the other hand, FOWLER and PERKINS, studied the gap length distri- 
bution, that is the function H(r) for variable r. They found an exponential 
curve for the whole r range of interest and for any value of g. Besides this, 
they found empirically the equation (6). 

A general prediction of models « A » and « 0 » is that the limit w® of w(0) 
for increasing ionization is zero. This is contradicted by the results of ex- 
periment. FOWLER and PERKINS made the assumption that the experimental 
minimum value w® of w corresponds to the exponential distribution of the 
gaps between crystals, before processing, along the path of the charged par- 
ticles. This picture is completely coherent. 

However, we are now facing the difficulty of explaining the initial bend 
of the H curves shown in Fig. 3, 4 and 5. This bend cannot be explained by 
any one of the models considered as yet. 

Models « À » and « € » give a pure exponential curve. One can try to fit 
it with the model « B », by considering that the « graininess » effect will be 
important only on the short distances: the fluctuation in crystal locations 
and in y are expected to round off the discontinuities of equation (18B) as soon 
as [> 3-4, so giving a continuous exponential curve. For smaller /’ we 
can expect some type of bending. 

However, the calculations made by F.P. for fluctuating distances of crystals 
as well as those made by us for fluctuating y (see Blatt « modified Herz-Davis 
model ») show that this is not the case. The growing of crystals during pro- 
cessing is so large that no detectable effect is expected with reasonable values 
of the crystal distance 6 and of the grain radius y, as soon as the fluctuation 
is supposed appreciable. Besides, the small calculated bend can be made to 
go towards the r axis as well as in the opposite direction as a consequence 
of small variations of the chosen parameters. 

All in all, the detailed calculations on any known model necessarily forces 
us to admit that the experimental distribution ought to be purely exponential. 

The only way out of these difficulties is to admit that all the models are 
oversimplified, and that some new effect must be considered. 

The first attempt one can make is to consider some short-distance effect 
which involves only the short gaps. However, such an effect would not explain 
the very good fit that we find with eq. (6) as shown by Table I, for all our 
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measurements in G-5 as well as in K-5 emulsions [note that eq. (6) is a par- 
ticular case of eq. (184) |. 

This result shows that the bend cannot be due to the disappearance of 
short gaps. Otherwise, we would expect for H(r) 

(r>») much smaller values than those experimentally È log U(r) 
found, or we should be obliged to take for G-5 emul- 
sions 2y ~ 0.4 um, which is unacceptably small. 

In other words, the experimental curves cannot be 
deduced from the pure exponential function (see curve 
a, Fig. 8) by bending it near the origin (curve b), but 
by a long range effect which gives a curve like e. 

Adoption of model « 0» gives no practical advan- 
tage because the constant factor exp [xg]/(1+xg) in 
eq. (18C) turns out to be not much larger than Fig. 8. Explaining the 
unity and almost irrelevant for practical purposes for meaning of r, and 6p. 
all reasonable values of «. 

So, we go back to model « A », and propose the following tentative ex- 
planation. It is by no means true that the centres of developed grains neces- 
sarily coincide with the crystal centres. During processing, we have consi- 

— derable mechanical displacements. 


Bae poro O29 200 Bersted cous The growing in the direction of 
Duning processing OOOKX)OO0O () growing of grains unaffected compact AgBr crystals 
alert processing cof&oo developed. grains is more difficult than the inter- 
n : : enetration of adjacent sponge- 
Fig. 9. — A very schematic picture of the P en Pong 


erowing of grains under the assumptions like developed grains. This effect 
| specified in the text. tends to reduce the blob lengths 
and to increase the gap lengths 
(see Fig. 9), while it does not alter appreciably either H(0) or g. 
So, eq. (6) is still valid, while (184) must be substituted with 


(20) H = B exp [— g{r — 6(r)}] = g exp [— g{[2y — d(r)] +r}], 


where d(r) is an increasing function of r, which becomes a constant, 0) Say» 
for r > r, (see Fig. 8) 


(21) 


Eq. (194) is now valid only for r> 7%. 

The r, value is (0.3 —0.4) um for G-5 plates and (0.2 +0.3) um for K-5 plates. 
Eqs. (18A) and (20) fit with all our results, as well as with those by CASTA- 
GNOLI et al. and with the Herz results reported by BLATT (see (?°), page 266). 
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All the quoted results give 6,~ 0.20 um in reasonable agreement with our 
present value, 6,=0.24 um, for G-5 plates. 

For K-5 plates we have 6,~ 0.10 um, which shows that 6, is strongly de- 
pendent on the crystal diameter. Of course, the processing can produce some 
variation of do. 

The preceding considerations make it reasonable to use—at least as a first 
approximation—the Blatt fluctuation formulae for pure exponential distri- 
butions. 

These are (9) and (10) for measurements made with »=0. For the more 
general case r 0, we have 


(22) ur=1—2(2y +r)g exp [—g(2y +7)], 
(23) ua = ni 1—[1—g(2y + r)] exp [— g(2y + r)]}. 


Here 4, is the fluctuation parameter for the variable #,(t) (or t-!2,(t)) so 
{havi > 2, tor r—> 0: 

The fact that the experimental values of u, are somewhat smaller than 
given by eq. (10) is in qualitative agreement with the initial bend of the H(r) 
curves. 

While in general the preceding formulae are only approximate, in the case 
r>ry exact expressions of u}, and y? can be immediately derived from (22) 
and (23), by substituting 2y with 2y — 6): 


(22) p= 1—2[(2y — 0) + 11g exp [— g{(2y — à) +3]; 


(23') pè = TU [1 — gf{(2y — do) + r}] exp [— g{(2y + do) +7}]}. 


5. — Methods of measurement. 


On the basis of their results which give pure exponential H(r) curves, 
FOWLER and PERKINS proposed a method of measurement based on eqs. (6) 
and (184). For thin tracks they measure only B and deduce g by means of (6). 
For thicker tracks the value of B becomes insensitive to variations of ioni- 
zation. Then the authors measure B and H(r) for a single selected value of r, 
and deduce g by means of (184). The deduced values of g, if normalized to 
the value g* assumed by g for plateau relativistic particles, turn out to be 
independent on the processing of the plate (27). 

The method does not require any special apparatus. For thin tracks (pro- 
tons with & of the order of some cm) it is very simple and, as our results con- 
firm the validity of eq. (6), in this case the method is not put in question. 
On the other hand, the method for thicker tracks seems to be subject to some 
dependence on the experimental counting technique. 
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If a measuring machine is available, the choice of the method must be 
made from a balancing of the following factors: 


) Time of measurement and involved calculations. 
) Sensitivity of measured parameters on the ratio R/M of range to mass. 
(iii) Fluctuations of the considered parameter. 


Stability of results with variation of processing conditions and with 
location in the emulsion volume. 


All these factors must be studied in any particular case, following the 
concrete experimental situation, and the aim of any given research (e.g., measu- 
rement of mass or discrimination between particles of known masses, ete.). 
For the particular case of our plates, we decided between the parameters di- 
rectly measured, to use {(0) for R< 10mm and w(0) for R>10 mm. 

However, the exposed results show the convenience of completely elimi- 
nating the short gaps (r < ») from the computations, and to consider for each 
r-gap only the length in excess of r. This means to measure the variables 
H(r,) and 2, (t) (see eq. (16)), and w’(r,) (see eq. (17)). 

Indeed, for r> 7, we have a pure exponential curve, so that we can make 
direct measurements of g—[w/(r,)f! without any difficulty or complication. 
Besides, just as a consequence of the bent curves H(r), the difference between B 
and H(r) is very small, so that the loss of information is negligible. 

These conclusions are particularly true for the new emulsions Hford K-5, 
where the bend is very small. 

A minor modification of our apparatus can allow one to measure ,(t). 
This modification has already been made, and we plan to exploit it for further 
work. 


6. — Conclusions. 


Summarizing, our most important conclusions are as follows: 


1) The gap length distribution H(r) is exponential for r>», with 
ro= (0.3-0.4) um in G-5 emulsions and n= (0.2--0.3) um for K-5 emulsions. 
For r< n the curves show a strong bend toward the 7 axis. 

2) The simple Fowler-Perkins equation (6) is confirmed (see Table I). 

3) None of the models proposed as yet for the track formation process 
is satisfactory. However, the pure exponential model supplemented by a 
simple assumption on mechanical displacements of grains (see page 717) gives 
account for all our results. 

4) The statistical errors of measurement can be easily calculated using 
eqs. (11), by means of the first approximation Blatt-O’Ceallaigh formulae (9) 
and (10) for r —0 and (22) and (23) in the general case. For r > r, eqs. (22) 


Usi 
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and (23') give exact expressions of the fluctuations. Caution must be used 
in taking account of other effects (O’Brien Z- and S-effects see page 712). 

5) Among the parameters as yet measured, in our processing conditions 
the best ones are: 1(0) for residual range R<10 mm; w(0) for 10 mm< # <50 mm. 
For larger R the Fowler-Perkins method would probabily be the simplest and 
the most convenient one. 


6) For future work, we propose to use «, (t) and wi) = 9 


Our thanks are due to Prof. E. PANCINI for many useful discussions, to 
Ing R. SANNA for constant help in handling the measuring apparatus, to Miss 
BASIGLIO, and to Dr. L. CrurroLoTTI for having performed the measurements 
and a part of calculations. 


APPENDIX 


On the fluctuation formulae. 


Eq. (22), with its particular case for "= 0 [eq. (9)], holds for pure expo- 
nential H(r) distribution, and, also, holds approximately for any known 
model (2°). For instance, the Blatt-Herz-Davis formula reduces at once to (22) 
for (2y + r)/B=T + 3. 

Eq. (23), with is particular case (10), seems also to be valid for any model. 
For instance, the Blatt-Herz-Davis formula can be put, with some algebra, 
in the form: 


» 


2 “ 


ua = —{1—[1—g(2y +1)] exp [—g(2y + r)]} + 


Q 


2 gp 
aa g(2y + 7)] 1} tgh? 3 


The second term difference from eg. (23) is rather small: not more than 
~10% of the first one. 


RIASSUNTO 


Si riferiscono i risultati di misure di ionizzazione eseguite in lastre Ilford G-5 e K-5 
mediante un apparecchio semiautomatico. Le fluttuazioni sperimentali risultano in 
buon accordo con i valori desunti da semplici formule, di facile uso pratico. Mentre 
tutti i modelli proposti finora circa il meccanismo di formazione delle tracce appaiono 
semplicistici, se ne propone una modifica che spiega i risultati. Si discutono i metodi 
di misura e si danno in proposito alcuni suggerimenti. 
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The Coefficients of Friction and Wear. 


P. SEN 
National Physical Laboratory - New Delhi, India 


(ricevuto il 12 Maggio 1958) 


Summary. — The combination of the Savage surface energy and the 
Bowden and Tabor junction formation mechanisms of production of 
frictional energy is suggested. The physical description of the experi- 
ments is sufficiently simplified and represented by a stable state which 
envelopes the whole field of friction. This leads to a qualitative and 
satisfactory agreement with the observed results for the values of co- 
efficients of friction between solid surfaces, between solid surfaces with thin 
liquid or solid soft metal films, for wear, and for the correlation between 
friction and wear. The coefficient of wear is defined analogously to the 
coefficient of friction. 


1. — Introduction. 


While an explanation, from which a semiquantitative agreement is avail- 
able, can be found for a restricted set of experiments in the field of friction, 
no comprehensive theory which covers the whole field seems to exist, and 
there are to be found mutual contradictions between the explanations that 
are offered for the different sets of experiments. The most widely accepted 
theory of BowpEN and TABOR (1) according to which friction is due to welded 
junctions is substantiated by the experiments of BOWDEN and HUGHES (?) 
and FINCH (3) in which such junctions are found to occur in high vacuum or 
by the joining of freshly cleaved surfaces respectively. But FINCH (?) also 
pointed out that such a process would lead to much higher rate of wear than 

(1) F. P. Bowpex and D. TaBor: The Friction and Lubrication of Solids (Oxford, 


1950). 
(2) F. P. BowpEN and T. P. Hucnes: Proc. Roy. Soc., A172, 263 (1939). 


(3) G. I. FincH: Proc. Phys. Soc., B 63, 465 (1950). 
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the observed one and ARCHARD and Hirst (*) have shown that this discre- 
pancy is of the order of 102 for soft metals and 10* for hard metals. 

Then SAVAGE (5) showed that for thin liquid films of molecular dimensions 
between two sliding surfaces, under loads for which the Amonton’s law is 
obeyed, the frictional energy can be derived from the surface energy of the 
liquid film and the experiments of GORODETSKAJA and KABONOV (1) on the 
variation of the coefficient of friction for an electrolytic liquid film due to the 
interfacial potential confirm his conclusions. But no attempt seems to have 
been made to extend this mechanism to the friction from metallic films whose 
significance in the reduction of wear was pointed out by FINCH (*). 

The experiments on wear, for instance by BURWELL and STRANG (°), AR- 
CHARD (7), ARCHARD and Hirst (*), and KERRIDGE and LANCASTER (5) indi- 
cate that while light wear follows Amonton’s law and is essentially a surface 
phenomenon, no specific trends are found for heavy wear, although welded 
junctions, reminiscent of the Bowden and Tabor theory of friction, are ob- 
served and perhaps an appreciable layer of the sliding surface partakes in 
its occurrence. 

In order to obtain a connected description of the above mentioned sets 
of experiments we shall follow the general principles of FINCH (3) and con- 
sider the surface energy and the welded junction mechanisms of friction to- 
gether by combining the SAVAGE (5) and BowDEN and TABOR (!) processes 
of production of frictional energy. Then we can allow the former process, 
which produces little wear, to account for the major part of the frictional 
energy and reserve the Bowden and Tabor mechanism to account for wear, 
for higher than expected coefficients of friction, and for the experiments of 
BowDEN and TABOR (1) where the frictional shear is found to be equal to the 
shear strength of metals at room temperature. This combination of surface 
and junction mechanisms of friction is made by postulating that the tempe- 
rature of the contacting surfaces in the steady state of uniform motion against 
one another is in the neighbourhood of the melting point of the more fusible 
surface. Then the area of contact becomes much larger than that obtained 
by BowDEx and TABOR (°°) from measurements of electrical conductivity 
across the contacting surfaces at rest and at room temperature. We note that 


(*) J. F. ArCHARD and W. Hirst: Proc. Roy. Soc., A 236, 397 (1956). 

(5) R. H. SAVAGE: Journ. Appl. Phys., 19, 1 (1948). 

(6) J. T. BurweLL and C. D. STRANG: Journ. Appl. Phys., 23, 18 (1952). 
(7) J. F. ARCHARD: Journ. Appl. Phys., 24, 981 (1953). 


(°) M. KERRIDGE and J. K. LANCASTER: Proc. Roy. Soc., A 236, 250 (1956). Other 
references also are given here. 

(*) F. P. BOWDEN and D. TaBor: Proc. Roy. Soc., A169, 391 (1939). 

(1°) R. HoLm: Electrical Contacts (Uppsala, 1946). 
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a larger area of contact has been obtained by ARCHARD (!!) in the direct meas- 
urements of the true area of contact between surfaces in motion with respect 
to each other, but an alternate explanation, which is in agreement with the 
Bowden and Tabor hypothesis is also possible in this case and has been 
given by him. This increase in the area of contact makes possible a corres- 
ponding decrease in the strength of the forces that cause friction and there- 
fore surface tension replaces welding and becomes the principal mechanism 
for production of friction. 


2. The surface energy mechanism of friction. 


To obtain the frictional energy consider a hollow cylinder which may con- 
sist of several separate components, whose inner surface 0, is in contact with 
the outer surface o, of a rotating solid cylinder 2 and let their smoothed 
sliding surfaces be formed on the average of minute hemispherical asperities. 
Due to their small curvature even a small load will deform the contacting 
asperities plastically until the area of contact between the surfaces is sufficiently 
enlarged to reduce the body stress to its yield point value. When the inner 
cylinder revolves the area of contact will change and the plastic deformation 
shall reoccur until, in the stable state, the area of contact becomes constant 
and supports the yield stress. Due to the continuous supply of frictional 
energy the contacting surfaces will be heated and we shall assume that their 
temperature is in the neighbourhood of the melting point of the more fusible 
surface, say o,. Then the area of true contact is given by (1) 


(1) o(1 


2) = GJY(L D), 


where G is the load applied at the sliding surfaces and Y(1,L,) is the yield 
stress of the substance 1 in the neighbourhood of its melting point L, °C. 
Since the shear stress s,(1, Z,) is equal to the surface energy y(1, L,) in the 
plastic state and under the yield pressure, using von Mise’s plasticity criterion (1) 


CI IE 


= 
bo 
ene 


where s,(1, T), s,(1, 7) and Y(1, 7) are the shear, body, and yield stresses 
respectively at the temperature 7 °C, we obtain the minimum value of 


(1) J. F. ArcHaRD: Proc. Roy. Soc., A 243, 190 (1958). 
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(3) Y(1, L,) & 2y(1, La). 


A surface may be considered to be a regular array of molecules a few 
angstroms apart and as the surfaces in true contact over a small region slide 
over one another the interactions between their molecules will be broken 
and reestablished periodically. By supplying the binding energy the surfaces 
will be separated from their ground state to their excited states which in 
turn, when they come in contact again, will fall back to their ground state 
and release the binding energy, which will be carried away in the form of 
heat energy. The energy needed for these operations is supplied by the force 
that produces the sliding and is a part of the total frictional energy supplied 
to the system. 

To obtain its estimate we note that the excited states can be either the 
« vapourized » states where the molten surface is separated from the solid 
surface and is independent of it or the «transferred » states where a mono- 
layer of the molten surface o, is torn from its bulk surface and sticks to the 
surface o,. The energies necessary for forming these states are 


(4) E,(1|2) = y(1, L,)(1 + cos 8(1|2)) (12), 
(5) E,(1]2) = 2y(1, Zy)o(1|2), 


respectively, where the equation (4) is the Dupré relation (ADAMS (!2)) and 
9(1|2) is the angle of contact between the surfaces o, and o, at L, °C. A sta- 
tistical estimate of the ratio 7 of the probability of the formation of the trans- 
ferred states and the probability of the formation of the vapourized states is 


(6) n = exp [(E,(1|2) — E;(1[2))/2.(1]2)] + E.(1|2)/E.(1|2), 


and is in accord with KERRIDGE and LANCASTERS (5) observed magnitude of 
less than unity. 


Then, in the equilibrium state, the frictional energy is 


F(1|2) = 2.(1|2)/1 + 9) + E(1|2)9/1 +), 


= 4o(1|2)y(1, L)(1 + cos 9(1]2))(3 + cos9(1|2)) 1, 


(12) N. K. Apams: The Physics and Chemistry of Surfaces (Oxford, 1941). 
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and the surface coefficient of friction w,(1|2) is, by definition, 
(8) (1 |2) = F,(1|2)/@ & 2(1 + cos 0(1|2))(3 + cos #(1|2))-2. 


It will have the maximum value of 1 when the angle of contact is zero or 
there is complete wetting and will have the minimum value 0 for 9=x 
which however is not found. We note that perhaps one of the physical modi- 
fications that occur when the Beilby layer is formed or when the hardness 
of the sliding surfaces is increased is the reduction of the wetting. The ob- 
served values of u(1|2) are generally of the order unity for metals with oxide 
layers, of the order 20 for pure soft metals, and of the order 5 for pure hard 
metals. Thus while we find qualitative agreement for metals with oxide layers, 
we shall explain the higher values in Sect. 4 by showing that then there is 
an appreciable probability for the formation of junctions. 

Now in accordance with the observations of KERRIDGE and LANCASTER (8) 
that light wear of surface o, occurs after it has been transferred to the sur- 
face o, we note that the transferred surface may also form a drop and detach 
itself and cause wear. Let K(1) be the probability of the formation of the 
drop. Then the number of molecules of wear per unit sliding distance is given by 


(9) W,(1|2) = M(1) K(1)y0(1]2) = GN) K(1)n/ VC, La), 


where N(1) is the molecular density of the surface o,. Defining the coefficient 
of wear w,(1|2) analogously to the coefficient of friction as the wear per unit 
load we obtain 


(10) o(1|2) = N1)K(1)m/Y(1, Li). 


KERRIDGE and LANCASTER (*) have estimated that the wear particles corres- 
pond to about six transferred layers and in that case K(1) & nf, which is 
perhaps not an unreasonable estimate. 


3. — Thin films. 


The mechanism of friction when there exists a thin liquid film, a few mole- 
cules thick, between two sliding surfaces has been given by SAVAGE (°) and 
indeed the mechanism described in Sect. 2 for friction between two solid 
sliding surfaces corresponds to it. But now while the true area of contact 
remains the same as before, to produce sliding, work must now be done to 
separate the solid and the liquid film surfaces. Let p,(1|3|2) denote the 


Ue) 
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coefficient of friction between the surfaces o, and o, in the presence of a thin 
film o,. Then we expect 


2) y(3,L,) (1 + cos D 
Bo(1|2)  y(1, Li)(1+4 cos O 


bet | CO 


where (3/2) is the angle of contact between the surfaces 03 and o, and the 
surface energy of the film o, is to be taken at the temperature L, °C and 
yield pressure Y(1, L,). Furthermore 


(12) Hall (3 12/1 (114]2)  y(3, La)/y(4, La) - 


and these ratios are in general agreement with the observed ratios. 

We can now extend the same mechanism to thin soft metal films between 
two hard sliding surfaces and we note that under loads at which Amonton’s 
law is obeyed 


tta(Cu|Pb|steel) _ y(Pb, Leu) _ 350 


de; uo(Cu|steel) © , 


cei 
y(Cu, Leu.) 1100 i 
while the observed values are 


po(Cu|Pb|steel)  .2 


14 - RE = 
ue) (Cu | steel) 8 


Again since y(steel|In|steel) is .04 and u(steel|steel) is .4 the expected value 
of y(In|Z,,..) is 150 dynes/em at L... and Y(steel, L 


steel par 


4. — Junction formations. 


Besides the surface excitation the formation and rupture of welded junctions 
contributes to frictional energy and wear. Junctions can be formed due to 
adhesion of freshly cleaved surfaces or to permanent plastic deformation due 
to excessive loading or to adhesion of pure metals in high vacuum, and we 
note that such junctions are found in every frictional process. A rupture of 
a junction will generally occur for the softer substance and require a yield 
shear at a temperature much below the temperature of the sliding surfaces 
and which may be approximately put equal to the room temperature and 
is denoted here by R °C. Although the thickness of a few molecular layers 
only is involved here, we may call it the «body » mechanism of friction in 
comparison with the surface mechanism described earlier and we can attri- 
bute the major portion of wear to it. 


2548 


NI 
bo 
| 


THE COEFFICIENTS OF FRICTION AND WEAR 


Let us consider the junctions formed by adhesion at low loads and let 
£,(1|2) be the probability of an area of contact to be an area of adhesion. 
Then the total frictional energy and the total wear per unit sliding distance are 


P(1|2) = F,(1|2) + F.(1|2) 
(15) = o(1|2) [4y(1, L.)(1 + cos #(1|2))- 
-(3+ cos #)1|2))-1(1 — &,(1|2)) + 8,0, B)E,(1|2)] , 


and 
(16) W(112) = W.(1]2) + W.(112) = 
= N(1) c(1|2)[E(1)m(1— €,(1|2)) +7(1)£(1]2)]; 


respectively. Here J(1) are the number of molecules which are detached per 
molecular contact, and we have assumed that the area of tear and the area 
of adhesion are the same. Furthermore if the adhesive areas are assumed to 
follow the laws of plasticity and Amonton’s law holds, we obtain 


(17) u(1|2) = 2-*[4(1 + cos &(1|2))(3 + cos 0(1|2))-!(1— €,(1|2)) + 
+ 8,(1, R)E (1|2)/y(1, La)] ; 


(18) o(1|2) = N(1)/Y(1, L)[EW)m(1-- E(L|2) +71)£(1)2)]. 


We note that for s (1, R)/y(1, La) & 10°, 2-104 and w(1|2) 20, 5 respectively 
we should have é£,(1|2) & 4:10-?, 4:10-*, which is in qualitative agreement 
with the observations of ARCHARD (7) and ARCHARD and Hirst (*). 

Adhesion can be considered to be due to the penetration of a rectangular 
potential barrier, whose width is a and height is equal to the yield stress 
Y(2, R) of the body at the room temperature R °C, by the molecules of the 
molten surface o, whose energy is represented by the stress s(1, Z,). Then 
,(1|2) is the probability of the penetration of the potential barrier and we 
find (7%) 


te ee Payee LL 


(19) (112) = 16 Y, R) 


where B is a constant. The quantities K(1), J(1) and a that occur in the 
relations (17) and (18) depend on the physical and chemical structure of the 
surfaces and to determine their nature we need more experimental results. 


(8) L. I. Sonirr: Quantum Mechanics (New York, 1949). 
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Generally since there is little surface wear we may expect K(1) to be small 
and we expect J(1) to be of the order unity. Then to explain the variation 
of &,(1|2) from 10-2 to 10-* for soft and hard metals we expect the width and 
height of the potential barrier to increase with the hardness of the surface. 

We note that a chemisorbed oxide layer is generally found on the surfaces 
of most metals. While the surface friction and wear are not appreciably 
atfected by it the probability for the formation of junctions is found to be 
much less and this may be attributed to the increase in the height and width 
of the repulsive potential barrier which must be penetrated for intermetallic 
contacts to occur. Thus under these conditions the coefficients of friction are 
of the order of unity or less and furthermore, as pointed out by FINCH (*), 
since the oxide layer which is rubbed away by the surface wear is immediately 
reformed due to the presence of oxygen in the atmosphere, a stable state 
where the friction is of the order unity and the wear is low, exists for the 
metals that are used in bearing engineering. 

By applying high loads BOWDEN and HUGHES (?) have obtained values 
of &,(1|2) of the order of unity and we shall now determine the coefficient of 
friction under these conditions. Let 2; be the total area of surface o, and let 
G, be the load such that 


(20) Zi YO, L) = Go. 
Then 
Veal) for. GC, 
(21) s,(1, L) = 


GYD) oe 


From (19) and (21) we note that while €(1|2) remains constant for @ < G, 
it is proportional to G for G> G, and we may write 

ert?) Lov! G= Gy, 

GE(1|2)/G, fontG = Gs 


Then for G> G, 


(23) w(1|2) = @o/(2@) [4(1 + cos 9(1|2))(3 + cos 9(1|2))-1(1 — G£,(1|2)/&) + 


+ Gs,(1, R)E,(1|2)/(Gy(1, La))], 
and 


(24) @(1]2) = GN(1)/(GY(1, L)) [K(1)y(1 — G£,(1]2)/G) +J1(1)GE,(1/2)/&], 
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and we note that under these simple assumptions while the contributions of 
the surface term become negligible the contributions of the body term remain 
the same as for (17) and (18) and Amonton’s law is followed. 


‘5. — Conclusion. 


To determine the nature of the mechanism which produces frictional energy 
the physical description of the experiments is much simplified and is formu- 
lated in terms of a stable state that envelopes the whole field of friction. 
Then the surface process of SAVAGE (5) and the junction process of BOWDEN 
and TABOR (') for the production of frictional energy are combined by pos- 
tulating that the areas of true contact are in a plastic and almost molten state 
under an yield stress and at a temperature in the neighbourhood of the melting 
point of the more fusible surface. We then find that a qualitative and satis- 
factory agreement is obtained with the observed results for the values of the 
coefficients of friction between solid surfaces, between solid surfaces with thin 
liquid or solid soft metal films, and for the correlation between the coefficients 
of friction and wear. 

To obtain more quantitative results some of the approximations made in 
order to get a simple statement of the problem will have to be discarded and 
we should like to enumerate the assumptions that have been made here. 
Frictional energy is generally accompanied by wear and therefore the area 
of contact of the sliding surfaces is constantly changing, specially when there 
is heavy wear, but we have assumed a constant area of contact. Then all 
the physical, structural and chemical properties of the sliding surfaces and 
of their interaction with each other have been collected together in their sur- 
face and intersurface energies. The physical and chemical reactions at the 
surface such as the formation of Beilby layer and oxide layer have not been 
fully treated although it is known that they influence much the coefficients 
of friction and wear. Then it is assumed that the temperature at the sliding 
surface is the melting point of the more fusible surface L, °C, but perhaps 
a lower temperature is found which can influence the results appreciably. 
The conductivity of the substances and the rate of supply of the thermal 
energy, which will determine the body and surface temperatures and the shear 
stresses, have also been neglected. Then the wear parameters K(1) and spe- 
cially J(1) are not fully defined and to do so a more explicit distinction between 
the body and surface layers in our model seems necessary. 


47 - Il Nuovo Cimento. 
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This work was done at the National Chemical Laboratory, Poona, and the 
author is grateful to Professor G. I. FincH for kind hospitality, for suggesting 
the problem and for helpful discussions. 


RIASSUNTO (*) 


Si propone la combinazione dell'energia di superficie di Savage col meccanismo 
di Bowden e Tabor della formazione di giunti per la produzione dell’energia di attrito. 
La descrizione fisica degli esperimenti risulta sufficientemente semplificata e rappre- 
sentata da uno stato stabile che avviluppa l’intero campo di attrito. Ciò porta ad un 
soddisfacente accordo qualitativo coi valori osservati dei valori dei coefficienti di attrito 
tra superficie solide e tra superficie solide con interposte sottili pellicole di liquidi o 
di metalli teneri, per quanto riguarda l’usura e la correlazione fra attrito e usura. Si 
definisce il coefficiente di usura in analogia al coefficiente d’attrito. 


(*) Traduzione a cura della Redazione, 


LETTERE ALLA REDAZIONE 


(La responsabilità scientifica degli scritti inseriti in questa rubrica è completamente lasciata 
dalla Direzione del periodico ai singoli autori) 


On Convergence to the Moller Wave Operators (*). 


M. N. Hack 


Argonne National Laboratory - Lemont, Illinois 


(Ricevuto il 23 Maggio 1958) 


It has been recognized by physicists for some time that the proper formulation 
of time-dependent scattering theory is in terms of wave packets, i.e., normalizable 
state vectors. Recently it was shown by J. M. Cook (!), by means of rigorous 
methods of the theory of linear operators in Hilbert space, that for potentials of 
integrable square 


(1) | | | Vir)? dady dz < co 


the limits of the U(0,¢) operators as t+ + co exist strongly and are in general, 
i.e., when the total Hamiltonian has bound states, properly isometric (defined every- 
where and norm-preserving, but non-unitary) (?). 

From the standpoint of the stationary-state formulation of scattering theory it 
is expected, and it will be proved here, that the results of Cook remain valid for 
potentials which fall off at large r faster than 1/r, and which are required to be 
square integrable only in the finite domain. That is, we impose the following con- 


ditions on the potential: 


(2) Vir)?dr < © for all finite R 
Irl<R 
and 
(3) in) = O0 er), as r-= oo, for some x> 1, 


4 


(*) Work performed under the auspices of the U.S. Atomic Energy Commission. 

(1) J. M. Cook: Journ. Math and Phys., 36, 82 (1957). 

?) M. RoSENBLUM: Pacific Journ. Math., 7, 997 (1957) has proved the following theorem. 
Let À and B be self-adjoint operators with weakly absolutely continuous spectral measures and 
let the « potential» P= B— A be a completely continuous operator such that the trace of |P | is 
finite. Then the limits as # — + co of U(0,{)= exp [/Bt] exp [— i4t] exist strongly and unitarily 
transform A into B. 


=] 
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For ç in a certain subset © of the Hilbert space © of all square integrable func- 
tions of r=(«, y, 2) it is shown in reference (1) that 


t 


(4) U(s,t)p= g + i exp tirs] dr exp [— iH(s —t)]V exp[—1Kt]p 


. 
s 


and that the linear manifold 9% spanned by 6 is dense in . Here A is the closure 
of — V2, i.e., it is the kinetic energy operator in suitable units, H—K+V, and 


U(s, t) = exp [ts] exp [— iH(s —t)] exp [At]. 


The set 6 consists of those elements y of whose Fourier transforms are functions 
of the form g(p)= exp [—a?(p—pp)?] with a? > 0. 
From the convergence of the integral 


co 


(5) | IV exp [— iKrt]p|] dr < co 


V 
— o 


it follows (*) for arbitrary e > 0 that 


U(0, thp — U(0, s)p| = || U(s, ep — | = 


t t 


= | | exp[—tH(s—t)|V exp [—iktlpdt| < | |Vexp [—iKt]oldr<e, 


. 
s 


for s and ¢ sufficiently large, positively or negatively. Consequently the strong 
limits of U(0,t)g as t— + co exist by virtue of the completeness of $, first of all 
for g in G and then by the density of 9% for all g in $. The isometry of the limit 
operators follows from the strong convergence and the unitarity of U(0,4). 

In this note, Equation (5) is established for potentials satisfying conditions (2) 
and (3). The existence of the strong limits and their proper isometry then follows 
from the argument of Cook (3). From the relation (4) 


exp [iKt]p(r) |? = (4(at +12) 3 exp [— ele —b}], 


where @ is in 6 and 
aie 


= 2 ; pre — el 
2 ‘Po i Data 


we have 


|V exp liÆt]e|?= | V(r)? exp [iKtjg(r)|?dr = (4(a4+t2))-3I, 


… 


(*) We remark that the strong continuity of V exp [— iKt]o as a function of { which is used 
in the proof of Eq. (4) holds, for 9 in 6, for potentials satisfying (2) and (3) by dominated con- 
vergence. We note also that 77 is self-adjoint in the present case, since potentials satisfying (2) and 
(3) fulfil the conditions of Karo: Trans. Amer. Math. Soc., 70, 195 (1951). 

(*) An inconsequential omission of a factor of 2 in ref. (4) has been taken into account in the 
formula for e. 


ON CONVERGENCE TO THE MOLLER WAVE OPERATORS 


~I 
w 
vs 


where 


7 


If [exp [—e(r — b)? Vir)? dr 


and the integrations are extended over all space. By (3) there exists an PR, > 0 
such that 


for r> i: 
Without loss of generality we can assume x< 3. For r > 2b, 
2 x PASS +2 
(r — b)? > (r —b)?>Hhr?. 


Therefore, letting R= max.{2b, R,}, 


2 


dr<J + yeta-3) 


pee 


1 [vepar al | exp [— }er?] 


r<R r>R 


where y is a constant independent of ¢ and 


n 


Je | VWAGAP lip 


r=R 
Suppose now that 2b > À, and set 
n= | V(r)? dr . 
r=R, 
Then 
r Î vi È M? DI 
J =|V(r)?dr +] V(r)?dr<J,+ | — dr = J,+ pb?-* . 
| ons 
v<R1 R,<r<2b r< 2b 


This inequality holds also when 2b< /,, since then J=J,. Thus 
| V exp [iKt]p|?< (4(a4 + t)) #(T,+ ub3-2x + ver (20- 9) E 


Recalling the expressions for b and c in terms of {, we see that 


ie 
| V exp [iKt]@]|= o +] | as |t|= co. 
Therefore 
V exp [iKt]g — 0 as t > + co 
and 
fir exp [— iKt]p||dt < co. 
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Branching Ratio of ThC (?!?Bi). 


D. PROSPERI 


C.N.R.N., Laboratorio Fisico della Divisone Geomineraria - Roma 


S. SCIUTI 


Istituto Nazionale di Fisica Nucleare - Sezione di Roma 


(ricevuto il 27 Giugno 1958) 


A knowledge of the branching ratio branching ratio are due to KOVARIK and 
of ThC is of considerable interest both ADAMS (!) who determined it by examin- 
for the nuclear spectroscopy of heavy ele- ing the « activities of ThC and ThC’ 


2800 


2400 Thc’ 


2000 


1600F 


Counts per channel 


1200 


Channels 


ments and for applicative questions con- are 

cerning the Thorium radioactive family. i . È 
sa - è () A. Kovarik and N. ADAMS: Phys. Rev., 
rhe first reliable measurements of the 34, 413 (1938). 


BRANCHING RATIO OF THC (??Bi) 


using a plane ionization chamber and 
layers of old Thorium oxide. 

Recently, SENFTLE, FARLEY and LaA- 
ZAR (7) have carried out new measure- 
ments using a proportional counter and 
very thin layers of ThB in equilibrium 
with ThC and ThC’. The results ob- 
tained by these authors are not in com- 
plete agreement with those obtained by 
Kovarik and ADAMS. 

For this reason we have undertaken 
a precise redetermination of the branch- 
ing ratio of ThC. 

A 2a gridded ionization chamber (3) 
connected to a 100 channel pulse height 
analyzer was used in the work. The 
technique used for preparing the source 
was analogous to that described by 
SENFTLE and coworkers. 

Nickel dises of 10 em diameter were 
suspended over about 500g of aged 
Thorium nitrate contained in a metal 
vessel. The Thoron decay products de- 
posited themselves spontaneously upon 
the dises which were kept at — 150 V with 


@) PE. E: SENFILE, T. A. FARLEY and 
N. Lazar: Phys. Rev., 104, 1629 (1956). 
(3) D. PROSPERI and S. SCIUTI: Precise 


Radiometric determination of Uranium and Tho- 
rium contents in Ores. Communication presented 
at Geneve Conference, September 1958. 


1 
w 
Or 


respect to the vessel. The time required 
for the deposition was about 20-hours. 
Each dise was then placed in the chamber 
which permitted the measurements of 
the « activities of ThC and ThCO' 

A typical pulse-height spectrum is 
given in the figure; the width of two 
x lines is about 2% at half height. In 
the table are given the final results of 
the above mentioned authors and our 
value which was obtained by taking the 
average of nine independent determin- 
ations. 


a x 
branching ratio —— 

> > 
+6 


| 
| Authors 
| 


KOVARIK 
and ADAMS 


0.362 + 0.006 


SENFTLE 
and coworkers | 


0.352 + 0.003 


Present work | 


As one can see from the table we 
obtained for the branching ratio a value 
which lies between those previously 
found. The expected over-all uncertainty 
is less than 0.9%. 
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Photoneutrons from Au. 


S. CAVALLARO, V. Emma, C. MILONE and A. RUBBINO 


Istituto di Fisica dell Università - 


Catania 


Centro Siciliano di Fisica Nucleare - Catania 


(ricevuto il 30 Giugno 1958) 


In order to extend the experimental 
knowledge on the photonuclear inter- 
action, the energy spectrum of the 
photoneutrons from Au has been studied 
with the experimental method described 
in a previous paper (1). 

A 60¢ Au target was irradiated by 
a collimated 30 MeV bremsstrahlung 
beam from the B.B.C. Turin betatron. 
The photoneutrons emitted at angles of 
90° with the photon beam were detected 
by means of the proton recoil tracks in 
Ilford C2 200 um thick emulsions. The 
plates were screened against the neutrons 
coming directly from the betatron and 
from the low energy photons coming 
from the target. 

Normally the plates were scanned 
with 55x oil immersion Koristka objec- 
tives and 8x eyepieces: («slow scan- 
ning »). Only for very short tracks the 
length was measured with a 100 X objec- 
tive. Following the conventional me- 
thod (*?), only the recoil tracks were 
considered which come to rest inside 
the emulsion and had an original di- 
rection into a «square piramid » of 15° 


(1) G. 'CORTINI, C. MILONE, A. RUBBINO 
and F. FERRERO: Nuovo Cimento (1958) in press. 

OAMAROSENCENMcleones MA AT PES 
38 (1953). 


semiaperture centred in the neutron 
beam direction. 

An additional «fast scanning » 
— gee (1) — was made with a 30x 


objective in order to improve the sta- 
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Fig. 1. — Energy spectrum of the photoneutrons 
from Au emitted at ~90° with the y-ray beam. 

= Experimental figures. For the high energy 
tail figures are multiplied by 10. The curve (a) 
is the evaporation spectrum calculated under 
the assumption indicated in the text. Curve (à) 
is a spectrum calculated for direct emission under 
the hypothesis that the residual nucleus after 
the neutron emission is left in the ground state. 

Y*= neutrons//?/0.5 MeV/mole/A2. 


PHOTONEUTRONS FROM Au 


tistics over 3.5 MeV neutron energy: the 
track lenght measurements were made 
with a 55x objective. 
Exposure and scanning 
summarized in Table I. 


data are 


agio Le 


Element | AU 
Weight (g) 60 
(y.n) threshold (MeV) 7.9 
ÆE max (MeV) | 30 


Dose (Roentgen) 


Slow scanning (p<15°) | 


Scanned volume | 54 mm? 


Number of tracks 


E,<5 MeV 774 
| Number of tracks 
E,>5 MeV [MON 45 
| Supplementary scanning 16 mm? | 
| with 100 for low | 
| energy | 
| Number of tracks | 
| E,<4 MeV | 186 
| Fast scanning 
| Scanned volume 210 mm? 
Number of tracks 
(p< 15°) | 
| 3.5 MeV<E,<5MeV | 62 
| E,>5 MeV | 145 
Number of tracks 
(15°<p<30°) | 
TEATINO | 135 


Taking into account the appropriate 
factors (1), from the experimental proton 
Spectrum the Au neutron spectrum shown 
in Fig. 1 is obtained. The background 
was 3%. 

The energy distribution may be 
compared in the conventional way with 
the evaporation theory (*). The expected 


(3) V. F. WeriIsSKOPF and D. H. Ewine: 
Phys. Rev., 57, 472 (1940). 


V3.0 


differential neutron energy spectrum 
EB) is (3) 
HE.) = KE;0(E,;)" 
#6 
| 0(Ex) 0,n(E,)1(E,, Es) dE, 
Ey Bn 
Thee Î E,o(E,)0(E,) dE, 
0 
where: 
E, — neutron energy; 
o(E,) — reaction cross section for 
neutrons of energy Æ, on 
residual nucleus (4); 
Be = binding energy 
(Ex) = energy level density in resi-- 
dual nucleus with excitation 
energy H,= H,—B,—E,; 
o(y,n) = cross section for the (+, n) 


process (°); 
I(E,, Eg) =bremsstrahlung spectrum with 
maximum energy Lg (57). 


Assuming for ©(#,) the form 


OG, OX Vol ; 


curve (a) is obtained for x=1,6(A-40)? (8). 
The experimental spectrum 
two points of disagreement with respect 
to the calculated one. 
First, we note a high energy tail 


shows. 


(*) As Prof. G. CorTINI pointed out to us, 
the integral in the denominator of this formula. 
cannot be ignored as was mistakenly done in 
a preceeding paper (*). Therefore the evaporation 
curves of fig. 2 and 3 of the quoted paper cor- 
respond to density level distributions rather- 
lower than indicated there (3.1). This gives no 
changement in the conclusions of paper (!). 

(4) H. FESHBACH and V. F. WEISSKOPF: 
Phys. Rev., 76, 1550 (1949). 

(®) R. MONTALS3ETTI, L. Katz and J. GoL- 
DEMBERG: Phys. Rev., 91, 659 (1953). 

(5) I. Katz and A. G. W. CAMERON: Can. 
Journ. Phys., 39, 437 (1957). 

(?) €. MILONE, S. MILONE-TAMBURINO, 
R. RINZIVILLO, A. RUBBINO and C. TRIBUNO: 
Nuovo Cimento, 7, 729 (1958). 

(8) B. C. DIVEN and G. M. ALMAY: 
Rev., 80, 407 (1950). 
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already observed in other elements by 
various authors (*!!), which must be 
interpreted, as is usual with a direct 
emission process, the Au tail is somewhat 
smaller than observed by us (! in Cr 
and Ta (Fig. 2). In order to get a rough 


N 
COR ow Ml I 
200; 
100 
0 
0 5 En (MeV) 10 19 
Fig. 2. — Photoneutrons spectra from Au, Cr 


and Ta for En > 1 MeV. N= number of photo- 

neutrons per 0.5 MeV interval. Arbitrary units. 

Figures are normalized to the same total number 
of neutrons. 


idea of the relative importance of the 
direct process we considered the very 
schematic assumptions that ~15% of 
the photoneutrons are emitted through 
direct processes and that the residual 
nucleus is left in its ground state. 

Under these assumptions and taking 
the experimental photonuclear cross 
section (°), we get the curve (b) in Fig. 1. 

The experimental points for #,>-2 MeV 
agree well with curves (a)+(b). 

We point out that under the same 
the theoretical spectrum 
from Pb agrees well with the experi- 
mental one obtained by Tous 
STEPHENS (1). 


assumptions 


and 


(©) P.UR- BYERLY jr. and W. E. 
Phys: Rev. 81, 473 (1951). 

@°) G. A. PRICE: Phys. Rev., 93, 1279 (1954). 

(1) M. i. ‘Toms and W. E: STEPHENS: 
Phys. Rev., 108; 77 (1957). 
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Turning now to the low energy part 
of the spectrum, we find that the low 
energy peak in the Au spectrum is 
obtained for neutron energy higher than 
that expected according to the evapo- 
ration theory. 

The position of the calculated max- 
imum may suffer some shift according 
to the assumed nuclear level density 
distribution but, for a coincidence of the 
calculated with the experimental max- 
imum one obtains a much higher dis- 
agreement in the other region of the 
spectrum. 

The resolution of the nuclear emul- 
sions in the low energy region (1--2 MeV) 
is rather poor. Besides the discrimi- 
nation of the tracks according the dip 
angle becomes very critical in the same 
region. Therefore the experimental fi- 
gures in this interval may have experi- 
mental errors much higher than the 
statistical ones. Therefore we do not 
attribute any special significance to the 
small disagreement between theoretical 
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Fig. 8. — Experimental energy spectra from 
Cu (), Pb (1) and Au, and theoretical curves 


calculated for evaporation. 
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and experimental position of low energy 
maximum. However we point out that 
a similar displacement may be observed 
in the photoneutron spectrum from 
Pb (1!) by Toms and STEPHENS (Fig. 3). 

On the contrary an agreement on 
the position of the maximum is found 
for the photoneutrons from Cu (?). 
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At altitude of 2550 m a.s.l., 
search is being carried out for the par- 
ticles of ~ 550 electronic re- 
ported by ALYHANTAN et al. (1). In about 
one year of measurements performed at 
an altitude of 3200 m a.s.l. these authors 
reported 11 examples of such a particle, 
at a rate of ~1/200 that of cosmic-ray 
muons stopping in the same range 
terval ~ 60 g/em? 
copper). 

Recently OREAR et al. (?) did not 
find any 550-mass particle in a sample 
of 550 muons stopped in an emulsion 
stack exposed to cosmic rays. An indi- 


an a 


masses 


in- 


(average range of 


cation for the existence of 550-mass 
particles would appear, nevertheless, 
from more recent results of KEUFFEL 


and coworkers (3) in their mass spectrum 


(1) A. I. ALIHANIAN ef al.: International 
Conference on High Energy Physics (Moscow, May 
1956); A. I. ALIHANIAN, N. V. SOSTAKOVIÙ, 
A. T. DADAIAN, V. N. FEDOROV and B. N. DE- 
RIAGIN: Zurn. Exp. Theor. Phyz., 4, 817 (1957), 
also for previous references. 

(2) J. OREAR, E. BIERMAN and S. ROSEN- 
DARFF: Postdeadline paper for the 1958 New 
York Meeting of the Am. Phys. Soc. 

(3) J. W. KEUFFEL, R. L. CALL and W. H. 
SANDMANN: Bull. Am. Phys. Soc., 3, 162 (1958). 


of cosmie under 15 em 
of lead (*). 

In the present experiment mass de- 
terminations are essentially based on 
measurements of energy loss and range. 
A fairly large solid angle of acceptance 
(10.5 ster em?) should make it possible to 
detect some 550-mass particles in a 
comparatively short time, if they exist 
in the abundance reported by the Soviet 
Group (1). 

The apparatus is shown in Fig. 1. 
The large bottom chamber (4) (C.R.) 
contains 17 plates of plexiglass each 
0.3 em thick and is used for range and 
scattering measurements. The change of 
the ionization in the various « gaps » of 


ray particles 


C.R. can be used to obtain further 
information. The top chamber (C.I.) 


is used essentially to recognize, from 


lack of alignment of the tracks in 
C.I. and C.R., possible interactions 


occurring between C.R. and the pro- 


(*) A few days ago other results have been 
reported, together with those contained in this 
letter, at the C.E.R.N. Conference (Geneva, 
July 1958), in the general report on strange 
particles by J. STEINBERGER. 

(*) E. FIORINI, R. GIaccoNI and C. SUCCI: 
Nuovo Cimento, 6, 943 (1957). 
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portional liquid scintillator D. The 


b) do not cross either the Geiger 


pulses of this counter, which allows 
measurements of energy loss of about 


Fig, 


15% accuracy, are displayed (after a 
convenient amplification of linearity 
better than 5% and gain stable within 
+2%) on the sweep of a «515 Tek- 
tronics » oscilloscope, whenever a trig- 
gering pulse occurs. 
Triggering pulses are produced by 
ionizing particles that: 
a) cross the solid angle defined by 
the Geiger counters A’, A”, B, 0 and 
the scintillator D: 
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counters A”, or the large liquid sein- 
tillator /; 


t vii 
\sertsinve die 
ECO 

\cloud chanbe 


| [Sensitive volume \ \ 
| oF cloud chamber \y 


sa CRI 


| t 
/ \ 
anni. 
j ] \ 
dh 
c) do not produce light crossing 
the Cerenkov counter C, (efficiency 


> 99.9% for relativistic particles). 


The cloud chambers are expanded by 
the triggering pulses and stereoscopic 
pictures are taken, simultaneously with 
the picture of the oscillogram giving the 
pulse height of the counter D. 

Mixed coincidences between counter 
B, and C,,B, and 0,, B, and C,, ensure 
that the solid angle of acceptance is 
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CONVERSI, E. 
entirely covered by the large scintil- 
lator F, used in anticoincidence. The 
top lead thickness, the anticoincidence 
counters A” and the anticoincidence 
Cerenkov counter C,, reduce to a small 
value the triggering background due to 
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the electrophotonic component of the 
cosmic radiation. 

The anticoincidence counter C, is 
such that electrons, muons and pions 
crossing counter D are all eliminated. 
On the contrary, 550-mass particles cer- 


CRE 


264 


2 


A SEARCH FOR PARTICLES OF 550 m, 


tainly do not produce light in C,, if they 
are not accompanied by fast ionizing 
particles and lose energy only by ioniza- 
tion before stopping in the plates of C.R. 

Up to date 1556 pictures have been 
collected in the nearly 610 h of measu- 
rement. This time must be reduced to 
— 430 h on account of the dead time 
(~7 min) of the chambers. We have 
analysed only those pictures showing a 


Fig. 


particle which, after entering the appa- 
ratus, stops (with or without interaction) 
between the 3 and the 16 plate of C.R. 
The rate of muons stopping in this same 
range interval after crossing the solid 
angle of the apparatus, is estimated to 
be about 3 per hour (the Cerenkov 
counter prevents light mesons from being 
actually recorded). This rate is derived 
from the well known sea level inten- 


sity (#8) and altitude dependence (9?) 
(®) B. Rossi: Rev. Mod. Phys., 20, 537 
(1948). 


(6) M. CONVERSI: Phys. Rev., 79, 749 (1950). 
(7) B. RossI, M. SANDS and R. SARD: Phys. 
Rev., 72, 120 (1948). 
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Hence, if 
550-mass particles exist in the abound- 
ance reported (!) we should have re- 
corded 6 of them. 

Each event analysed is represented 
in Fig. 2 by a dot or by a small circle. 

A dot corresponds to a particle which 
stops after having lost energy in C.R. 
only by ionization. On account of the 
small thickness of the plates of C.R., 


of slow cosmic-ray muons. 


3. 

these events can be recognized un- 
ambiguously by visual estimate of their 
change in ionization along the track. 
A small circle corresponds to a stopping 
particle which clearly interacts in C.R. 

The shaded region of Fig. 2 is that 
in which 550-mass particles should fall 
(if they lose energy only by ionization) 
with a probability of ~ 90%. This can 
be seen directly from the distribution of 
the identified protons around the cor- 
responding curve of Fig. 2. 

All but one of the 9 events 
the shaded region of Fig. 2 correspond 
to protons (or K-particles) clearly inter- 
acting in one of the plates of C.R. In 


inside 
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each of these 8 events the track of the 
triggering particle is characterized by 
lack of appreciable scattering and change 
in ionization in the gaps before the plate 
in which it ends. 

In Fig. 3a is reproduced the picture 
(taken while C.I. was unfortunately not 
operating) of the event corresponding to 
the only dot falling in the shaded region 
ot Fig. 2. For purpose of comparison 
the picture of an identified proton 
stopping in the same plate of C.R. is 
reproduced in Fig. 3b. The absence of 
an appreciable scattering through the 
plates and the small change of ionization 
in the gaps of C.R., as well as the 
similarity of the track of Fig. 3a to 
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that of Fig. 3b, strongly suggests that 
the only dot falling in the shaded region 
of Fig. 2 corresponds to a proton which 
lose part of its energy interacting in 
the material between counter D and C.R. 

If this interpretation is correct we 
must conclude that no 550-mass particles 
have been recorded in a time interval in 
which we expected 6 of them, on the 
basis of their reported abundance (1). 


* OK OK 
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